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Abstract 

The present knowledge of QCD confining forces between static test charges is 
summarised, with an emphasis on lattice results. Recent developments in relating 
QCD potentials to quarkonium properties by use of effective field theory methods 
are presented. The validity of non-relativistic QCD and the adiabatic approxi- 
mation with respect to heavy quark bound states is explored. Besides the static 
potential and relativistic correction terms, the spectra of glueballs and gluinoballs, 
hybrid excitations of the QCD flux tube between fundamental colour sources, po- 
tentials between charges in various representations of the SU (3) gauge group, and 
multi-particle interactions are discussed. Some implications for quarkonia systems 
and quark-gluon hybrid mesons are drawn. 
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1 Motivation 



The phenomenology of strong interactions contains three fundamental ingredients: the 
confinement of colour charges, chiral symmetry breaking and asymptotic freedom. The 
latter requirement culminated in the invention of quantum chromodynamics (QCD) 
some 25 years ago. Predicting low energy properties of strongly interacting matter still 
represents a serious theoretical challenge. This is particularly disappointing since non- 
perturbative techniques are not only important in QCD but also for an understanding of 
physics beyond the standard model or perturbation theory. For instance a rigorous proof 
is still lacking that shows QCD as the microscopic theory of strong interactions to give 
rise to the macroscopic properties of chiral symmetry breaking and quark confinement. 

So far Lattice Gauge Theory |IJ constitutes the only known entirely non-perturbative 
regularisation scheme. By numerically simulating gauge theories on a lattice, one can 
in principle predict properties of interacting QCD matter without any non-QCD input 
(except for the quark masses) . Such simulations have provided convincing evidence not 
only for quark confinement but also for chiral symmetry breaking. Moreover, at finite 
temperature, pure gauge theories are found to undergo a confinement-deconfinement 
phase transition |3], [|, while chiral symmetry is restored at high temperature [§, , in 
QCD with sea quarks. The accuracy of these results has been tremendously improved 
during the past decade with the availability of more powerful computers and advanced 
numerical techniques. 

Unfortunately, the speed and memory of present day computers still allows only for 
"solving" relatively simple QCD problems to a satisfactory precision. One particular 
weakness that the standard lattice methodology shares with, for instance, the QCD sum 
rule approach || ||] is the difficulty in calculating properties of radially excited hadrons. 
In simple potential models, however, the spectrum of such excitations can easily be 
computed. Such models have been successfully applied in quarkonium physics since the 



discovery of the J/ip resonance more than two decades ago [10, 11, 0, |13[ 14, 15, 16, 

A Hamiltonian representation in terms of functions of simple dynamical variables 
such as distance, angular momentum, relative momentum and spin allows for an un- 
derstanding of the underlying system that is rather transparent and intuitive. One 
would like to clarify what component of the success of this simple picture results from 
the freedom of choice in constructing a phenomenological Hamiltonian and what part 
indeed reflects fundamental properties of the underlying bound state dynamics. Not 
long ago, a semi-relativistic Hamiltonian that governs heavy quarkonia bound states 



has been directly derived from QCD |22], ||, ||, ||, ||, |7|, |8|, g|, B H • Starting 
from a non-relativistic expansion of the QCD Lagrangian (NRQCD) P3L |33], |3"5| , the 
gluonic degrees of freedom have been separated from the heavy quark dynamics into 
functions of the canonical coordinates (the potentials) and integrated out by means of 
lattice simulations []29| . The resulting Hamiltonian incorporates many properties of the 
previously proposed purely phenomenological or QCD inspired models. 

Heavy quarks closely resemble static test charges which can be used to probe mi- 
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croscopic properties of the QCD vacuum, in particular the anatomy of the confinement 
mechanism. Indeed, from charmonium spectroscopy and even more so from bottomonia 
states, a lot has been learned about the nature and properties of QCD confining forces. 
Either motivated by experimental input or by QCD itself, many effective models of low 



energy aspects have been proposed, in particular bag models [Bq, 37, pq, KM HQ], strong 



coupling and flux tube models Ell H3, H3l H4|, bosonic string models [H3, 46, 47], the 



stochastic vacuum model [48], [49], [!>QJ, dual QCD |5l], [52], [53|, [54| and the Abelian Higgs 



model |55[, instanton based models ][5(], [57], |58|] and relativistic quark models |59|. Many 
of these models are either expected to apply best to a non-relativistic setting or can 
most easily be solved in the situation of slowly moving colour charges. 

In view of the fact that many problems like properties of complex nuclei are unlikely 
ever to be solved from first principles alone, to some extent modelling and approxima- 
tions will always be required. Recently, using the stochastic vacuum model as well as 
dual QCD and the minimal area law, that is common to the strong coupling limit and 
string pictures, the potentials within the quarkonium bound state Hamiltonian have 
been computed |60|, |6lJ , and compared to lattice results to test the underlying assump- 
tions in the non-relativistic setting |)2|, It is a challenge for lattice simulations to 
realise simple QCD situations in which low energy models can be thoroughly checked. 

Predictions of low energy quantities like hadron masses and form factors are the 
obvious phenomenological application of lattice QCD methods. In view of the new b 
physics experiments Babar, Belle, HERA-B and LHCb, precise non-perturbative QCD 
contributions to weak decay constants are required to relate experimental input to the 
least well determined CKM matrix elements. Heavy-light systems are also thought to 
be sensitive towards CP violations. In view of the proposed linear electron colliders 
NLC and TESLA a calculation of the top production rate, e + e~ — > ti, near threshold is 
required to precisely determine the top quark mass and even in this high energy regime 
non-perturbative effects might turn out to play an substantial role. Therefore, develop- 
ing heavy quark methods and verifying their accurateness against precision experimental 
data from quarkonium systems is of utmost interest. Even quarkonia themselves con- 
tain valuable information. For instance, one would expect cleaner discriminatory signals 
for heavy quark-gluon hybrid states, that should exist as a consequence of QCD, than 
for their light hybrid counterparts. Moreover, the first B c mesons have recently been 
discovered and it is a challenge to predict their spectrum. Last but not least, quarko- 
nia systems contain information on the c and b quark masses that are fundamental 
parameters of the Standard Model. 

This report is organised as follows: in Section 0, phenomenological evidence for 
linear confinement from the spectrum of light mesons and quarkonia is presented. In 
Section ||, a brief introduction to the lattice methodology is provided before the present 
knowledge on the static QCD potential will be reviewed in Section In view of latest 
results from lattice simulations including sea quarks, particular emphasis is put on the 
"breaking" of the hadronic string in full QCD. Subsequently, in Section |5| static forces 
in more complicated situations, in particular hybrid potentials, bound states involving 
static gluinos, potentials between charges in higher representations of the SU(N) colour 
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group, and multi-body forces are discussed. In Section [| attention is paid to relativistic 
corrections to the static potential and the applicability of the adiabatic approximation. 
The results are then applied to quarkonium systems in Section [7]. 



2 The hadron spectrum 

The discovery of asymptotically free constituents of hadronic matter in deep inelastic 
scattering experiments gave birth to QCD as the generally accepted theory of strong in- 
teractions. However, the most precise experimental data to-date, the hadron spectrum, 
have been obtained in the low energy region and not at the high energies necessary to 
resolve the quark-gluon sub-structure of hadrons. While perturbative QCD (pQCD) 
should be applicable to high energy scattering problems to some extent, solving QCD 
in the low energy region poses a serious problem to theorists: not only does one have to 
deal with a strongly coupled system but also with a relativistic many-body bound state 
problem. Moreover, unlike in the prototype gauge theory, QED, even on the classical 
level the QCD vacuum structure is non-trivial, giving rise to instanton induced effects 
for example. 

It is instructive to consider the historical developments that culminated in the dis- 
covery of QCD, in particular since the pre-QCD era was dominated by concepts that 
were almost exclusively inspired by non-perturbative phenomenology, such as the reso- 



nance spectrum. General S'-matrix properties and dispersive relations |64], [)5| formed 
the formal basis of such pre-QCD developments. A serious conceptual problem of the S- 
matrix approach (also known as the bootstrap) is the fact that the unitarity of tree level 
scattering amplitudes is broken as soon as one allows for virtual point-like quanta of spin 
larger than one to be exchanged between external particles. This observation was one of 
the motivations for Veneziano's duality conjecture |j6| and the dual resonance model of 



the late 60s which finally culminated in the invention of string theories |67, |68], [E| |69| . 

While the ^-matrix framework addressed dynamical issues of strong interactions, the 
naive SUp(3) quark model |7D|, [7T| served well in classifying all known hadronic states, 



in particular after it had been extended by the colour SU(3) degrees of freedom |72|, [73[ 



However, the quark model alone did not relate to any dynamical questions of the un- 
derlying interaction. For instance, no explanation was provided for the alignment of 
particles of mass m and spin J along almost linear Regge trajectories in the m? — J 
plane 173], |63|. Bosonic string theories finally did not only resolve the unitarity puz- 
zle of the S'-matrix theory but also offered an explanation for the linearity of Regge 
trajectories pR ER [69j. However, string theories encountered internal inconsistencies 



when formulated in four space-time dimensions [7J| and were also incompatible with 



the Bj0rken scaling observed in e p collisions |f7S |. An explanation for the latter was 



provided by the invention of partons [77], |78], [79] and asymptotic freedom. 

With the advent of QCD dynamics |5(], BTj], these partons were identified as the 
quarks of the eightfold way and became the accepted elementary constituents of hadronic 
matter: the string theory of strong interactions that had been developed in parallel sur- 
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vived only as a possible low energy effective theory, in four space-time dimensions. While 
QCD — unlike all preceeding suggestions — certainly explains asymptotic freedom, it 
is still unproven that it indeed results in collective phenomena such as the confinement 
of quarks and gluons or chiral symmetry breaking. However, lattice simulations provide 
convincing evidence. 

It is legitimate to speculate whether QCD really contains all low energy information: 
is the set of fundamental parameters that describes the hadron spectrum compatible 
with the parameters needed to explain high energy scattering experiments or is there 
place for new physics? For example a (hypothetical) gluino with mass of a few GeV 
would affect the running of the QCD coupling between mz and typical hadronic scales 
that are smaller by two orders of magnitude. Is QCD the right theory at all? If so, 
quark-gluon hybrids and glueballs should show up in the particle spectrum. Although 
these general questions are not central to this article they motivate continued phe- 
nomenological interest in QCD itself from a general perspective. 

The discovery of states composed of heavy quarks, namely charmonia in 1974 and 
bottomonia in 1977, enabled aspects of strong interaction dynamics to be probed in a 
non-relativistic setting. By means of simple potential models a wealth of data on energy 
levels and decay rates could be explained. The question arises: if these models yield 
the right particle spectrum, can they eventually be derived from QCD? What do such 
models tell us about QCD and what does QCD tell us about such models? 

Before addressing these questions in later Sections, here some aspects of hadron 
spectroscopy that relate to flux tube and potential models are summarised. 



2.1 Regge trajectories 



Since the early sixties it has been noticed that mesons as well as baryons of mass m 
and spin J group themselves into almost linear, so-called Regge trajectories [|4], |65], [74[ 
in the m 2 — J plane up to spins as high as J = 11/2. In Table ^TT] the light meson 
spectrum is summarised. Only resonances that are confirmed in the Review of Particle 
Properties [32 have been included. The it, K*, K% and K triplets have been replaced 



by their weighted mass averages. The second column of the Table represents the J 
assignment. Each increase of the orbital angular momentum by one unit results in a 
switch of both, parity and charge assignments. 



The data of Table |2.1| is displayed in Figure 2T], together with linear fits of the form 



J(m) = a(0) + a'm z 



(2.1) 



Similar plots can be made for the baryon spectrum. a(0) is known as the Regge intersect 
and, 

(2.2) 



a 



27TCT' 



as the Regge slope. The resulting values for the "string tension", a, are displayed in 



Table |2.2| . While statistical errors on the data points increase with J, the applicability of 
the relativistic string model that, as we shall see below, predicts the linear dependence 



6 



Table 2.1: Light meson masses. 



state 


JP(C) 


m/MeV 


71 

6i(1235) 
vr 2 (1670) 


o-+ 

1+- 

2-+ 


138 
1229(3) 
1670(20) 


p(770) 
a 2 (1320) 
p 3 (1690) 
a 4 (2040) 


1— 

2++ 

3~ 
4 ++ 


770(1) 
1318(1) 
1691(5) 
2020(16) 


a;(782) 
/ 2 (1270) 
cj 3 (1670) 
/ 4 (2050) 


1" 

2++ 

3- 
4++ 


782 
1275(1) 
1667(4) 
2044(11) 


0(1020) 
£(1525) 
03(1850) 


1~ 

2++ 
3~ 


1019 

1525(5) 
1854(7) 


V 

^(1170) 


o-+ 

1+- 


547 
1170(20) 


K 

#i(1270) 
X 2 (1770) 


0" 
1+ 
2" 


495 
1273(7) 
1773(8) 


AT*(892) 
fT 2 *(1430) 
X*(1780) 
fT 4 *(2045) 


1- 

2+ 
3" 
4+ 


893 
1428(2) 
1776(7) 
2045(9) 



Table 2.2: String tensions from Regge trajectories. 



trajectory 


v^/MeV 


AJ 


7T,6i,... 


469(6) 


0.06 


p,a 2 , ■ ■ ■ 


429(2) 


0.03 


u, / 2 , . . . 


436(8) 


0.12 




437(5) 


0.06 


K,K U ... 


480(4) 


0.04 


K*,K*,... 


424(5) 


0.07 
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Figure 2.1: Regge trajectories. 

is expected to improve with J. Therefore, in the fits we have decided to ignore the 
experimental errors and give all points equal weight. AJ denotes the root mean square 
deviation between fitted angular momenta and data points, normalised by the root of 
the degrees of freedom (i.e. the number of data points minus two) and reflects the overall 
quality of a fit. 




Figure 2.2: Rotating quarks, connected by a string of constant energy density 

A simple explanation of the linear behaviour is provided by the relativistic string 
model [^, |45| : imagine a rotating string of length 2d with a constant energy density 
per unit length, a (Figure |2.2|) . If this string spans between (approximately) massless 
quarks, we might expect those quarks to move at (almost) the speed of light, c = 1, 
with respect to the centre of mass. The velocity as a function of the distance from the 
centre of the string, r, in this set-up is given by, v(r) = r/d. From this, we calculate 
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the energy stored in the rotating string, 



rd (£f (j 

m = 2 . = nda, (2.3) 



and angular momentum, 



J = 2 [ d d ; arv{r) = ^ = (2.4) 
Jo i - v *( r ) 2 2na 



which results in the relation of Eq. (|2.2| ) between Regge slope, a', and string tension, 
a. This crude approximation can of course be improved. For example, one can allow 
for a rest mass of the quarks. Velocities smaller than c will result in a slight increase 
of the Regge slope. The assumption that the string energy entirely consists of a lon- 
gitudinal electric component in the co-rotating frame yields predictions for spin-orbit 



splittings ||2[ etc.. 

For the two Regge trajectories starting with a pseudo-scalar (jr and K), one finds 
values, 470 MeV < ^fo < 480 MeV, while all other numbers scatter between 424 and 
437 MeV. The value extracted from the p,a,2, ■ ■ • trajectory, which is the most linear 
one, is y/a = (429 ± 2) MeV. 

2.2 Quarkonia 

Soon after the discovery of the J / ip meson in e + e~ annihilation, the possibility of a non- 
relativistic treatment of such states, in analogy to the positronium of electrodynamics, 



was suggested JTTJ] . Quarkonia, i.e. mesonic states that contain two heavy constituent 
quarks, either charm or bottomF], owe their name to this analogy. Within the quark 
model, the quark anti-quark system can be characterised by its total spin, S = Si + S2 
(s = or s = 1), the relative orbital angular momentum, L, and the total spin, 
J = L + S. Within the standard spectroscopic notation, n 2s+1 lj, n denotes the radial 
excitation while I = is labelled by the letter S, I = 1 by P, I = 2 by D etc.. The 
parity of a quark anti-quark state is given by, P = (— while the charge conjugation 
operator (if quark and anti-quark share the same flavour) has eigenvalue, C = (— l)' +s . 

In making the above J PC assignments, we ignore the possibility of the gluonic de- 
grees of freedom contributing to the quantum numbers. This simplification results in 
certain combinations to be quark model forbidden (or spin-exotic), namely, J PC = 
H , , 1 h , 2 + ~, 3~ + , . . .. Another aspect is that some 

jPC 

assignments can be gen- 
erated in various ways. For instance, 3 Si and 3 Di states both result in J PC = 1 . 
As soon as gluons are introduced, the relative angular momentum, L, is not conserved 
anymore and physical vector particles will in general be superpositions of excitations 
from these two channels: strictly speaking, only the number of nodes of the wave func- 
tion, n, the spin J, parity P, charge C (in the case of flavour singlet mesons), and the 

: Due to the large weak decay rate, t — > bW + , the top quark does not appear as a constituent in 
bound states (see e.g. Ref. Q). 
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Table 2.3: Classification of charmonium and bottomonium states. 



n 2s+1 lj 


jPC 


cc 


bb 


1% 


o-+ 


Vc 


Vb 


1 3 S 1 


1— 


J/i> 


T 


2 3 S 1 


1— 


iP(2S) 


T(2S) 


l'Pi 


1+- 


h c 


h 


l 3 Po 


0++ 


XcO 


Xbo 


l 3 Px 


1++ 


Xcl 


Xbi 


1 3 P 2 


2++ 


Xc2 


Xbi 



constituent quark content (neglecting annihilation processes and weak decays) represent 
"good" quantum numbers. 

In Table JO, we have compiled quantum numbers and names for some members of 
the J /if) and T families. Little is known experimentally about B c mesons, which are 
bound states of a b and a c quark. For these particles an additional peculiarity has to be 
considered: charge and total spin are no longer "good" quark model quantum numbers. 
For I > 1 this results in mixing between the J = I would-be singlet and would-be triplet 
states. 



1200 


ii ii 


1000 


BB 
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400 




200 









-200 


ii ii 



Y Tl c >J/¥ %b0,1,2 Xc0,1,2 

Figure 2.3: Energy splittings, AE, for charmonia and bottomonia with respect to the 
T and J /if) triplet S wave states. 

In Figure |2.3| , all experimentally determined splittings with respect to the l 3 ^ state 
for the T and J /if) families are depicted. We have restricted ourselves to states, listed 
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in the Review of Particle Properties |82| , that are below the DD and BB thresholds 
(dashed horizontal lines) for charmonia and bottomonia, respectively, with the exception 
of the T(4S). While the mass of the J/ip (3.097 GeV) considerably differs from that of 
the T (9.46 GeV), indicating a substantial difference in the quark masses, m& ~ 3m c , 
both 2 3 Si — l 3 Si splittings agree within 5 % (589 and 563 MeV). We define the spin 
averaged x mass by, 

rrisp = ~ (ma Po + 3ms Pl + 5m 3 p 2 ) mi Pl . (2.5) 
9 

Again, within a few per cent, the 1 3 P — l 3 Si splittings agree (429 MeV vs. 440 MeV). 
Unfortunately, while the r\ c has been discovered, no pseudo-scalar bb 1 S , meson has 
yet been seen, such that a consistent comparison with respect to spin averaged S state 
masses, 

m s = ^ ( mi s + 3m 35l ) , (2.6) 

is not possible. 

While the 2S — IS and IP — IS splittings seem to agree within a few per cent, the 
fine structure splittings between the P states come out to be almost three times as large 
in the charm case, compared to that for the bottom, 

This is consistent with the expectation that in the limit of infinite quark mass, fine 
structure splittings will eventually completely disappear, in analogy to hydrogen-like 
systems. However, for the ratio between the respective m X2 — m Xl splittings one finds 
a different numerical value, 0.47(2), indicating a more complicated dependence on the 
inverse quark mass than mere proportionality. 

For sufficiently heavy quarks, one might hope that the characteristic time scale 
associated with the relative movement of the constituent quarks is much larger than 



that associated with the gluonic (or sea quark) degrees of freedom [11]. In this case 
the adiabatic (or Born-Oppenheimer) approximation applies and the effect of gluons 
and sea quarks can be represented by an averaged instantaneous interaction potential 
between the heavy quark sources. Moreover, the bound state problem will essentially 
become non-relativistic and the dynamics will, to first approximation, be controlled by 
the Schrodinger equation, 



P 2 

f- + V{r) 



ip n ii 3 (r) = E nl ip nU3 (r), (2.8) 



with a potential, V(r) (r = |r|), or, if spin effects are taken into account, semi- 
relativistic Pauli-Thomas-like extensions. In the adiabatic approximation quarkonia 
are the positronium of QCD. However, unlike in QED where the interaction potential 
can be calculated perturbatively and the spectrum predicted, we are faced with the 



11 



Table 2.4: The scaling of level splittings, AE, inverse length scales, (r x ), and the 
squared velocity, (v 2 ), with reduced mass and coupling. 



potential 


AE 


(r- 1 ) 


(v 2 ) 


Coulomb, Eq. Q 






e 2 


logarithmic, Eq. (|2.10|) 






<W 


linear, Eq. (|2.11) 






- 2/3 ^ 4/3 



inverse problem of determining or guessing the interaction potential and the reduced 
quark mass, = m/2, from the observed spectrum, E n i, and decay rates. The latter 
can be related to properties of the wave function at the origin |Tj|. If the adiabatic 
approximation is justified we would expect, to leading order in a semi-relativistic ex- 
pansion, the same potential to explain cc as well as bb spectra since QCD interactions 
are flavour blind. On the other hand it is clear that the adiabatic approximation will 
at least fail for unstable excitations like the ^(SS') or T(4S) since decays cannot be 
accounted for by a one channel Hamiltonian with a real potential. 

In Appendix we derive general properties of the spectrum for power law and 
logarithmic potentials. The main results for a Coulomb potential, 

V(r) = --, (2.9) 
r 

a logarithmic potential, 

V(r) =Cln(—), (2.10) 
\r J 

and a linear potential, 

V(r) = ar, (2.11) 

are displayed in Table |2~4| . 

From the spin-averaged quarkonia spectra it is evident that the underlying potential 
cannot be purely Coulomb type. Otherwise, the 25* — 15 splitting would be approxi- 
mately degenerate with the lowest lying nP — IS splitting and, moreover, T splittings 
would be enhanced with respect to J/ip splittings by the ratio of the quark masses, 
mb/m c pa 3. However, a logarithmic potential that would explain the approximate mass 
independence of spin-averaged splittings is incompatible with tree level perturbation 
theory, i.e. Eq. ( ^9"|) , with e = (4/3)a s . 

The Cornell potential fI2fl , 

V(r) = --+<rr, (2.12) 
r 

contains the perturbative expectation plus an additional linear term. The parameters e 
and a can be adjusted such that within the range of charm and bottom quark masses, the 
linear dependence of the Rydberg energy on fi R is compensated by the 1 / fi]( 3 behaviour 
expected from the large distance linear term: within the distance scales relevant for the 
quarkonium bound state problem, the Cornell potential looks effectively logarithmic. 
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At quark masses larger than nib, the Coulomb term will eventually dominate and split- 
tings will diverge in proportion with Note that the Cornell potential predicts the 
average velocity, (v 2 ) oc AE/ fi R , to saturate at the value (v 2 ) = e 2 for large quark mass 
while from the approximate equality of bottomonia and charmonia level splittings one 
would expect {v 2 ) / (v 2 ) ~ m c /m b . (v 2 ) quantifies the quality of the non-relativistic ap- 
proximation while the applicability of the adiabatic approximation is more complicated 
to establish from a QCD perspective. 

Before the discovery of the T(2S), fits to the spin averaged quarkonia spectra re- 
sulted in parameter values |R|, e ~ 0.25 and yfa ~ 455 MeV. After inclusion of the 



T states, that probe the potential at smaller distances, values like |p.6|, e ~ 0.51 and 



y/a « 412 MeV, and jT7|, e & 0.52 and yfa « 427 MeV, emerged. However, within 
the region, 0.2 fm < r < 1 fm, which is effectively probed by spin-averaged quarko- 
nia splittings, the e ~ 0.25 parametrisation only marginally differs from the e ~ 0.5 
parametrisations; the higher value of the Coulomb coefficient is compensated for by a 
smaller slope, a. Interestingly, the slope of the Cornell potential is in qualitative agree- 
ment with yfa ~ 430 MeV, the estimate of the string tension from Regge trajectories 
of light mesons, discussed in Section 

While the spin averaged spectrum probes the potential at distances r > 0.2 fm, 
fine structure splittings are sensitive towards the Lorentz and spin structure of the 
interacting force as well as to the functional form of the potential at short distances. 
We shall discuss this in detail in Section [7|. 



3 Lattice methods 

Lattice QCD was invented by Wilson [|lj shortly after QCD emerged as the prime 
candidate for a consistent theory of strong interactions. The main intention was to define 
an entirely non-perturbative regularisation scheme for QCD, based on the principle of 
local gauge invariance. Besides regulating the theory, the lattice lends itself to strong 
coupling expansion techniques in terms of the inverse QCD coupling, 

_ 2N _ 2N ^ 
g 2 47ra s ' 

Such techniques complement the conventional perturbative weak coupling expansion 
and have, in particular in the Hamiltonian formulation of lattice QCD H84| , stimulated 
the flux tube model of Ref. |[44||. However, so far nobody has managed to analytically 



relate the strong coupling limit of QCD to weak coupling results. For instance, in U(l) 
as well as in SU (N) gauge theories one obtains an area law for Wilson loops [Eq. ((4.1|)1, 
i.e. confinement, in the strong coupling limit. While in (3 + l)-dimensional U(l) lattice 
gauge theory the strong coupling regime is separated from a non-confining weak coupling 
region by a phase transition, in SU (N) one would hope that no such phase transition 
at finite (3 exists and confinement survives at weak coupling. 

Besides offering new analytical insight and techniques, the lattice approach to QCD 



lends itself to treatment on a computer [85, pq, m. To allow for a numerical evaluation 



13 



of expectation values it is convenient to work in Euclidean space-time in which a path 
integral measure can be defined. Moreover, the time evolution operator becomes anti- 
Hermitian which results in n-point correlation functions decaying exponentially, rather 
than exhibiting oscillatory behaviour. Most results that are obtained in Euclidean space 
can be related to the space-like region of the Minkowski world and can in principle be 
analytically continued into the time-like region of interest. With results that have been 
obtained on a discrete set of points with finite precision, however, such a continuation 
is anything but straight forward. Fortunately, unless one is interested in real time 
processes like particle scattering, this is in general not required. In particular the mass 
spectrum remains unaffected by the rotation to imaginary time as long as reflection 



positivity holds |p7[ , which is the case at least for the lattice actions discussed in this 
article. 

In what follows, the aspects of lattice simulations that are relevant for our discussion 
are summarised. For a more detailed introduction to Lattice Gauge Theories the reader 



may consult several books [§8], |89|, |90], and review articles |92|, p6|, |93|, |94|, p5|, |96| . 
The conventions that are adapted throughout the article are detailed in Appendix ||. 

3.1 What can the lattice do? 

The lattice allows for a first principles numerical evaluation of expectation values of 
a given quantum field theory that is defined by an action 5" in Euclidean space-time. 
However, the accessible lattice volumes and resolutions are limited by the available 
(finite) computer performance and memory. 

The obvious strength of lattice methods are hadron mass predictions. Only recently 
computers have become powerful enough to allow for a determination of the infinite 
volume light hadron spectrum in the continuum limit in the quenched approximation^ 
to QCD within uncertainties of a few per cent |)7|. To this accuracy the quenched 



spectrum has been found to differ from experiment. Some collaborations have started to 
systematically explore QCD with two flavours of light sea quarks and the first precision 
results indeed indicate deviations from the quenched approximation in the direction of 
the experimental values [Q. Even if one is unimpressed by post-dictions of hadron 



masses that have been known with high precision for decades such simulations allow 
fundamental standard model parameters to be fixed from low energy input data, like 
quark masses |99|, |100| , |101| , |102|1 and the QCD running coupling |[L03I - Of course, as we 



shall see, a wealth of other applications of phenomenological importance exists. 

Unfortunately, only the lowest radial excitations of a hadronic state are accessible 
in practice. Lattice predictions are restricted to rather simple systems too. Even the 
deuteron is beyond the reach of present day super-computers. Therefore, it is desirable 
to supplement lattice simulations by analytical methods. The computer alone acts as 
a black box. In order to understand and interprete the output values and to predict 



2 In the quenched approximation, vacuum polarisation effects due to sea quarks are neglected by 
replacing the fermionic part of the action by a constant. In the language of perturbative QCD this 
amounts to neglecting quark loops. 
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their dependence on the input parameters, some modelling is required. Vice versa, the 
lattice itself is a strong tool to validate models and approximations. Unlike in the "real" 
world, one can vary the quark masses, m$, the number of colours, N, the number of 
flavours, n^, the temperature, the spatial volume, the space-time dimension and even 
the boundary conditions in order to expose models to thorough tests in many situations. 



3.2 The method 

In a lattice simulation, Euclidean space-time is discretised on a torus| with L\L T lattice 
points or sites, x = na, rii = 0, 1, . . . , L a — 1, n 4 = 0, 1, . . . , L T — 1, separated by the 
lattice spacing^], a, that provides an ultra-violet cut-off on the gluon momenta, q < ir/a, 
and regulates the theory. Two adjacent points are connected by an oriented bond or 
link, While Dirac quark fields, q l x , are represented by 4 x N tuplesQ at lattice 

sites, x, gauge fields, 



U 



x,ii 



V 



cxp 




x+afj, 



dx„ A^{x ) 



G SU{N), 



(3.2) 



are "link variables". V denotes path ordering of the argument and (i is a unit vector 
pointing into /i direction. We further define, U x -^ = . The transformation prop- 



erty of a lattice fermion field under gauge transformations, Vt x e SU(N), is [Eq. ([B.ll[) 1 



q(x) — > {l(x)q(x), q(x) — > g(x)f2^(x) 



(3.3) 



From Eq. ( |rTTU|) 



one can derive the the transformation property of links, 



(3.4) 



U 



X,/l 



U n = Q II - 



(3.5) 



It is easy to see that the trace of a product of links along a closed loop is gauge 
invariant. Other gauge invariant objects are N gauge transporters whose colour indices 
are contracted by completely antisymmetric tensors of rank iV at a common start and 
a common end point, a quark and an anti-quark field that are connected by a gauge 
transporter or a state of N quarks whose colours are transported to a common point, 



where they are anti-symmetrically contracted. The situation is depicted in Figure [TT 
for N = 3. 

The simplest non-trivial gauge invariant object that can be constructed is the prod- 
uct of four links, enclosing an elementary square, 



TT TT . TT^ TT^ 

u x,n u x+afi,u v x +av,fi u x,v ' 



(3.6) 



3 For fermions anti-periodic boundary conditions are chosen in the temporal direction. 
4 For simplicity, we assume a\ = a-i = = = a, L\ = = L3 = L a . 

5 The superscript, i = 1, . . . n/, runs over the flavours. The factor "4" is due to the Dirac components. 
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Figure 3.1: Examples of gauge invariant objects. Lines correspond to gauge trans- 
porters. 




A A 

x + |^ + v 



A 



X + \i 



Figure 3.2: The plaquette, U XIXV . 



the "plaquette" (Figure |3T2|) . The plaquette determines the local curvature of the gauge 
fields within the group manifold, i.e. it is related to the field strength tensor, 



U. 



exp (ia 2 ^ 



X,flU I 1 



(3.7) 
(3.8) 



where we denote the normalised trace of an element in a Z)-dimensional representation 
of the gauge group by Trp or Tr, 



1 1 D 

Tr l^, = Tr^l^ = — tr 1 D = — ^2 8u = 1. 



(3.9) 



For the fundamental representation above, we have D = N. a, (3 = 1, . . . , iV label 
the colours and = F^ u T a , where the N x N matrices T a denote the gauge group 
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generators in the fundamental representation. Note that = and that = —J 7 ^ 
is anti-hermitian, as a consequence of U^ u = = Ul . 

Discretised lattice actions are formulated in a manifestly gauge-invariant way and 
should approach the continuum action in the limit, a —>■ 0. Since the action depends on 
couplings rather than directly on the lattice spacing, it is not a priori clear if this limit 
can be realised. We shall discuss the approach to the continuum limit below. For the 
moment, we remark that from the asymptotic freedom of perturbative QCD we expect 
a to approach zero as g — > 0, i.e. /3 — > oo. 

The simplest gluonic action is the so-called Wilson action, 



Sw[U]=P E [l-ReTr(^)] 



(3.10) 



X,fl>U 



where Tr denotes the normalised trace of Eq. ( |3.9|) . From Eqs. QB.14 ) and (|3.7|) it is 
easy to see that Sw = SVm[1 + 0(a 2 )]. The constant term in the action is irrelevant 
as it cancels from expectation values. The choice of the action is far from unique. For 
instance an alternative form, suggested by Manton | |104|| , has been used in the glueball 
studies of Refs. [ |105| , |106|| . The action can in principle be systematically improved to 

This Symanzik 
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approximate the continuum action to a higher order in a ||107| , 
improvement programme has first been applied to Yang-Mills lattice gauge theory by 
Liischer and Weisz 



109, 110, 111 



In a classical theory, all the coefficients of higher 
dimensional operators that are added to the plaquette of the Wilson action can easily 
be determined. However, in the quantum field theory case of interest, the coefficients 
are subject to radiative corrections, and have to be determined non-perturbatively to 
fully eliminate the 0(a 2 ) lattice artefacts of the Wilson action. Although this has not 



been achieved yet, impressive results on static potentials ||112 



TT4], the glueball 



spectrum [|115| , |116| , |117|| and thermodynamics ||118|| have recently been obtained with 
Symanzik improved gluonic actions with coefficients, approximated by a mean field 
("tadpole") estimate ||119| , |120| . An alternative improved gluonic action that has been 



used in recent lattice studies [|121| is the renormalisation group improved Iwasaki ac- 
tion ||122j, 123||. The renormalisation group approach towards an improved continuum 



limit behaviour has been systematised in the work of Hasenfratz and Niedermayer [|124j| 
on "perfect" lattice actions. Approximately perfect actions have been constructed for 
example in Refs. |25[ [L2|, |I27| . 



A naive discretisation of the Dirac fermionic action of Eq. ( |B.12| ) suffers under the 
fermion doubling problem (cf. Refs. [BR 91]). The simplest way to remove the un- 



wanted modes is to give them extra mass by adding an irrelevant term, —aqD^D^q, to 
the action. This results in Wilson fermions ||128 |, 



s f\U> Q^]=J2 qxM xy (U)q y , 



x.y 



where 



M xy = S xy -Kj2 [(1 - In) U Xtlt 5 x+ h V + (1 + 7/J Ul_^5 x ^ y 

II. 



(3.11) 



(3.12) 
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One of the disadvantages of this solution is that continuum fermions are only approxi- 
mated up to 0(a) lattice artefacts. Remember that the gauge action was correct up to 
0(a 2 ) errors. The parameter, k, is related to the inverse bare quark mass, 

2 V, Kj t\jp 



ma 



(3.13) 



where k c > 1/8 approaches the free field (U X>II = 1) limit, k c = 1/8, as (3 — > oo. Note 



that the quark fields in Eq. ( 3.11|) have been rescaled 




(3.14) 



Another popular alternative is the Kogut-Susskind action ||129|| which is correct up 
to 0(a 2 ) lattice artefacts. However, it requires four mass degenerate quark flavours. 



The Sheikoleshlami-Wohlert action [|130|1 is an O(a) Symanzik improved variant of 
the Wilson fermionic action. The coefficient of the additional term is known non- 
perturbatively |[L31|| . Other suggestions of Symanzik improved fermionic actions have 
been put forward for instance by Naik [ |132| ] and Eguchi | |133| |. Domain wall fermions 
have been suggested [|134j , 135|| , in order to realise (approximate) chiral symmetry in 
the lattice theory. These fermions have received renewed attention since they have been 



found to fulfil the Ginsparg- Wilson relation ||136| , |137|| . They share this feature with 
other fermionic actions like the "perfect" action of Ref. ||138|| and the action derived by 
use of the overlap formalism ||139| . |140|| in Ref. ||141|| . However, we are interested in quite 
the opposite of massless fermions, namely heavy quarks, such that these exciting new 
developments are of limited interest in the present context. 

Expectation values of operators, O, are determined by the computation of the path 
integral, 

(O) = \ f [dU][dq][dq] 0[U}e- s[u ' q ^. (3.15) 



The normalisation factor, or partition function, Z, is such that (1) = 1. The short- 
hand notation, q, represents {q l x } and U stands for all gauge fields, {U x ^}. The high- 
dimensional integral is evaluated by means of a (stochastic) Monte-Carlo method as an 
average over an ensemble of n representative gauge configurations^, 



, n: 



(O) 



n i=l 



(3.16) 



Therefore, the result on the expectation value is subject to a statistical error, AO, 
that will decrease like 1 / \fn: the more measurements are taken, the more precise the 
prediction becomes. For this reason one might speak of lattice measurements and lattice 
experiments, in analogy to "real" experiments. The method represents an exact approach 
in the sense that the statistical errors can in principle be made arbitrarily small by 
increasing the sample size, n. 

6 The basic numerical techniques employed to generate these configurations are e.g. explained in 
Ref. [glj and references therein. 
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3.3 Getting the physics right 

In general, the action that is simulated depends on n/ quark masses, rrii, as well as on 
a bare QCD coupling, g. By varying g and rrii the lattice spacing, a(g,mi), is changed. 
Lattice QCD is a first principles approach in that no additional parameters are intro- 
duced, apart from those that are inherent to QCD, mentioned above. In order to fit 
these rif + 1 parameters, n/ + l low energy quantities are matched to their experimental 
values: the lattice spacing, a(g, rrii), can be obtained for instance by fixing m p as deter- 
mined on the lattice to the experimental value. The lattice parameters that correspond 
to physical m u pa can then be obtained by adjusting m n /m p ; the right m s can be 
reproduced by adjusting mx/m,, or m^/m p to experiment etc.. 

If the right theory is being simulated all experimental mass ratios should be repro- 
duced in the continuum limit, a — > 0, which will be reached as g — > 0, such that it 
becomes irrelevant what set of experimental input quantities has been chosen initially. 
In practice, the available computer speed and memory are finite and simulations are 
often performed within the quenched approximation, neglecting sea quark effects, or at 
un-physically heavy quark masses. Therefore, unless controlled extrapolations to the 
right number of flavours, n/, and masses of sea quarks, rrii, are performed, residual scale 
uncertainties that depend on the choice of experimental input parameters will survive 
in the continuum limit. Once the scale and quark masses have been set, everything else 
becomes a prediction. 

Lattice results in general need to be extrapolated to the (continuum) limit, a — > 0, 
at fixed physical volume. The functional form of this extrapolation is theoretically well 
understood and under control. This claim is substantiated by the fact that simulations 
with different lattice discretisations of the continuum QCD action yield compatible 
results after the continuum extrapolation has been performed. For high energies, an 
overlap between certain quenched lattice computations and perturbative QCD has been 
confirmed too ||103| , |142|| , excluding the possibility of fixed points of the /3-function at 



finite values of the coupling, other than g = 0. After taking the continuum limit, an 
infinite volume extrapolation should be performed. In most cases, results on hadron 
masses from quenched evaluations on lattices with spatial extent, L a a > 2 fm, are 
virtually indistinguishable from the infinite volume limit within typical statistical errors 
down to pion masses, m n pa m p /3. However, for QCD with sea quarks the available 
information is not yet sufficient for definite conclusions, in particular as one might 
expect a substantial dependence of the on-set of finite size effects on the sea quark 
mass(es). The typical lattice spacings used in light hadron spectroscopy cover the 
region 0.05 fm < a < 0.2 fm. 

The effective infinite volume limit of realistically light pions cannot be realised at 
a reasonable computational cost, neither in quenched nor in full QCD. Therefore, in 
practice another extrapolation is required. This extrapolation to the physical light 
quark mass is theoretically less well under control than those to the continuum and 
infinite volume limits. The parametrisations used are in general motivated by chiral 
perturbation theory and the related theoretical uncertainties are the dominant source 
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of error in latest state-of-the-art spectrum calculations |J7|. Ideally, the Monte Carlo 
sample size n is chosen such that the statistical precision is smaller or similar in size 
than the systematic uncertainty due to the extrapolations involved. 



3.4 Mass determinations 

In order to extract the ground state mass of a state with quantum numbers a, one starts 
from a connected gauge invariant correlation function, 

C a (t) = (0|*t(f)¥ a (0)|0) - |(0|vI/ Q |0)| 2 , (3.17) 

where |0) denotes the vacuum statefl. a contains the momentum and the J PC quantum 
numbers of the state of interest as well as the constituent quark content, i.e. isospin, 
strangeness etc.. In most cases, one is interested in the rest mass. Therefore, usually 
involves a summation over all spatial positions, x, within a time slice to project onto 
spatial momentum, p = 0. Any other lattice momentum can be singled out by taking 
the corresponding discrete Fourier transform. Due to the translational invariance on the 
lattice, it is sufficient to project only either source or sink onto the desired momentum 
state. 

In what follows, we will for simplicity assume, L T a — > oo. At finite L T a additional 
contributions arise from the propagation into the negative time direction around the 
periodically closed temporal boundary. Such effects can easily be taken into account 
whenever they turn out to be numerically relevant. By inserting a complete set of 
eigenstates of the Hamiltonian, \$ a ,n), into Eq. (|3.17| ), one obtains, 

C a (t)=Y / \c n (a)\ 2 e- E ^ t . (3.18) 

n 

with 

cn(a) = ($ a ,„|* a (0)|0). (3.19) 

E n (a) is the energy eigenvalue of the state | e~ Ht \Q a ,n) = e~ £ ' n< - Q ' )t |$ Q , in ), and, 
^a(t) = e Ht ^ a (0)e~ Ht . In the limit, t — > oo, the ground state mass, 

Eo(a) = - lim j In C a (t), (3.20) 

can be extracted. The above formula converges exponentially fast and is, therefore, 
suitable for numerical studies. In general, ty a can be any linear combination of § a ,n 
and its choice is not unique. This observation is exploited in iterative smearing or 



fuzzing techniques [143, 144, |145|, 146, 147, 148, 149] that seek to prepare an initial 



state with optimised overlap to the level of interest. This will then allow the infinite 
time limit of Eq. (|3.20j) to be effectively realised at moderate temporal separations, t. 
In principle, not only a single correlation function but a whole cross-correlation matrix 

7 In Eq. ( |3.15| ) we have employed the short-hand notation, (O) = (0|O|0), for the vacuum expectation 
value of the operator O. 
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between differently optimised ^'s can be measured. In doing so, there is the chance 
that by diagonalising the matrix and employing sophisticated multi-exponential fitting 
techniques not only the ground state energy can be extracted but also those of the 



lowest one or two radial excitations 
In Eq 



151, 116, 152, 153 



\cn(a)\ 



mni 

1. The 

2 = 1, 



(|3.18| ) we have adapted the normalisation convention 
J2 n \&a,n) (&a,n\ = 1- This results in < |c n (a)| 2 < 1 and Y. 
deviation of |co(a;)| 2 , the ground state overlap, from the optimal value, |co(a)| 
determines the quality of the smeared operator, ty a . It should be noted that if 
contains Dirac spinors, e.g. if it is a pion creation operator, the standard normalisation 
condition would be, (& n ,m\&iT,n) = 2m ;r]m (5 mn , instead. As a consequence, Eq. ( |3.18| ) is 
replaced by, 



2m„ 



(3.21) 



For the manipulations yielding Eq. (|3.18|) we have assumed the existence of a positive 
definite self-adjoint Hamiltonian. Liischer [154] has shown that the Wilson gluonic 
and fermionic lattice actions fulfil both, reflection positivity [|155| , |87j with respect to 
hyperplanes going through lattice sites and through the centre of temporal lattice links 
(see also Ref. ||156|| ). This feature implies the existence of a positive transfer matrix and 
the possibility of analytical continuation to Minkowski space-time. Another important 
consequence of reflection positivity is that the coefficients of the series in Eq. ( |3.18| ), are 
non-negative and that, therefore, the limit of Eq. ( |3.2U| ) is approached monotonically 
from above. General properties of the transfer matrix for continuum limit improved 
actions are discussed in Ref. f|157| . 



3.5 The continuum limit 

A continuum limit of the lattice theory can be defined at fixed points associated to 
phase transitions of second or higher order in the space spanned by the bare couplings 
of the action. In the vicinity of such a phase transition any correlation length, £/a, 
diverges which implies, a — ► 0, if we associate £ to a physical distance or mass, £ = 
1/m. Moreover, universality sets in, i.e. the behaviour of different correlation lengths 
is governed by one and the same critical exponent. This results in ratios between two 
correlation lengths, or masses, to saturate at constant values: the system forgets the 
lattice spacing, a. One refers to this behaviour as "scaling" . In the case of the Wilson 
gluonic action, the leading order violations of scaling are expected to be proportional 
to a 2 while for the Wilson fermionic action, they are only linear in a. 
The Callan-Symanzik /3-function, 

PM = = -(3 a 2 s - ftaj - (3.22) 

parameterises the variation of the QCD coupling, a s = g 2 /(47r), with a scale /i. Per- 
turbative QCD tells us, (3$ > and f3\ > 0, which implies asymptotic freedom: the 
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limit a s = is reached with /i — > oo, i.e. the continuum limit of lattice QCD, a — > 0, 
corresponds tc0 (3 — > oo. Far away from the phase transition, no unique /3-function can 
be defined; due to the occurrence of power corrections, different masses will in general 
run differently as a function of the bare coupling. Lattice results seem to imply that in 
zero temperature SU (N) gauge theory no fixed point other than a s = exists. 
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Figure 3.3: The ratio, m p ro, in quenched QCD, extrapolated to the continuum limit. 
The p masses in lattice units are taken from Refs. |J7[] (CP-PACS) and ||158|| (GF11). 



While the coefficients /3q and fii within Eq. ( |3.22|) are universal, higher order coeffi- 



cients depend on the renormalisation scheme. Integrating Eq. ( |3.22| ) yields, 

/ r a (v) da \ . 
/x = Aexp / — — , 3.23 
\Ja(A) 2f3(a)J 

where we define the integration constant, the so-called QCD A-parameter, via the two 
loop relation, 

A=lta^e X p(-^l M )ft aM r4. (3.24) 

In Appendix |CJ, we display results on the coefficients $ of Eq. ( |3.22| ) for reference and 
detail how to translate between different schemes. 

In QCD with sea quarks, the lattice cut-off, a, will not only depend on the cou- 
pling but also on the bare quark masses of the Lagrangian. This dependence can be 
parameterised into quark mass anomalous dimension functions. The continuum limit 
of a theory with rif different quark masses will be taken along a trajectory on which rif 



Here, fi represents the inverse lattice coupling of Eq. (3.1) and not the (3 function 
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physical mass ratios are kept fixed. In the approximation to QCD with two degener- 
ate light quark masses for instance the physical curve m^jm p « 2/11 would serve this 
purpose. 
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Figure 3.4: The scalar glueball mass in units of r$ as a function of a 2 . The glueball 
masses in lattice units are taken from Refs. p| (UKQCD), |L6| (GF11), |6l], [162 
(MT) and BTB3 (FSST). 



In Figure |3.3| , we show a continuum limit extrapolation of the quantity m p r , where 
r is a length scale implicitly defined through the static potential ||164|| , V(r), 



dV{r) 



dr 



1.65. 



(3.25) 



r=r 



From bottomonium phenomenology [|164J , |30|| , we can assign the experimental value 



r l = (394 ±20) MeV, while m p 



770 MeV. The data on m p has been obtained in the 



quenched approximation to QCD, by use of the Wilson fermionic and gluonic action 



by the GF11 and CP-PACS collaborations [158, 97]. The corresponding r values have 



been obtained from the interpolating formula of the ALPHA collaboration ||165|| for 
5.7 < (3 < 6.57, 



a/r = exp {- [d + d x {(5 - 6) + d 2 ((3 - 6) 2 + d 3 (/3 - 6) 3 ] } 



(3.26) 



with d = 1.6805, di = 1.7139,4 = -0.8155, d 3 = 0.6667. 

The leading order scaling violations of m p r are expected to be proportional to the 
lattice spacing, a. The data points cover the range, 5.7 < (3 < 6.47, or, 0.17 fm > 
a > 0.047 fm. Only the CP-PACS results have been used in the linear fit. In the 
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continuum limit the ratio m p ro deviates from the phenomenological estimate by about 
15 %, indicating the limitations of the quenched approximation. In Ref. |97| deviations 
of some quenched ratios between masses of light hadrons from experiment of up to 10 % 
have been observed. 

Due to the substantial slope of the extrapolation, the result obtained on the finest 
lattice with a resolution of about 4 GeV still deviates by almost 10 % from the continuum 
limit extrapolated value. This is different from the situation regarding the glueball 
spectrum where leading order lattice artefacts are proportional to a 2 . In Figure |3. 4 
we display the continuum limit extrapolation for the lightest quenched glueball mass 
that has scalar quantum numbers, J PC = ++ . The (3 range covered in the Figure, 
5.7 < (3 < 6.4, is about the same as that of Figure |3.3| . However, within statistical 
errors, the f3 = 6.4 results are compatible with the continuum limit and this despite the 
fact that the scalar glueball behaves rather pathologically ||116|| in the sense that the 
slope of this extrapolation is much larger than in any other of the glueball channels. 
The continuum limit extrapolated mass comes out to be m(0 ++ ) = 1.485(35) GeV or 
m(0 ++ ) = 1.720(50) GeV, depending on whether the scale is set from the p-mass or r , 
respectively; clearly, the dominant source of uncertainty is quenching. 
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Figure 3.5: The scale, r , in lattice units against the coupling (3. 
Refs. pi (ALPHA), |J|, |167| (BS) and [168[ (EHK). 



The data are from 



In Figure |3T5| , we plot r 1 a obtained from quenched Wilson action simulations [ |165 . 



166| , [167], |168|| versus the bare coupling, (3. The results are also displayed in Table |37T . 



Within the range, 5.5 < (3 < 6.8, the lattice spacing varies by a factor of about 7. The 
interpolating curve for 5.7 < (3 < 6.57, Eq. (|3.26| ), is included into the plot as well as an 
estimate obtained by converting the result [|103|| , A^r = 0.294(24), into the bare lattice 
scheme [p.69|l at high energy (1000 r^ 1 ) and running the coupling down to lower scales 
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Table 3.1: The scale tq in lattice units, obtained from 577(3) gauge theory simulations 
with Wilson action. 



p 






a 






Bra (bs) 




168H ( E BK) 




165| (ALPHA) 


Eq. (|3726|) 


5.5 


2.005(29) 








5.54 




2.054(13) 






5.6 


2.439(62) 


2.344(8) 






5.7 


2.863(47) 


2.990(24) 


2.922(9) 


2.930 


5.8 


3.636(46) 




3.673(5) 


3.668 


5.85 




4.103(12) 




4.067 


5.9 


4.601(97) 






4.483 


5.95 






4.808(12) 


4.917 


6.0 


5.328(31) 


5.369(9) 




5.368 


6.07 






6.033(17) 


6.030 


6.2 


7.290(34) 




7.380(26) 


7.360 


6.3 


8.391(72) 






8.493 


6.4 


9.89(16) 




9.74(5) 


9.760 


6.57 






12.38(7) 


12.38 


6.6 


12.73(14) 






12.93 


6.8 


14.36(8) 









via Eq. (|3.23|) , using the three loop approximation of the /3-function, Eqs. (|3.22f) , (|C.2|) , 
( p.3| ) and ( |C.5| ). Taking into account the logarithmic scale, deviations from asymptotic 
scaling are quite substantial, at least for (3 < 6.4. One of the reasons for this failure of 
perturbation theory at energy scales of several GeV are large renormalisations of the 
lattice action [ |119|| , due to contributions from tadpole diagrams [|120|| . One might hope to 



partially cancel such contributions by defining an effective coupling | |119| , |170| , |171| , |172j , 
|173| , |120|| from the average plaquette value, measured on the lattice and, indeed, such a 
procedure somewhat reduces the amount of violations of asymptotic scaling ||172j, |173| . 



4 The static QCD potential 

We shall introduce the Wegner- Wilson loop and derive its relation to the static potential. 
Subsequently, expectations on this potential from exact considerations, strong coupling 
and string arguments as well as perturbation theory and quarkonia phenomenology are 
presented. Lattice results are then reviewed. Finally, the behaviour of the potential 
at short distances, the breaking of the hadronic string and aspects of the confinement 
mechanism are discussed. 
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4.1 Wilson loops 



The Wegner- Wilson loop has originally been introduced by Wegner ||174|| as an order 
parameter in Z 2 gauge theory. It is defined as the trace of the product of gauge variables 
along a closed oriented contour, 6C, enclosing an area, C, 



W{C) = Tr \V 



cxp 



8C 



Tr 



n u x> , 



(4.1) 



While the loop, determined on a gauge configuration, {U X) ^}, is in general complex, 
its expectation value is real, due to charge invariance: in Euclidean space we have, 
(W(C)) = (W*(C)) = (W(C))* = 0. It is straight forward to generalise the above 
Wilson loop to any non-fundamental representation, D, of the gauge field, just by re- 
placing the variables, U X)fl , with the corresponding links, . The arguments below, 
relating the Wilson loop to the potential energy of static sources go through, indepen- 
dent of the representation according to which the sources transform under local gauge 
transformations. In what follows, we will denote a Wilson loop, enclosing a rectangular 
contour with one purely spatial distance, r, and one temporal separation, t, by W(r,t). 
Examples of Wilson loops on a lattice for two different choices of contours, 5C, are 
displayed in Figure [lT| 




A 




Figure 4.1: Examples of rectangular on- and off-axis Wilson loops with temporal extent, 
t = 5a, and spatial extents, ri = 3a, and, r 2 = 2^/2 a, respectively. 

In Wilson's original work the Wilson loop has been related to the potential 
energy of a pair of static colour sources, by use of transfer matrix arguments. However, 
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it took a few years until Brown and Weisberger attempted to derive the connection 
between the Wilson loop and the effective potential between heavy, not necessarily 
static, quarks in a mesonic bound state |]175||. Later on mass dependent corrections to 



the static potential have been derived along similar lines p2| , In Section |6.3| , we will 
discuss these developments in detail. Here, we derive the connection between a Wilson 
loop and the static potential between colour sources which highlights similarities with 
the situation in classical electrodynamics and which is close to Wilson's spirit. 
For this purpose we start from the Euclidean Yang-Mills action, Eq. ( |B.14| ), 

S = i/<ftrPX- (4-2) 



Ag 2 J w w 

The canonically conjugated momentum to the field, A%, is given by the functional 
derivative, 

<=mkr 7- F '- = >■ (4 - 3) 

The anti-symmetry of the field strength tensor implies, 7rf = 0. In order to obtain a 
Hamiltonian formulation of the gauge theory, we fix the temporal gauge, A\ = 0. In 
infinite volume, such gauges can always be found. On a toroidal lattice this is possible 
up to one time slice t', which we demand to be outside of the Wilson loop contour, 

a > t. 

The canonically conjugated momentum, 

vr" = -i * (4.4) 

now fulfils the usual commutation relations, 

[A a j y il ]=t5 m 5 a \ (4.5) 
and we can construct the Hamiltonian, 

H = Jd 3 x ^d 4 A; - ^F%F*\ =\jd 3 x {El El - B?B?) , (4.6) 

that acts onto states, ^[A^]. In Euclidean metric, the magnetic contribution to the 
total energy is negative. 

A gauge transformation, Q, can for instance be represented as a bundle of SU(N) 
matrices in some representation D, flo(x) = e >ua ^ T o. We wish to derive the operator 
representation of the group generators, TJj, that acts on the Hilbert space of wave 
functionals. For this purpose we start from, 



tfW(x)T£(x) + 



y = y + 5V. (4.7) 



From Eq. ( |3.4j ) one easily sees that, 5Ai = A^ — Ai = —(diU + i[Ai,uj]). We obtain, 
6 ^= f d3x TTT^ 5 Mx)= tew(x)A^T = -- tea; a (x)(A^) a (x)^, (4.8) 
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where we have performed a partial integration and have made use of the equivalence, 

" -~E U (4.9) 



SAi g 

of Eqs. ( fO|) and ( [4.4|) . Hence we obtain the representation 



T£ = -^-(D l E l ) a : (4.10) 

the covariant divergence of the electric field operator is the generator of gauge transfor- 
mations! 

Let us assume that the wave functional is a singlet under gauge transformations, 
i?(ft)$L4 M ] = #L4 M ]. This implies, 

(A30°* = o, (4.ii) 

which is Gaufi' law in the absence of sources: ^ lies in the eigenspace of D^Ei that 
corresponds to the eigenvalue zero. Let us next place an external source in fundamental 
representation of the colour group at position r. In this case, the associated wave 
functional, ty a , a = 1, . . . , N, transforms in a non-trivial way, 

[R(n)*]a = 0^/3, (4.12) 

This implies, 

(A#i)°* = -g5 3 (r)T a V, (4.13) 

which again resembles Gaufi' law, this time for a point-like colour charge at position^ 
r. For non-fundamental representations, D, Eq. fl4.13| ) remains valid under the replace- 



ment, T a -^T%. 

Let us now place a fundamental source at position and an anti-source at position 
r. The wave functional, \l/ r , which is an N x N matrix in colour space will transform 
according to, 

K a p = ^ 7 (0)^(r)* r , 75 . (4.14) 

One object with the correct transformation property is a gauge transporter (Schwinger 
line) from to r, 

1 „ + , . 1 



= -j=U*(r,t) = -j=V exp (i J' dxA(x,t) 



(4.15) 

which on the lattice corresponds to the ordered product of link variables along a connec- 
tion between the two points. Since we are in temporal gauge, A±(x) = 0, the correlation 
function between two such lines at time-like separation, t, is the Wilson loop, 

(W(r,t)) = l^M^M)), (4.16) 



9 Of course, on a torus, such a state cannot be constructed. Note also that in our Euclidean space- 
time conventions Gaufi' law reads, D,Ej(x) = — p(x), where rho denotes the charge density. 
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which, being a gauge invariant object, will give the same result in any gauge. Other 
choices of \l/ r , e.g. linear combinations of spatial gauge transporters, connecting with 
r, define generalised (or smeared) Wilson loops, W^(r,t). 

Following the discussion of Section |3.4| , we insert a complete set of transfer matrix 
eigenstates, \$ r ,n), within the sector of the Hilbert space that corresponds to a charge 
and anti-charge in fundamental representation at distance r, and expect the Wilson 
loop in the limit, L T a ^> t, to behave like, 

(W*(r,t)) =O $ r,n|* r |0)|V^, (4.17) 

n 

where the normalisation convention is such that, ($ n |$„) = (\Edv|/) = 1 ; anc l the com- 
pleteness of eigenstates implies, J2 n l(^n|^|0)| 2 = 1. Note that no disconnected part 
has to be subtracted from the correlation function since \l/ r is distinguished from the 
vacuum state by its colour indices. E n (r) denote the energy levels. The ground state 
contribution, Eq(t), that will dominate in the limit of large t can be identified as the 
static potential. 

The gauge transformation properties of the colour state discussed above, which 
determine the colour group representation of the static sources and their separation, r, 
do not yet completely determine the state in question: the sources will be connected by 
an elongated chromo-electric flux tube. This vortex can for instance be in a rotational 
state with spin A ^ about the inter-source axis. Moreover, under interchange of the 
ends the state can transform evenly (g) or oddly (u). Finally, in the case of A = 0, it can 
transform symmetrically or anti-symmetrically under reflections with respect to a plane 
containing the sources. It is possible to single out sectors within a given irreducible 
representation of the relevant cylindrical symmetry group ||176|| , D^, with an adequate 



choice of A straight line connection between the sources corresponds to the -Dooh 
quantum numbers, £+. Any static potential that is different from the £+ ground state 
will be referred to as a "hybrid" potential. Since these potentials are gluonic excitations 
they can be thought of as being hybrids between pure "glueballs" and a pure static- 
static state; indeed, high hybrid excitations are unstable and will decay into lower lying 
potentials via radiation of glueballs. We will address the question of hybrid potentials 



in detail in Sections 5.2 and 5.3 



4.2 Exact results 



We identify the static potential, V(r), with the ground state energy, E (r), of Eq. ( |4.17| ) 
that can be extracted from the Wilson loop of Eq. (|4.1|) via Eq. fl3.20| ). By exploiting 
the symmetry of a Wilson loop under an interchange of space and time directions, it 
can be proven that the static potential cannot rise faster than linearly as a function of 
the distance r in the limit, r — > oo [ 177| . Moreover, reflection positivity of Euclidean 
n-point functions ||155| , p7| implies convexity of the static potential [ |178|| , 



V"(r) < 0. 



(4.18) 
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The proof also applies to ground state potentials between sources in non-fundamental 
representations. However, it does not apply to hybrid potentials since in this case the 
required creation operator extends into spatial directions orthogonal to the direction of 
r. Due to positivity, the potential is bound from below[^]. Therefore, convexity implies 
that V(r) is a monotonically rising function of r, 



V'(r) > 0. 



(4.19) 



In Ref. [ 179 |, which in fact preceded Ref. | 178 |, somewhat more strict upper and 
lower limits on Wilson loops, calculated on a lattice, have been derived: let a a and 
a T be temporal and spatial lattice resolutions. The main result for rectangular Wilson 
loops in representation D and d space-time dimensions then is, 

(W{a a , a T )) rt/(a ^ ) < (W{r, t)) < (1 - c y 1^1^-2^ ( 420 ) 
with c = exp[— 4(d — 1)D(3\. The resulting bounds on V(r) for r > a a read, 

- ln(l - c) < a T V(r) < -— \a(W(a a , a T )); (4.21) 

in consistency with Ref. yi77[| , the potential (measured in lattice units, a T ) is bound 
from above by a linear function of r and it takes positive values everywhere. 



4.3 Strong coupling expansions 

Expectation values, Eq. (|3.15|) , can be approximated by expanding the exponential of 
the action, Eq. ( p,10|) , in terms of /3, exp(— f3S) = 1 — (3S + • • ■. This strong coupling 
expansion is similar to a high temperature expansion in statistical mechanics. When 
the Wilson action is used each factor, /3, is accompanied by a plaquette and certain 
diagrammatic rules can be derived |T], |180| , |181| , |182| , |183|| . Let us consider a strong 



coupling expansion of the Wilson loop, Eq. (|4.1| ). Since the integral over a single group 
element vanishes, 

"dUU = 0, (4.22) 

to zeroth order, we have, (W) = 0. To the next order in (3, it becomes possible to 
cancel the link variables on the contour, 5C, of the Wilson loop by tiling the whole min- 
imal enclosed (lattice) surface, C, with plaquettes. Hence, one obtains the expectation 



value 182, Ba M 



[(3/4}- aTca{sc) + ---,N = 2 
[/3/2N 2 ]- ama(SC) + ---,N >2 



(w(c)) = { rJ:^ m -™w)" „ > ( 4 - 23 ) 



10 The potential that is determined from Wilson loops depends on the lattice cut-off, a, and can 
be factorised into a "physical" potential V(r) and a (positive) self energy contribution: V(r,ci) = 



V(r) + V sc \f(a). The latter diverges in the continuum limit (see Section 4.5). While the "physical" 



potential, V{r), will become negative at small distance, V(r,a) is indeed non-negative. 
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for SU (N) gauge theory. area(<5C) denotes the area of the minimal lattice world sheet 
that is enclosed by the contour SC. 

If we now consider the case of a rectangular Wilson loop that extends r/a lattice 
points into a spatial and t/a points into the temporal direction, we find the area law, 



(W(r,t)) =exp[-a d rt] + 



with a string tension, 



a d a 



-d In — . 

18 



(4.24) 



(4.25) 



The numerical value of the denominator applies to SU(3) gauge theory; the poten- 
tial is linear with slope, ad, and colour sources are confined at strong coupling, d = 
(|ri| + | 1 + l r 3|)/ r > 1 denotes the ratio between lattice and continuum norms and 
deviates from d — 1 for source separations, r, that are not parallel to a lattice axis. 
The string tension of Eq. ( 4.25|) depends on d and, therefore, on the lattice direction; 
0(3) rotational symmetry is broken down to the cubic subgroup Oh- The extent of 
violation will eventually be reduced as one increases (3 and considers higher orders of 
the expansion. Such high order strong coupling expansions have indeed been performed 
for Wilson loops ||184|| and glueball masses 



185] 



Unlike standard perturbation theory, 
186| , |187|| , the strong coupling 



[156] and, therefore, has a finite radius of conver- 



whose convergence is known to be at best asymptotic 
expansion is analytic around j3 
gence. 

Strong coupling SU(3) gauge theory results seem to converge for |88[] (3 < 5. One 
would have hoped to eventually identify a crossover region of finite extent between the 
validity regions of the strong and weak coupling expansions [|188|| , or at least a tran- 
sition point between the leading order strong coupling behaviour, a 2 oc — ln(/3/18), of 
Eq. ( |4.25[ ) and the weak coupling limit, a 2 oc exp[— 27r/3/(3/3 )], of Eq. ( |3.24j ). However, 
even after re-summing the strong coupling series in terms of improved expansion param- 
eters and applying sophisticated Pade approximation techniques |189| ], nowadays such 
a direct crossover region does not appear to exist, necessitating one to employ Monte 
Carlo simulation techniques. One reason for the break down of the strong coupling 
expansion around j3 ~ 5 seems to be the roughening transition that is e.g. discussed in 
Refs. [19U, 191 1 ; while at strong coupling the dynamics is confined to the minimal area 
spanned by a Wilson loop (plus small "bumps" on top of this surface), as the coupling 
decreases, the colour fields between the sources can penetrate over several lattice sites 
into the vacuum. 



We would like to remark that the area law of Eq. ( |4.24j ) is a rather general result for 
strong coupling expansions in the fundamental representation of compact gauge groups. 
In particular, it also applies to U(l) gauge theory which we do not expect to confine in 
the continuum. In fact, based on duality arguments, Banks, Myerson and Kogut ] |192j 
have succeeded in proving the existence of a confining phase in the four- dimensional 
theory and suggested the existence of a phase transition while Guth ] |193|] has proven 
that, at least in the non-compact formulation of U(l), a Coulomb phase exists. Indeed, 
in numerical simulations of (compact) U (1) lattice gauge theory two such distinct phases 
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were found [ 194 , |195|| , a Coulomb phase at weak coupling and a confining phase at strong 
coupling. The question whether the confinement one finds in SU(N) gauge theories in 
the strong coupling limit survives the continuum limit, (3 —>■ oo, can at present only be 
answered by means of numerical simulation. 



4.4 String picture 



The infra-red properties of QCD might be reproduced by effective theories of interacting 
strings. String models share many aspects with the strong coupling expansion. Origi- 
nally, the string picture of confinement has been discussed by Kogut and Susskind [34 



as the strong coupling limit of the Hamiltonian formulation of lattice QCD. The strong 
coupling expansion of a Wilson loop can be cast into a sum of weighted random deforma- 
tions of the minimal area world sheet. This sum can then be interpreted to represent a 
vibrating string. The physical picture behind such an effective string description is that 
of the electric flux between two colour sources being squeezed into a thin, effectively one- 
dimensional, flux tube or Abrikosov-Nielsen-Olesen (ANO) vortex | |196| , |197| , |198] , |199| . 
As a consequence, this yields a constant energy density per unit length and a static 
potential that is linearly rising as a function of the distance. 

One can study the spectrum of such a vibrating string in simple models ||200| , |19C| , |47| . 
Of course, the string action is not a priori known. The simplest possible assumption, 
employed in the above references, is that the string is described by the Nambu-Goto 
action |)8[ |69| in terms of (d — 2) free bosonic fields associated to the transverse degrees 
of freedom of the string. In this picture, the static potential is | |20(J| , |201| (up to a 
constant term) given by, 



V{r) 



ar\ 1 



(d- 2)vr 



(d — 2)tt (d-2) 2 vr 2 



12 a r 2 



ar 



24 r 



1152a r 2 



(4.26) 



while for a fermionic string ||202|| one would expect the coefficient of the correction term 
to the linear behaviour to be only one quarter as big as the Nambu-Goto one above. In 
the bosonic string picture, excited levels are separated from the ground state by, 



y n 2 ( r ) = V 2 {r) + {d-2)ima 



V(r) + 



(d - 2)wn 



n 2 



2r 



(4.27) 



with n assuming integer values. It is clear from Eq. Q4.26j ) that the string picture at 
best applies to distances, 



r r c 



\d-2)n 
12 a 



(4.28) 



In four dimensions one obtains, r c « 0.33 fm, from the value, \fo ~ 430 MeV, from the 
p, a 2 , . . . Regge trajectory. 

The expectation of Eq. ( 4.26Q has been very accurately reproduced in numerical 
simulations of Zi gauge theory in d — 3 space-time dimensions ||203|| . In a recent 
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study of d — 4 SU(3) gauge theory ||114j| the hybrid potentials have been found to 
group themselves into various bands that are separated by approximately equi-distant 
gaps at large r. However, up to distances as large as 3 fm these gaps seem to be 
inconsistent with 7r/r, the expectation of Eq. ( |4.27|) . These newer data contradict 



earlier findings in SU (2) gauge theory ||146|| where good agreement with the Nambu- 
Goto string picture has been reported, such that we do not regard this issue as finally 
settled. The consistency of lattice data with Eq. ([4.271) at large separations would 
support the existence of a bosonic string description of confining gauge theories in the 
very low energy regime ||204j , [205| , [206| , |207| , |208|| . Of course, in d < 26, the string 
Lagrangian is not renormalisable but only effective and higher order correction terms 
like torsion and rigidity will in general have to be added ||209 . 

It is hard to disentangle in d = 4 the (large distance) 1/r term, expected from string 
vibrations, from the perturbative Coulomb term at short distances. Therefore, three- 
dimensional investigations (where perturbation theory yields a logarithmic contribution) 
have been suggested ||210|| . Another way out is to determine the mass of a closed string, 
encircling a boundary of the lattice [p 1 l|j with a spatial extent, / = L a a (a torelon Ell 



for details see Appendix |B|), which is not polluted by a perturbative tail. The bosonic 
string expectation in this case would be |[2i0|| , 



E n (l) = al 



{d-2)n 
61 



(4.29) 



0.66 fm. The numerical 



The naive range of validity of the picture is / ^> l c 
value applies to d — 4. An investigation of the finite size dependence of the torelon 
mass in d = 4 SU(2) gauge theory has been done by Michael and Stephenson [213 



who found excellent agreement with the bosonic string picture already for distances, 
lfm < I < 2.4 fm, quite close to l c , on the 3 % level. Qualitative agreement has 
also been reported by Teper |214|1 from simulations of SU(2), SU(3), S , f/(4) and SU(5) 
gauge theories in three dimensions. 

The bosonic string picture for r 3> /3 = aL T predicts a behaviour similar to Eq. (|4.29|) 



for the finite temperature potential, calculated from Polyakov line correlators ||215 



1 



-ln(P*(r)P(0)> 



a((3)r + ■■■, 
The Polyakov line is defined as [Eq. ( p.3| )1, 



a((3) 



a 



(d-2] 



7T 



6/3 2 



+ 



P(x) = Tr ^ T 



/ raL T 

exp I i J dx^A^x) 




(4.30) 



(4.31) 



where T denotes time ordering of the argument. The dependence of the effective string 
tension on the temperature has recently been checked for rather low T _1 = j3 < 
1.13 P c ~ 0.85 fm in a study of SU(3) gauge theory ||216|| . Although the sign of the 
leading correction term to the zero temperature limit is correct, the difference comes 
out to be bigger than predicted. It would be interesting to check whether the result will 
converge towards the string expectation at lower temperatures. 
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4.5 The potential in perturbation theory 

Besides the strong coupling expansion, which is specific to the lattice regularisation, the 
expectation value of a Wilson loop can be approximated using standard perturbative 
techniques. 

We will discuss the leading order weak coupling result that corresponds to sin- 
gle gluon exchange between the static colour sources which, although we neglect the 
spin structure, we will call "quarks" for convenience. From the Lagrangian, Cym = 
2^tr FpvF^v, one can easily derive the propagator of a gluon with four-momentum, q, 

GtM)=/-^F, (4-32) 

where /i, v are Lorentz indices and a, b = 1, . . . Na = N 2 — 1 label the colour generators. 
The same calculation can be done, starting from a lattice discretised action. The Wilson 
action, Eq. ( |3.10|) , yields the result of Eq. (|4.32| ), up to the replacement, 

^-^sinpfi). (4.33) 

Other lattice actions yield slightly different results but they all approach Eq. ( |4.32| ) 
in the continuum limit, a — > 0. Up to order a 2 , the momentum space potential can 
be obtained from the on-shell static quark anti-quark scattering amplitude: the gluon 
interacts with two static external currents pointing into the positive and negative time 
directions, A a ^ af3 = d^T®^ and A'^ a/3 = —S U ^T^ S . Hence, we obtain the tree level 
interaction kernel, 

K aP7 s(q) = ~^W 5 . (4.34) 

For sources in the fundamental representation, the Greek indices run from 1 to N 
and the quark anti-quark state can be decomposed into two irreducible representations 
of SU(N), 

N <g> N* = 1 © N A . (4.35) 
We can now either start from a singlet or an octetQ initial $/3 7 = QpQ* state, 

$J 7 = Sfr, (4.36) 

< A = ^7-^<W (4-37) 

where the normalisation is such that ^^^^ = 5^. A contraction with the group 
generators of Eq. ( |4.34|) yields, 



$J 7 T^T 7 a , = (4.38) 
1 

07 - t a/3- t 7<5 — ~ 2iV" 



K a TZbT? s = --rr^/, (4-39) 



1 We call the state Na an "octet" state, having the group SU (3) in mind. 
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where Cp = Na/(2N) is the quadratic Casimir charge of the fundamental representa- 
tion. 

We end up with the potentials in momentum space, 

V s (q) = -C F g 2 \, V (q) = = -±-V.{q), ( 4 - 4 0) 



g 2 ' 2N q 2 N 



A 



governing interactions between fundamental charges coupled to a singlet and to an octet, 
respectively: the force in the singlet channel is attractive while that in the octet channel 
is repulsive and smaller in size. 

How are these potentials related to the static position space inter-quark potential, 
defined non-perturbatively through the Wilson loop, 

V(r) = - lim 4- ln(W(r, ( 4 - 41 ) 

The quark anti-quark state creation operator, \l/ r , within the Wilson loop contains a 
gauge transporter and couples to the gluonic degrees of freedom. Thus, in general, it 
will have overlap with both, QQ* singlet and octet channels^. Since the singlet channel 
is energetically preferred, V s < V , we might expect the static potential to correspond 
to the singlet potential. 

To lowest order in perturbation theory, the Wilson loop is given by the Gaussian 
integral, 

(W(r,t)) = exp |-± / d*xd*yj;(x)G a Jl(x - y)J b u (y)] , (4.42) 
where = ±T a if (x, y) G 5C and = 0, elsewhere^- Eq. ( |4.42| ) implies for t > r, 



(W(v,t)) 



exp (c F gh j^dt' [G(r, t') - G(0, t')]j . (4.43) 



We have omitted gluon exchanges between the spatial closures of the Wilson loop from 
the above formula. Up to order a 3 s (two loops), such contributions result in terms whose 
exponents are proportional to r and r/t and, therefore, do not affect the potential of 
Eq. ( 4.41 ). Gz^^x), the Fourier transform of G^ u (q), contains the function, 

r d 4 a e iqx r°° 11 
After performing the t-integration, we obtain, 

QJ 

V(r, l 2) = -C F ^ + V SC M, (4.45) 
r 



12 Of course, for quark and anti-quark being at different spatial positions, the singlet-octet classifica- 
tion should be consumed with caution in a non-perturbative context. 

13 Note that this formula that automatically accounts for multi-photon exchanges is exact in non- 
compact QED to any order of perturbation theory. However, in theories containing more complicated 
vertices, like non-Abelian gauge theories or compact lattice U(l) gauge theory, correction terms have 
to be added at higher orders in g. 
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where a s = g 2 /(An). The piece, 



V*a(n) = C F g 2 \ 

Jq 



d 3 q 1 
9 <m (2^F? 



C F a s -fi, 



(4.46) 



that linearly diverges with the ultra-violet cut-off, /x, results from self-interactions of 
the static (infinitely heavy) sources. Beyond tree level, g 2 will depend on q, such that 
a s in momentum space has to be replaced by a s (q*) with some effective ?*(//). We find, 



where 



V{q) = V s (q), 



V(q, 0) = / d 6 r e l * r V(r), V{r) = V(t, fi) - V sel{ (fM). 



(4.47) 



(4.48) 



This self-energy problem is well known on the lattice and has recently received at- 
tention in continuum QCD, in the context of renormalon ambiguities in quark mass 
definitions gT§ . 

At order a£ a class of diagrams appears in a perturbative calculation of the Wilson 
loop that results in contributions to the static potential that diverge logarithmically with 
the interaction time ||219|| . In Ref . ||220|| , within the framework of effective field theories, 
this effect has been related to ultra-soft gluons due to which an extra scale, V Q — V s , is 
generated. Moreover, a systematic procedure has been suggested to isolate and subtract 
such terms to obtain a finite interaction potential between heavy quarks. However, one 
would wish to understand and regulate such contributions not only for heavy quarks but 
also in the static case. At present it is not clear whether the interaction potential within 
a heavy quark bound state whose effective Hamiltonian contains a kinetic term will, in 
the limit of infinite quark masses, approach the static potential that is defined through 
the Wilson loop. Hence, one should carefully distinguish between the static and heavy 
quark potentials. We shall discuss a physically motivated reason for the breakdown of 
standard high order perturbative calculations of the Wilson loop in Section L8 . In our 



opinion the presence of a low energy scale, which we shall identify with the gap between 
ground state potential and hybrid excitations, results in problems within perturbation 
theory in the limit of large t. 

That something in the position space derivation of the perturbative potential might 
be problematic is reflected in Eq. ( J4.43| ) that contains an integration over the interaction 
time. We know for instance from the spectral decomposition of Section !0 that for 
any fixed distance r, Wilson loops will decay exponentially in the limit of large t. 
However, the tree level propagator in position space is proportional to, (r 2 + 1 2 ) -1 , i.e. 
asymptotically decays with t~ 2 only. We notice that the integral receives significant 
contributions from the region of large t as demonstrated by the finite t ^> r tree level 
result, 



-\n(W(r, t)) 



C F a s 2 ( 

-t— < arctan 



7T [ 



r 
2f 



In 1 + 



+ (r + t)V selt . (4.49) 
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Ignoring this problem for the moment, one finds the weak coupling equality, Eq. ( f4.47| ), 
to hold up to two loops (order a 3 ) in perturbation theory. Some of the hybrid potentials 
of Section E72] that can be extracted from generalised Wilson loops, (Wy), in which the 
wave function, transforms non-trivially under the cylindrical rotation group D^, 
however, receive leading order octet contributions. This is because the creation operator, 
explicitly couples to the gluonic background. 

The tree level lattice potential can easily be obtained by replacing by g M and (in 
the case of finite lattice volumes) the integrals by discrete sums over lattice momenta, 



2tt rii L a L a 

Qi = , ni = hi,..., — • 

L a a' 2 ' ' 2 



The lattice potential reads, 



where 



V(r) = Keif (a) - C F a s 



An 



E 



Lla? ^ Ei Mi 



(4.50) 



(4.51) 



(4.52) 



and V se if(a) = Cpa s [1/0]. We have neglected the zero mode contribution that is sup- 
pressed by the inverse volume, (aL a )~ 3 . In the continuum limit, [1/r] approaches 1/r 
up to quadratic lattice artefacts whose coefficients depend on the direction of r while 
Keif (a) with nf flavours of Wilson fermions diverges like ||221| , |222|| , 



Keif(a) = Cpa^ 1 [3.1759115 . . . + (16.728 ... - 0.423 . 



O.s 



(4.53) 



The numerical values apply to the limit, L a — ► oo and, in the case of the one loop 
coefficient, iV = 3. Note that under the substitution, fi « 1.5879557 ir/a, the tree level 
term of Eq. ( |4.53| ) is identical to Eq. ( |4.46| ) . A one loop computation of on-axis lattice 
Wilson loops in pure gauge theories can be found in Ref. [221 1. The tree level form, 



Eq. ( f4.51p , is often employed to parameterise lattice artefacts. 

Besides defining the static potential from Wilson loops, on a volume with temporal 
extent, (3 = aL T , and periodic boundary conditions it can be extracted from Polyakov 
line correlators^, 



V(r) 



— lim 

/3->oo 



d 
dp 



<P*(r)P(0)> 



(4.54) 



at any given time the pair of Polyakov lines has the gauge transformation properties of a 
static quark anti-quark pair and, thus, the ground state is the same as that of a Wilson 
loopf5 In the Polyakov line correlator, no projection is made onto the ground state 
of the flux tube. Therefore, one might expect 



(P*(r)P(0)) 



1 

iV2 



+ N A e 



-/3V (r) 



(4.55) 



The Polyakov line is denned in Eq. ( [4.31 ) 



5 This statement is not entirely correct on a finite spatial volume as we shall see in Section 4.7.3 
However, for distances, r, with n < aL a /2, the ground state is indeed the same. 
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where the "octet" potential, V , can be thought to be related to hybrid excitations of the 
inter-quark string. At small j3 (high temperature) the exponentials can be expanded and 
the term proportional to g 2 vanishes due toV s = —NaV : the leading order r dependent 
contribution to the correlation function requires two gluons to be exchanged, 

(p*(r)p(o)> = (1 + §^ a2 $) e ~ PVsclc - (456) 

The above result can also be produced by a direct perturbative evaluation of the 
Polyakov line correlator in position space: the correlation function contains two disjoint 
colour traces, therefore, single gluon exchanges only contribute to the self-energy. The 
colour factor that accompanies two gluon exchanges is, -^tr (T a T b )S ac 5 bd j^ti (T c T d ) = 



jj^. Hence, we indeed reproduce Eq. (|4.56Q . By assuming the singlet channel (V s < V a ) 



to dominate Eq. (|4.55|) in the asymptotic limit of large (3 one obtains the result of 
Eq. ( |4.45| ), i.e. the same potential as from Wilson loops. However, if we insist on 
perturbation theory to hold for the correlation function itself at large (3, i.e. at low 
temperature, a misleading (and divergent) result is obtained. We have demonstrated 
that extra information how to treat the limit j3 — > oo has to be provided to obtain the 
correct zero temperature tree level potential from Polyakov line correlation functions. 
We take this as an indication that in three loop calculations of the Wilson loop the 
t —>■ oo limit should be performed with caution too. 

4.6 Potential models 

Several parametrisations of the QCD potential have been suggested in the past, either 
QCD inspired or purely phenomenological. One should keep in mind that one would not 
necessarily expect a potential that reproduces the observed quarkonia levels to coincide 
with the static potential calculated from QCD, due to the approximations involved, 
namely the adiabatic and non-relativistic approximations. 

A purely phenomenological logarithmic potential, V(r) = Cln(r/ro), has been sug- 
gested as an easy way to produce identical spin-averaged charmonia and bottomonia 
level splittings ||224[| . This idea has been incorporated into the Martin potential |18|. |I~9|j . 



V(r) = C + (r/ro) Q , with a ~ 0.1. Potentials that have QCD-like behaviour built in at 
small distances have been suggested for instance in Refs. JT2], [T3[ . We have already 
discussed the prototype Cornell potential |12| , V(r) = — e/r + ar, that interpolates 
between perturbative one gluon exchange for small distances and a linear confining be- 
haviour for large distances. Another elegant interpolation between the two domains, 
containing the one loop running of the QCD coupling, 



av(q)= Wv' tv=h \k)' (4 - 57) 

has been suggested by Richardson |T5[] : in momentum space, ty is substituted by, 
t' v = ln(l + q 2 /A v ), which does not affect the perturbative ultra-violet domain since 
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t'y — > ty as q 2 — ► oo. However, the Landau pole at q 2 
energy behaviour of the resulting potential, 



V{r) 



d 3 q 



3 «qr 



Pq J (27r)3g2 ln (i + g 2/ A 



Ay is regulated and the low 



(4.58) 



V) 



is given by, 



V(r) 



or 
2Po 



oo 



A 2 



(4.59) 
(4.60) 



i.e. the ansatz connects the QCD scale parameter, Ay to the string tension, o. We 
have neglected an infinite additional constant from Eq. ( |4.59| ) that can be eliminated 
by adding an appropriate counter term to the integrand of Eq. ( f4.58|) . 
From Eq. ( 4.60|) and the relation, 
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ax/(2j9o) 



(4.61) 
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■/VF « 0.70253(0.70048) for n ) 



we find Ajjg/^/a m 0.71639 for rif = or A^j^ 
3(4), respectively. This has to be compared to the value, A-jj^/^fo = 0.52 ± 0.05, 
determined by lattice simulations [|103j , |166|| for rif = 0. Experimental results from 
e + e~ scattering experiments at LEP and SLAC indicate somewhat bigger ratios [82], 
A^/v^ 7 = 0.88(12) MeV, for n f = 4, where we have assumed, Jo = (430 ± 20) MeV: 
while the Richardson potential overestimates the A-parameter in the quenched case 
it might approximate the experimental nj = 3 situation quite well. However, this 
coincidence is rather accidental. 

Many so-called QCD potentials have been suggested that incorporate two loop per- 
turbation theory at short distances, with varying interpolation prescriptions to different 
assumptions on the large distance behaviour. The most popular potential within this 
class is probably the Buchmuller-Tye parametrisation |20| that, like the Richardson po- 
tential, is formulated in momentum space. For collections of various parametrisations, 
we refer to Refs. [p26| , |227| |. While phenomenological potentials like a logarithmic as 
well as the Martin potential are ruled out at large and intermediate distances by lattice 
data and at short distances by pQCD, such parametrisations may still serve to explore 
the sensitivity of the heavy quark spectrum towards QCD. 

Basically, all potentials that yield a correct description for the spin-averaged quarko- 
nia spectra are only slight variations around the Cornell potential in the relevant region, 
0.2 fm < r < 1 fm. Unfortunately, the top quark is too heavy to form stable hadronic 
states and basically only the production rate of ti in e + e _ or collisions as a func- 

tion of the energy will directly depend on the potential at very short distances. Decay 
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rates and fine structure splittings of quarkonia in principle can probe the potential at 
short distances too and predictions of these quantities indeed depend very sensitively on 
the underlying ansatz [ |228| . As we will see in Sections |^ and |7^, the predictive power of 
quarkonium physics on the short range potential is reduced by theoretical uncertainties 
in the matching of an effective field theory to QCD. A big part of the (multiplicative) 
uncertainty in the fine structure, however, cancels from ratios of such splittings. 



4.7 Lattice results 



The static QCD potential has been determined to high accuracy in quenched lattice 
studies with Wilson |229|, |23(| £33], [232], |23j, £34], [23j, |2TJ, [236], £37|, [23H], [146], [23|, [M], 
147j m, [Tn| |16|, HI HH as wel1 as various improved lattice actions |43J [m], 

351 |1 14| | in SU(2) and SU(3) gauge theories. Results for QCD with sea quarks have 
been obtained in Refs. |44], |t|, |46|, |167|, g47|, |48|, |4|. After discussing the methods 



most commonly used we will present results on the potential in QCD, without and with 
sea quarks. 



4.7.1 Evaluation method 



The relative statistical errors of Wilson loop expectation values turn out to increase 
exponentially fast with the Euclidean time extent, t, of the loop. Therefore, after 

23%, [233|, |234||, replacement of the straight spatial 
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some pioneering studies ]2: 
connection within the Wilson loop by operators with improved overlap to the physical 
ground state turned out to be essential for a reliable determination of the potential at 
large distances from data at moderate t separations. For this purpose, in Refs. H231 



235| , |236|| , linear combinations of certain spatial paths connecting quark and anti-quark 
were employed. Subsequently, iterative smearing techniques [|143j, |144|, |146j, |147|| turned 



out to be extremely successful in optimising the ground state overlap. 

Among all algorithmic and technical tricks employed in lattice simulations smearing 
is certainly the most important one. The underlying concept somewhat resembles that 
of cooling techniques that are applied to extract classical properties of quantum field 
configurations | |250|| with the difference that, since smearing is a purely spatial procedure, 
the spectrum of the theory remains unaffected: fat links are iteratively constructed by 
replacing a given link by the sum of itself and the neighbouring six (in d — 3 + 1 
dimensions) spatial staples with some weight parameter, a > 0, 



U, 



SU{N) 



U, 



- -TP 



• -J 



(4.63) 



Psu(N) denotes a projection operator, back onto the SU(N) manifold. One possible def- 
inition is P7| , U = P SU (N)(A) E SU(N), RetrfMt = max. The procedure, Eq. ( ^63l) 



can be iterated several times over the whole lattice. The number of iterations and a 
represent free parameters which can be varied to optimise the overlap of an operator, 
constructed from the fat links, with the physical ground state in question. Several 
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variations of the algorithm exist. For example, all links within a given timeslice can 
be replaced at once or several subgroups can be replaced, subsequently. Blocking or 
fuzzing algorithms can be used, in which a fat link of smearing level n extends over 
more lattice sites than the previous links of level n — 1. Smearing and fuzzing can be 
combined etc.. All smearing and fuzzing methods have in common that the expectation 
value of a plaquette built from fat spatial links is increased during the iterations, similar 
to cooling, which means that the contribution to the gauge action from spatial links 
is reduced: the movement of the magnetic field through colour space under a change 
of the spatial position is minimised. Operators, built from such fat links, are likely to 
effectively decouple from excitations since the ground state wave function is always the 
smoothest wave function within any given channel. Smearing or fuzzing methods can 
be combined with variational minimisation techniques when determining a correlation 
matrix between a set of different operators [|146j| , to achieve further improvement. 

The potential is finally extracted from expectation values of smeared Wilson loops, 
W(r,t), where the spatial transporters are constructed from fat links, 

V(t) = limV(r,t), (4.64) 

t— >oo 

V(r, t ) = -lh W r,.) )s ,-h J^ )) . (4.65) 

a T denotes the temporal lattice spacing. On the lattice, the limit of large temporal 
separation is approximated by a single- or multi-exponential fit to Wilson loops for 
a range, t > t m in( r )- Positivity of the transfer matrix implies that V(r,t) converges 
towards the asymptotic value, V(r), monotonically from above, a feature that is essential 
for the reliable detection of saturation of effective masses, V(r,t), into a plateau. In 
general, within given statistics, t m m(r) will depend on the distance, r. Within the typical 
window of lattice spacings, 0.2 fm > a > 0.05 fm, in pure gauge theories and standard 
smearing and simulation techniques, this dependence happens to be weak. However, 
this does not necessarily have to be so but depends very much on the interplay between 
the dynamics of the underlying theory, smearing methods and statistical errors. 

In order to illustrate the importance of a careful analysis of the t dependence of the 
lattice data we consider the case of an unsmeared on-axis Wilson loop on an isotropic 
lattice, (W(r,t)) = (W(t,r)). For t > r, we expect (W(r,t)) oc e _y( ^ r ^. The symmetry 
under interchange of r and t implies, (W(r,t)) oc e~~ v ^ r for r t. This means, 

V(r, t) = a eS (t)r, a eS (t) = V\t) (4.66) 

Thus, approximating V(r) by an effective potential, V(r, t min ), with an r-independent 
value of t m i n automatically implies a linear rise |[242j , [251|| within the region, r>f 



mm; 



for any potential with non-vanishing derivative. This illustrates the importance of 
separately investigating the approach to the plateau for each distance. Let us examine 
closely the situation for the Cornell potential, V(r) = V se a + ar — e/r. In this case, 
taking one and the same t-value for all separations we find, 

a eS (t) = a + (4.67) 
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even a pure Coulomb potential, a = 0, implies a non-vanishing <7 e fr at finite i < r. 
Of course, the symmetry of the Wilson loop under interchange of r and t also implies 
that no plateau in V(r,t) can be found, unless i>r. For smeared Wilson loops, one 
would still expect a similar 1/t 2 approach (with a different coefficient) of cr c g towards 
the asymptotic limit, while effective masses, V(r, t), will approach V(r) exponentially 
fast at any r. 

4.7.2 The quenched potential 



o 

> 




Figure 4.2: The quenched Wilson action SU(3) potential, normalised to V(r ) = 0. 



In Figure [12[ we display the quenched potential, obtained at three different (3 values 
in units of r$ « 0.5 fm from the data of Refs. [|173| , P9f| . The lattice spacings, determined 
from ro, correspond to a ~ 0.094 fm, 0.069 fm and 0.051 fm, respectively. The curve 
represents the Cornell parametrisation with e = 0.295. At small distances the data 
points lie somewhat above the curve, indicating a weakening of the effective coupling 
and, therefore, asymptotic freedom. We will discuss this observation later. All data 
points for r > 4a collapse onto a universal curve, indicating that for (3 > 6.0 the scaling 
region is effectively reached for the static potential. Moreover, continuum rotational 
symmetry is restored: in addition to on-axis separations, many off-axis distances of the 
sources have been realised and the corresponding data points are well parameterised by 
the Cornell fit for r > 0.6 tq. Prior to comparison between the potential at various (3, 
the additive self-energy contribution, associated with the static sources, that diverges 
in the continuum limit has been removed. This is achieved by the parametrisation- 
independent normalisation of the data to V(r ) = 0. 
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Figure 4.3: Comparison between the tree level expectation of violations of rotational 
invariance and lattice data at (3 — 6.4. 



The lattice potential at (3 < 6.5 is well described by the functional form ||241]| , 



V(r) = Kont(r) - l5V(r), V cont (r) = V + ar-- + %, (4.68) 

for separations as small as r > ^a. The 1/r 2 term is not physically motivated but 
effectively parameterises the weakening of the coupling with the distance while the 
difference between tree level lattice and continuum perturbation theory results, 



SV(r) 



-, (4.69) 
r 



is used to quantify lattice artefacts. In Figure we compare the theoretical differ- 
ence 5V(r) to 5V(r) = [V con t(r) — V(r)]/l, as calculated from the lattice data after 
determination of the fit parameters, Vq, <j, e, / and I, at /3 = 6.4, in lattice units ||166| . 



The Figure demonstrates that at the level of precision achieved, deviations from the 
continuous fit curve are statistically significant for r < 4a. Moreover, deviations from 
Kont( r ) are qualitatively indeed very well parameterised by a multiple of the tree level 
difference. 



4.7.3 Finite size effects 

In lattice simulations the potential is determined on a torus with finite volume and 
the question of finite size effects (FSE) arises. Obviously, the ground state potential 
is affected by the infra-red cut-off. For instance, by exploiting the r «-> t symmetry of 
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Euclidean space-time, it is clear from Appendix [D] that for the extreme case of a spatial 
extent smaller than the critical temperature of the deconfmement phase transition, 
any asymptotic string tension will disappear. The other source of FSE on Wilson 
loops is related to un-wanted interactions of source and anti-source around the periodic 
boundaries that will become negligible as L a a — ► oo: by unwrapping the spatial torus 
onto an infinite hyper-cubic lattice of cells with spatial periods, L a a, it becomes obvious 
that each charge at position r is accompanied by an infinite set of mirror charges at 
r n = r+nLo-a, integer. For on-axis geometries, r = ri, for instance, the closest mirror 
charge will be separated by a distance L a a—r from the origin, followed by another charge 
at L a a + r. Therefore, in this case one would naively expect V(r) = V(L a a — r). The 
symmetry, V(r) = V(r n ), is indeed reflected in the tree level weak coupling expansion 
results of Eqs. ( |4.51| ) and (|4.52|) and in fact holds to any order of perturbation theory 
for the singlet and octet potentials. 

General considerations, based on the centre symmetry of the action which is dis- 
cussed in Appendix [D], however, lead us to expect the potential to be more robust 
against FSE than perturbation theory suggests. We will assume the source to be at 
position and the anti-source to be at r. The spatial connection, built into a 
(smeared) Wilson loop is a linear combination of products of link variables along paths 
that have trivial winding number around the periodic boundaries. has definite eigen- 
values, Z{ G Zjy, with respect to the centre transformations associated with the three 
spatial directions. If we place the hyperplanes at which centre transformations, z, are 
applied at position x with xi > r*, i.e. such that they do not interfere with the shortest 
connection between the two test charges, we have Zi = 0. However, a gauge transporter 
to a mirror charge at r n has the centre transformation property, = z Ui . Since the cen- 
tre symmetry is both, a symmetry of the action as well as of the path integral measure, 
Zi are conserved quantum numbers: the creation operator, \l/ r , only couples to mirror 
charges at distances in which all n^'s are multiples of N. For on-axis separations in 
SU(N) this means that the closest mirror charge contributing to the Wilson loop will 
be at a distance NL a a — r |[236 , |242|| , rather than L a a — r as the geometric argument 



alone or perturbation theory would have suggested. This suppression of FSE does obvi- 
ously not work for Polyakov loop correlators in which the state of the gluonic flux tube 
remains unspecified. While in the standard weak coupling expansion the gauge group 
only influences the group theoretical pre-factors the centre charge affects the zero mode 
sector [|252f , q = (that is suppressed by a power of the volume). 



Numerical simulations of SU(2) and SU(3) gauge theories ||242| , |147|| suggest that 



FSE on the static potential determined from Wilson loops, even at distances as big as 
a a, are numerically undetectable on the 1 - 2 % level for L a a > 3r . A 
reason besides the protection due to the centre symmetry for this finite size friendliness 
is the rather rapid on-set of the deconfmement phase transition which is first order 
in SU (3) gauge theory. Full QCD, however, is less well explored yet and one might 
expect somewhat bigger FSE, at least for light sea quarks since the fermionic part 
of the action explicitly breaks centre symmetry. In particular, it might be hard to 
discriminate between breaking of the flux tube due to screening by sea quarks and FSE 
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as reasons for an eventually flattening potential at large distances. 



4.7.4 Sea quark effects 
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Figure 4.4: Comparison between quenched (/3 = 6.2) and un-quenched n/ — 2 (/3 
5.6, n = 0.1575) ground state and U u potentials (from Ref. |[248|| ). 



When including sea quarks, one would expect two physical effects, one at large 
distances and one at small distances. While within the quenched approximation the 
number of quarks and anti-quarks are separately conserved, with sea quarks, only the 
difference (the baryon number) is a conserved quantity. Light quark anti-quark pairs 
can be created from the vacuum and in general transitions between a colour "string" 
state, spanned between two static sources, and two static-light mesons can occur. If the 
energy stored in the colour string between the sources exceeds a certain critical value 



at some distance, 



the string will "break" and decay into two static-light mesons, 



separated by a distance, r. Therefore, in the limit, r — > oo, the ground state energy will 
stop rising with the distance and saturate at a constant level: the static sources will be 
completely screened by light quarks that pop up out of the vacuum. 

The other effect will change the potential at short distances. While the vacuum po- 
larisation due to gluons has an anti-screening effect on fundamental sources, sea quarks 
result in screening. Therefore, the running of the QCD coupling with the distance is 
slowed down with respect to the quenched approximation. This is for instance reflected 
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in the factor, 11 — 2nf, within the perturbative /3-function coefficient, (3q, of Eq. ( |C.2| ). 
When running the coupling from an infra-red hadronic reference scale down to short 
distances, the effective Coulomb strength in presence of sea quarks should, therefore, 
remain at a higher value than in the quenched case. It should be possible to detect this 
effect in the coefficient, e, within the Cornell parametrisation. 
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Figure 4.5: Sea quark mass dependence of the effective Coulomb strength of the static 
inter-quark potential. 



In Figure |4.4| , a recent comparison between the quenched ((3 = 6.2) and rif = 2 
static potential ((3 = 5.6, k = 0.1575) by the T\L collaboration at a sea quark mass, 
fn u d ~ rn s /2, is displayed ||253| , [248|| . Besides the ground state potential, £+, the 
lowest lying hybrid potential, U u , is shown. Estimates of masses of pairs of static-light 
mesons (2m ps and m ps + m s ) into which the static-static systems can decay are also 
included into the Figure. The potentials have been matched to each other at a distance, 
r = ro ~ 0.5 fm. Around r m 2.3 ro ~ 1.15 fm, both un-quenched potentials, E+ and 
H u , are expected to become unstable. However, the data are not yet precise enough to 
resolve this effect. A similar comparison between the static potential and 2m ps has first 



been performed by the UKQCD collaboration [254]. 

At small r the un-quenched data points are found to be somewhat below their 
quenched counterparts: the effective Coulomb force indeed remains stronger. To quan- 
tify this effect, we fit the potentials [|46|, |55], [167|, |24§ (quenched at P = 6.0 and (3 = 6.2 
and un-quenched at (3 = 5.6 and various quark masses) for identical fit range in physical 
units, 0.4 r < r < 2r , to the parametrisation of Eqs. ( |4.68| ), ( }4.69|) , with / = 0. The 
resulting effective e values from these four-parameter fits are displayed in Figure |4.5| . 
Larger k values correspond to smaller quark masses. With two flavours of sea quarks 
of masses slightly larger than that of the strange quark, down to m U( i rs m s /3, the 
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effective Coulomb strength is increased by 17 to 22 % which is not too far off from the 
most naive expectation of 14 %, from the ratio, /^o" 7-0 ^//?^ -2 '' = 33/29. Within given 
errors the dependence on the sea quark mass cannot be resolved. However, one would 
expect this to be rather weak as the simulated quark masses are all much smaller than 
infra-red reference scales like r^ 1 . Similar results have been reported by the CP-PACS 
collaboration [247||. 



4.8 Beyond perturbation theory at short distances 

The singlet potential, 



(4.70) 



has been calculated to one loop long ago ||225|, [219|| and now the two loop result is also 
knownQ |58], |259], g6(J. It is, 



ay, P ert(g) = aj^iq) (l + aia^(g) + a 2 a^-\ ) . 



ai is defined in Eq. ( |4.62|) and 
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C(3) = 1.2020569 . . . denotes the Riemann (^-function. For rif = 5 the numerical values 
are, ai = 0.3802034 . . ., a 2 = 0.9868211 . . .. Since a 2 3> a±, perturbation theory seems 
to be rather slowly or badly convergent. 

Naively, one might expect the perturbative calculation of the static QCD potential 
to be reliable in the limit of large energy scales, q ~ 1/r, i.e at short distances. However, 
unlike in QED, the QCD potential is the ground state energy of a bound state composed 
of the two static colour sources and gluons. Bound state properties are associated with 
a characteristic scale, A, which plays the role of an inverse gluonic coherence length. 
We identify A with the gap between ground state and first excitation. As we shall 
see in Section 5.2 , for large r this gap corresponds to the difference between a hybrid 



state and the ground state potential which, from the bosonic string picture, we expect 
to decrease at large distances like n/r. However, in the limit r — > the gap will not 
diverge but saturate at a constant level that corresponds to the scalar glueball mass, 
A ~ 1.7 GeV. Note that in QCD with light sea quarks it will be even smaller, of the 
order of the mass of two pions. The presence of such a low energy scale, affecting the 
short distance behaviour, can result in differences between the perturbatively calculated 
singlet potential and the static potential. 



16 The lead ing log contribution to the three loop result has been derived by Brambilla and collab- 
orators [22C ] and confirmed in R ef. [ 256 1 while a two loop result for the case of massive quarks has 



recently been obtained by Melles |257| 



47 



Non-perturbative A 4 / g 4 power corrections to ay that are due to the gluon condensate 
are indeed expected from the standard operator product expansion [|261|, [262|, |263|, |264 



Recently, this picture has been challenged by several authors | |265| , |266| , |267| , |268| , |269| , 
|270| , [271| , |272| , |273| , |274|| who found various arguments in support of a term, proportional 

to A 2 /q 2 , 

A-, 



a v (q) = a VtPeTt (q) + c v - 



MS 

q 2 



+ 



(4.73) 



Lorentz invariance implies that no power law corrections of even lower order in 1/q 
exist. For the quenched case, where precise data exist down to lattice spacings as small 
as a" 1 ~ 5.5 GeV, the lattice potential has been compared to perturbation theory and, 



indeed, a non- vanishing value, cy = 4.8 ± 1.4, has been found in Ref. [[275|1 for rif = 0, 
after subtracting one loop perturbation theory. We briefly summarise this analysis 
below. 

A Fourier transform of the momentum space potential yields, 



V(r) 



-a 



oc R {l/r) 



(4.74) 



with [EM 260, 277, 275 
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(4.75) 



where fi = exp(— 7_e)/V and 7# = 0.5772156. . . denotes the Euler constant. While in 
the ultra-violet the effect of the 1/q 2 power correction to ay on V(q) is suppressed by 
a factor, 1/g 4 , this suppression is proportional to r only in position space [Eqs. ( f4.74| ) - 

By employing a recursive lattice finite size technique ||278|, |142||, the ALPHA col- 



laboration has recently obtained a value for the running coupling in quenched SU(3) 
QCD [ |103| ]. They quote the result, 



A^ = 0.602(48)/r , 



(4.76) 



in units of the Sommer scale r . We use this result as an input to determine ajjg at 
a high energy scale (1000 r^ 1 ) and run it down from there to scales using the four 
loop renormalisation group equation, Eq. (|3.22|) . Subsequently, the resulting ajjg(fi) is 
converted into aR(e 7B /i) to one and two loops via Eqs. ( [4.71|) and ( [4.75|) (with cy = 0), 
and the perturbative potential is determined. In Figure |4.6| , we compare the result to 
lattice data. The disagreement is increased when going to higher order perturbation the- 
ory. Furthermore, the difference between perturbation theory and the non-perturbative 
determination is consistent with a linear term ||275|| , as expected from Eqs. (^4.74|) and 

(EH)- 

Because of the significant size of the coefficient, ct2, different ways of re-summing the 
series or performing the Fourier transform can result in somewhat different results and, 
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Figure 4.6: The quenched lattice potential and perturbation theory. The lattice data 
have been normalised such that V(ro) = 0. The values 0.94 Tq 1 and 0.77 r^ 1 have been 
subtracted from the one and two loop formulae, respectively, to allow for comparison. 



therefore, in a different coefficient of the linear term, — 2CpCyAj^ = (3.4 ± 1.0)cr. Due 
to the non-convergent character of the perturbative series, subtracting the perturbative 
tail from a physical operator is never a well-defined procedure anyway. The ambiguity 
involved is related to renormalons that result from the interplay between perturbation 
theory and non-perturbative contributions — between the ultra-violet and the infra-red. 
By subtracting two loop perturbation theory, we find the linear term to have a slope 
about six times as big as the string tension. In contrast, tree level perturbation theory, 
with the coupling being treated as a free parameter, is compatible with the Cornell 
potential, i.e. results in the same linear slope, a, for small and large distances. 



4.9 String breaking 

Having discussed the potential at very short distances, we shall re-examine the large 
distance behaviour. While a linear rise is expected in pure SU (N) gauge theory, in full 
QCD the coupling of gluons to fundamental matter fields will result in a screening of 
inter-quark forces at large distances. However, this behaviour has not been detected 



so far in simulations involving sea quarks (cf. Figure |4.4|) . One reason might be that 
smeared Wilson loops are highly optimised to achieve enhanced overlap with the low- 
est lying string state and might, therefore, almost completely decouple ||255|| from the 



physical ground state at large r that consists of two disjoint static-light mesons. Ar- 
guments based on the strong coupling expansion | |279| , |280| as well as on the bosonic 
string picture |[281|| support this suggestion. Investigating string breaking in full QCD 
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is computationally very expensive as high statistics are required to resolve the potential 
in the large r region of interest. 

While string breaking has not been detected in the Wilson loop, the finite temper- 
ature potential, extracted from Polyakov line correlators at temperatures close to the 
deconfinement phase transition, exhibits a flattening, once sea quarks are included into 
the action ||282| , |283|| . First indications of this effect have been reported as early as 
in 1988 ||284|1 . Unlike Wilson loops, Polyakov line correlators automatically have non- 
vanishing overlap with any excitation, containing static quark and anti-quark, separated 
by a distance, r; in particular the static quarks can be accompanied by two disjoint sea 
quark loops, encircling the temporal boundaries, while in the Wilson loop case, co- 
propagating sea quarks are terminated by the spatial transporters at X4 = and 24 = t. 
The difference becomes perhaps most obvious in the loop expansion of Ref . |[281|| . Since 
for non-zero temperatures the finite temperature potential (or free energy) extracted 
from Polyakov line correlators will in general differ from the potential, extracted from 
Wilson loops, the situation is not yet satisfying. 





Figure 4.7: Graphical representation of the correlation matrix that is relevant for an 
investigation of string breaking in full QCD. 



The first ambitious studies of string breaking in QCD using operators with better 
projection on the broken string state are at present being performed (see e.g. [ 254 , 
286| , |287| , |248|| ). Such calculations involve diagonalisation of a two by two correlation 
matrix between string states and two pairs of static-light states. This matrix is visu- 
alised in Figure £0], where straight lines correspond to gauge transporters while curved 
lines represent light quark propagators that are obtained by inverting the fermionic ma- 
trix, M, of Eq. ( 3.11 ). The off-diagonal elements encode transitions between a string 
state with fixed ends and two static-light mesons while the bottom-right element repre- 
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sents two interacting static-light mesons, separated by a distance r. The second of the 
two contributions to this element corresponds to the exchange of a light meson. The 
situation becomes slightly more involved when the Dirac spin structure is taken into 
account. 

Light quark propagators are normally just calculated for one source point to reduce 
the effort in terms of computer time and memory to a tolerable size. However, pairs 
of quark propagators emanating from different sites are required to calculate the bot- 
tom right element of the correlation matrix for various distances and times. Moreover, 
for a precise determination of expectation values of Wilson loops one usually exploits 
self-averaging: the average of the Wilson loop is not only taken over the Monte-Carlo 
generated ensemble of gauge configurations but also within each configuration; Wilson 
loops with different corner point coordinates are averaged, exploiting translational in- 
variance. This practice is essential to reduce statistical fluctuations to an acceptable 
level. By use of refined stochastic estimator techniques ||288| | to calculate the all-to-all 
light quark propagators required for this purpose one will eventually be able to confirm 
string breaking at distances r rs 2.3 Tq. 
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Figure 4.8: The breaking of the adjoint string in 2 + 1-dimensional SU(2) gauge theory 
at a gauge coupling, (3 = 12, in lattice units (from Ref. [|151|| ). 



In addition to the potential, breaking of closed strings (torelons) in QCD with sea 
quarks has been investigated ||167| , |255| , [248|| . Although the results suggest an effect, its 
statistical significance of 2.5 standard deviations is not yet entirely convincing. While 
the situation in the case of interest is not settled yet, toy models have been investigated 
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in three and four space-time dimensions. Similar to the situation of fundamental QCD 
colour sources being screened by sea quarks, one expects the string between adjoint 
sources in pure gauge theories to decay into a pair of gluelumps (or glueballinos) , bound 
states between a static adjoint source (that can be thought to approximate a heavy 
gluino) and gluons (for a detailed discussion see e.g. Ref. |p89jl ). Until recently, the 
situation was controversial: while in some early studies indications of the breaking 
of the adjoint string in four-dimensional SU(2) and SU(3) gauge theories have been 
reported |290| , |291 |, in a later simulation of three-dimensional SU(2) the broken string 



state has not been detected within Wilson loop correlators [ [292|| . 
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Figure 4.9: The breaking of the fundamental string in the 3 + 1-dimensional SU(2)- 
Higgs model in physical units (from Ref. ||293|| ). fi corresponds to the mass of a Higgs 
particle, bound to the static source. 



In other studies a correlation matrix similar to that of Figure fO] between the string 
state and a two gluelump basis has been investigated. This was first done in four- 
dimensional SU(3) gauge theory ||294|| , followed by simulations of three-dimensional [[29 5| . 
151 and four-dimensional ||296| , |297|| SU (2) gauge theories. The main result of the d = 3 
study of Ref. ||151|| is depicted in Figure |4~8|. In addition to the ground state the first 



three excitations are included into the Figure. At small r, these resemble the first radial 
excitation of the string, two ground state gluelumps and one ground state gluelump 
plus an excited state gluelump, respectively. The horizontal lines indicate masses of 
pairs of isolated gluelumps where Eg stands for the ground state and E' G for the first 
excitation. A similar breaking pattern of the adjoint string has also been confirmed in 
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four dimensional SU(2) gauge theory ||296| , [297|| . Studies of the breaking of the adjoint 
string have been preceded by investigations [|298|| of the breaking of the fundamental 
string in SU(2) gauge theory with a Higgs field in the fundamental representation. The 
latest results from such simulations obtained in three ||299|| and four dimensions [|300| , 
293 1 . confirm the expected screening of the static sources 



see Figure [1 

In QCD one would expect a cross-over of the string state and the two static-light 
meson levels similar to the observations within the toy models considered above, even 
in the quenched approximation ||30 1| , |285|| . How then can one distinguish the quenched 



scenario from the un-quenched one where the string is allowed to break? In the quenched 
case, the separate conservation of baryon and anti-baryon numbers implies that an open 
string state creation operator is orthogonal to a creation operator for the two static-light 
state. Therefore, each operator has zero overlap with the respective other state and only 
the assignment of the ground state to a particular operator will become interchanged 
around the would-be string breaking distance at which the two energy levels cross. 
Unlike the behaviour depicted in Figures |4.8| and fO, no energy gap at this distance 



will occur. Both, the string breaking distance and the associated energy gap, which is 
related to a phase shift in the mixing matrix of the un-quenched relevant for 

an understanding of decay rates such as that of the T(4S) or T(5S) mesons into pairs 
of BB mesons. 



4.10 Colour confinement 

We have established the linearly rising potential in pure Yang-Mills gauge theories as a 
numerical fact and have made this behaviour plausible from strong coupling and string 
arguments. However, the dynamical question of how SU (N) gauge theory as the theory 
of asymptotic freedom results in the formation of colour flux tubes with constant energy 
density per unit length remains un-answered. In the past decades, many explanations 
of the confinement mechanism have been proposed, most of which share the feature that 
topological excitations of the vacuum play a major role. These pictures include, among 
others, the dual superconductor scenario of confinement | |302| , |303j , |304| | and the centre 
vortex model ||305| , |211| , |306| , [307], |308| , |309|| . Depending on the underlying scenario, the 



excitations giving rise to confinement are thought to be magnetic monopoles, instantons, 
dyons, centre vortices, etc.. Different ideas are not necessarily exclusive. For instance, 
all fore-mentioned excitations are found to be correlated with each other in numerical 
as well as in some analytical studies, such that at present it seems to be rather a matter 
of personal preference which one to consider as more fundamental. 



Recently, the centre vortex model has enjoyed renewed attention 1 310 ]. In this 
picture, excitations that can be classified in accord with the centre group provide the 
disorder required to produce an area law of the Wegner- Wilson loop and, therefore, 
confinement. One striking feature is that — unlike monopole currents — centre vortices 
form gauge invariant two-dimensional objects, such that in four space-time dimensions, 
a linking number between a Wegner- Wilson loop and a centre vortex can unambiguously 
be defined, providing a geometric interpretation of the confinement mechanism |3 1 1[ | . 
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Figure 4.10: Electric field distribution between two static SU(2) sources in the MA 
projection, in lattice units, a « 0.081 fm. The sources are located at the coordinates 
(-7.5, 0)a and (7.5, 0)a. 



We will only discuss the superconductor picture, which is based on the concept 
of electro-magnetic duality after an Abelian gauge projection that has originally been 
proposed by 't Hooft and Mandelstam ||302| , |303|| . The QCD vacuum is thought to behave 



analogously to an electrodynamic superconductor but with the roles of electric and 
magnetic fields being interchanged: a condensate of magnetic monopoles expels electric 
fields from the vacuum. If one now puts electric charge and anti-charge into this medium, 
the electric flux that forms between them will be squeezed into a thin, eventually string- 
like, Abrikosov-Nielsen-Olesen vortex which results in linear confinement. 

In all quantum field theories in which confinement has been proven, namely in com- 
pact U(l) gauge theory, the Georgi-Glashow model and N = 2 SUSY Yang-Mills theo- 
ries, this scenario is indeed realised. However, before one can apply this simple picture 
to QCD or SU(N) gluodynamics one has to identify the relevant dynamical variables: 
it is not straight forward to generalise the electro- magnetic duality of a U(l) gauge 
theory to SU(N) where gluons carry colour charges. How can one define electric fields 
and dual fields in a gauge invariant way? 

In the Georgi-Glashow model, the 5*0(3) gauge symmetry is broken down to a 
residual U(l) symmetry as the vacuum expectation value of the Higgs field becomes 
finite. It is currently unknown whether QCD provides a similar mechanism and various 
reductions of the SU(N) symmetry have been conjectured. In this spirit, it has been 
proposed | |304| to identify the monopoles in a [/(l)^ -1 Cartan subgroup of the SU(N) 
gauge theory after gauge fixing with respect to the off-diagonal SU (N) /^(l)^ -1 degrees 
of freedom. After such an Abelian gauge fixing QCD can be regarded as a theory of 
interacting photons, monopoles and matter fields (i.e. off-diagonal gluons and quarks). 
One might assume that the off-diagonal gluons do not affect long range interactions. 
This conjecture is known as Abelian dominance ||312|| . Abelian as well as monopole 
dominance are qualitatively realised in lattice studies of 577(2) gauge theory ||313| , pi If 



in the maximally Abelian (MA) gauge projection pi5|, which appears to be a suitable 
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Figure 4.11: Longitudinal electric field, E = E\\ = |E|, and density of superconducting 
magnetic monopoles, /, in the centre plane between the sources. 



gauge fixing condition. 

In Figure |4.10| , the electric field distribution between SU{2) quarks, separated by 
a distance, r = 15a ~ 1.2 fm, is displayed ||316| . This distribution has been obtained 
within the MA gauge projection. The physical scale, a ~ 0.081 fm, derived from 
the value, \fa = 440 MeV, for the string tension, is intended to serve as a guide 
to what one might expect in "real" QCD. Clearly, an elongated Abrikosov-Nielsen- 
Olesen vortex forms between the charges. In Figure |4.11 , a cross section through the 
centre plane of this vortex is displayed. While the electric field strength decreases 
with the distance from the core, the modulus of the dual Ginsburg-Landau (GL) wave 
function, /, i.e. the density of superconducting magnetic monopoles, decreases towards 
the centre of the vortex where superconductivity breaks down. In this study the values 
A = 0.15(2) fm and £ = 0.25(3) fm have been obtained ||316| , |317|| for penetration depth 
and GL coherence length, respectively. The ratio A/£ = 0.59(13) < l/y/2 classically 
corresponds to a type I superconductor very close to the border of type II behaviour, i.e. 
QCD flux tubes appear to weakly attract each other. However, for a final settlement on 
which side of the Abrikosov limit SU (2) gauge theory lies, quantum corrections should 
be considered. A recent analysis of the same lattice data in terms of the classical four- 
dimensional Abelian Higgs model has resulted in similar conclusions [31S]. For a more 
detailed discussion the reader is referred to Ref. |[317|| . 



55 



5 More static potentials 



We will discuss a variety of excitations of the pure gauge vacuum, such as hybrid po- 
tentials, glueballs, gluelumps, potentials between charges in non-fundamental represen- 
tations and three-body potentials. In particular hybrid potentials, whose short range 
behaviour is related to the glueball and gluelump spectra, turn out to be relevant for 
quarkonia as they give rise to extra states that are not expected from the quark model. 
They are also related to relativistic correction terms to the static potential and deter- 
mine the validity range of the adiabatic approximation as we shall see in Sections |6.3.5 
and |6.5-4 Prior to discussing hybrid potentials, we shall introduce hybrid mesons. 

5.1 Hybrid mesons 



At the same time that QCD was invented it has been noticed |80[ that the spectrum 
of this theory should in principle contain bound states without constituent quark con- 
tent, the so-called glueballs, in addition to the mesons and baryons of the quark model. 
The question, however, arises what constitutes the difference between a flavour singlet 
meson that contains "sea" gluons and a glueball that contains sea quarks. In general 
such hypothetically pure states will mix with each other to yield the observed particle 
spectrum. Still, the possibility of gluonic excitations will result in extra levels within 
certain mass regions that would not have been expected from simplistic pure constituent 
quark model arguments. Moreover, glueballs with exotic, quark model forbidden, quan- 
tum numbers should exist. While in QCD the difference at least between non-exotic 
glueballs and flavour singlet mesons is somewhat obscured, the quenched approximation 
contains "pure" glueballs and the spectrum of such states may be used as an input for 



mixing models ||319 | 



Another non-trivial spectroscopic consequence of the QCD vacuum structure are so- 
called hybrid mesons |320| , |38], [321|| , i.e. mesons with "constituent" glue; by considering 



excitations of the glue, mesons can acquire exotic quantum numbers tocQ. There is a 
slight problem with the notion of "constituent" glue. Neither the number of gluons is 
conserved, nor do they have a non-vanishing rest mass. How then can one define the 
difference between "constituent" and "sea" glue? Do not all mesons include a gluonic 
component? Even in the quenched approximation, where a glueball is a perfectly well 
defined object, we cannot easily switch off "sea" gluons to identify hybrids. What a 
"hybrid" is can only be understood within certain models like bag models J36], |38| . 



the strong coupling lattice model |37|, |41| or the flux tube model |E| that distinguish 
between hybrids and standard quark model states. Such models offer extensions of the 
quark model that help in classifying the observed hadron spectrum and can guide lattice 



simulations as well as sum rule calculations. In Section 6.5.4 we shall also see that in the 



framework of a semi-relativistic expansion the classification can be made more precise. 



17 The possibility of mixing with such exotic hybrids as well as four quark (qqqq) molecules in fact 
renders the notion even of a spin-exotic glueball fuzzy in full QCD. 
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In full QCD an operator that is bilinear in the quark fields with given J PC content 
and flavour related quantum numbers isospin, I and -Z3, strangeness, charm and beauty 
will in general couple to all mesonic states within the given channel: in particular QCD 
makes no clear distinction between states with identical quantum numbers such as the 
flavour singlet states 77, rf, r\ c and r\ h and, for instance, pseudo-scalar glueballs: in the 
flavour singlet sector even the notion of a valence quark as opposed to a sea quark is, 
strictly speaking, ill defined. However, the r] c is experimentally clearly distinct from an 
i], containing light constituent quarks; almost no mixing between the would-be pure cc 
and the corresponding pure light quark state of the naive quark model occurs: while 
this is not an exact symmetry the assumption that the number of valence charm quarks 
and anti-quarks are separately conserved is a very good approximation of the physical 
situation. In this sense, one can still assign a flavour to individual constituent quarks. 

We intend to create a meson, i.e. a state containing a quark, q, and an anti-quark, 
<f , with given 

jPC 

assignment. The most general creation operator that is bilinear in 

the quark fields is, 

X 

where we have chosen the coordinates such that the quark is at the origin, a, (3 = 1,2,3 
are colour and fi, v = 1,...,4 Dirac indices. While T determines the internal spin 
symmetry of the state, the function of the gauge fields, O, generates both, relative 
angular momentum of the quarks as well as excitations of the gluonic degrees of freedom. 
For trivial O^P = 8 a p5(x.), qj^q creates a pseudo-scalar meson (J PC = CT + ) and q^y^q a 
vector (1 ): in a colour and flavour singlet state without relative angular momentum 
of the quarks, Fermi statistics implies, P = — 1 and C = (— l) s (in this special case 
J = S). When one allows q to differ from q', the creation operator is no longer a charge 
eigenstate and scalars q'q (0 + ) as well as axial- vectors q'j^j^q (1 + ) can easily be created 
too. 

In general, each O x is a combination of gauge connections between and x. If one 
allows O to have a non-trivial spatial distribution, angular momentum, L, can be in- 
troduced. This can be achieved by the choice, O x = Yu 3 (9, </>)f/(x), where t/(x) denotes 
a Schwinger line connecting with x and Yu 3 are the familiar spherical harmonics^. 
By combining such angular excitations with the pseudo-scalar creation operator, all qq 
states, J PC = 0-+ 2~+ . . ., can be created while combination with a vector results 
in, J PC = 0++, 1— , 1++, 2—, 2++, . . .. Note that P = C = (-l) l+s . 

Let us now investigate the case, O^f = 5(x)£>q^. The chromo-magnetic field, B = 
J2 a ^ a T a , transforms like an octet under gauge transformations and is traceless. It has 
the internal quantum numbers of an axial vector, while the electric field E is a 



18 On the lattice the continuum 0(3) rotation group is broken down to the discrete point group, Oh, 
associated with the cubic symmetry plus inversions, and the spherical harmonics will be replaced by 
functions that are designed to project onto irreducible representations of the latter subgroup, rather 
than onto continuum I. The necessary group theory has been worked out in Refs. [322, 323] for glueballs 
and in Ref. [324] for hybrid mesons. Since Oh is a subgroup of 0(3), irreducible representations 
of the point gr oup can be subduced from a spin representation of the continuous group. See also 
Refs. |7| |25). 
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Figure 5.1: Lattice paths, O, with Oh®C quantum numbers, and 7\ + , respectively. 

vector, 1 . Obviously, the combination, q ar y 5 B°' l3 q 13 , in which the quarks couple to a 
colour octet and interact with the magnetic gluon, shares the vector J PC = 1 spin 
assignment with the quark singlet state, q a, Jiq a , while both, q a ^iB^q^ and q a ^q a are 
pseudo-scalars: It appears plausible to assume that the colour singlet operators have a 
better overlap with the physical ground state while the colour octet operators show an 
improved coupling with would-be hybrid excitations. We finally note that £ijk<I a ljB^ q^ 
results in a spin-exotic 1 _+ assignment. 

In general, one will employ a spatially extended creation operator. Two examples of 
such lattice operators, O, that incorporate bended gauge transporters (staples) which 
result in a non-trivial gluonic state are depicted in Figure |5.1| . The first one corresponds 
to a lattice spin content, , while the second one is within the T{~ + representation 
of Oh ® C. The lowest lying continuum spin from which T\ can be subduced is, I = 



1. In combining the above paths with various possible quark bilinears ||324|| , the first 
operator projects onto mesons with J PC = _+ , 1 h , 1 ,2~ + ,... while the second 
operator yields, J PC = +_ , 1 ++ , 2 + ~, . . .. Spin-exotic states have been indicated 
in bold. The lightest spin exotic mesons come out to have J PC = 1~ + in studies of 
both, quenched QCD and QCD with two flavours of sea quarks ||326| , |327| , [328]. As a 



next step mixing effects with possible nfi spin-exotic four-quark molecules should be 
considered. 
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Figure 5.2: Creation operator for the E u hybrid potential. 



5.2 Hybrid potentials 

While the distinction between a hybrid meson and an ordinary meson is not well defined, 
a hybrid potential with quantum numbers other than E+ between static colour sources, 
separated by a distance r is clearly distinct from the ground state potential or its radial 
excitations of the ground state potential. Hybrid potentials can be classified in analogy 
to excitations of homonuclear diatomic molecules ||1 76[ |325|| . The relevant symmetry 
group is -Dooh in the continuum and D ih on a cubic lattice (for on-axis separation of the 
sources). An angular momentum Af about the molecular axis can be assigned to the 
state. In addition, the state might transform evenly (gerade, g) or oddly (ungerade, u) 
under the combined parity of a charge inversion and a reflection about the midpoint 
of the axis, rj. Finally, reflections with respect to a plane that includes the axis can 
be performed. For A = |A?| 7^ such reflections just transform one state within a 
A-doublet into the other: Aj — > — Af. However, for A = 0, the transformation property 
under this reflection gives rise to an extra parity index, a v . Conventionally, the angular 
momentum is labelled by a capital Greek letter, A = 0, 1, 2, 3 . . . = E, n, A, <3> . . .. The 
straight line connection transforms in accord with the representation, E+. In Figure |5~2| 
we have visualised a creation operator for the lattice state, E u , that can be subduced 
from the continuum representation, U u . The fact that staples pointing into positive and 
negative directions are subtracted from each other reflects the spin one nature of the 
state. Note that the combinations of Figure [5J] contain similar elementary paths. The 
necessary group theory and lattice operators have been worked out in Ref. ||231| . 



Lattice results for hybrid potentials have been obtained in SU(2) |P31| , |146| , [24C|| and 
SU(3) H, |330|, 1 



33 1 , |114| , p32j| gauge theories as well as in QCD with two flavours 



of sea quarks ||167| , [248|| . For a recent review, see Ref. ||333|| . Employing the adiabatic 



and non-relativistic approximations for heavy quarks, one can estimate possible hybrid 
charmonia and bottomonia levels by solving the Schrodinger equation with such hybrid 
potentials. The only peculiarity is that the angular momentum, K = L + S s , that 
couples to the spin of the quarks, S = Si + S 2 , to produce the total spin, J = K + S, 
differs from the angular momentum due to the relative motion of the quarks, L. S g 
denotes the spin of the gluonic flux tube whose projection onto the axis is, Af = S g r. 
Thus, (AA|S^|AA) > A(A + 1) and k > A. Within the leading order Born-Oppenheimer 
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Figure 5.3: Hybrid excitations of the static SU(3) potential (from Ref. ||114 



approximation, K and A are conserved, but not L or S 9 . The centrifugal term, 1(1 + 1) 
that appears in the radial Schrodinger equation, Eq. ( |A.3| ), has to be substituted by the 
correct factor |7|, (L 2 ) = k(k + 1) - 2A 2 + (S 2 ). 

Mass estimates of hybrid bottomonia, obtained in this way from hybrid potentials, 
can be found in Refs. [[239| , |331| , [332|| . Like in the case of light mesons the 1 h , + ~ and 
2 H quarkonium spin-exotica, that are governed by the U u potential in the adiabatic 
approximation, turn out to be the lightest ones. Within the quenched, non-relativistic 
and leading order Born-Oppenheimer approximations bottomonia hybrids come out to 
lie only slightly above the BB threshold. To this order in the semi-relativistic expansion, 
which does not yet incorporate spin sensitive terms, the masses of hybrid H , 0~ + , 1 h , 
1 , 1 H , 1 ++ , 2 H and 2~ + states are degenerate. It is clear, however, that for the 
non-exotic hybrids the use of an excited state potential within the Born-Oppenheimer 
approximation is at best dubious. 

In Figure |5.3| , the spectrum of hybrid potentials from the most comprehensive study 
so far | |114| | is displayed. Continuum limit extrapolated lattice results are indicated by 
pairs of solid curves while dotted curves correspond to the classical Nambu-Goto string 
expectation in four dimensions, Eq. Q4.27Q . Dashed curves indicate rni/r gaps, added to 
the ground state potential, the leading order contribution of the bosonic string picture. 
To guide the eye, the lowest lying states, and Il n , are included into both plots. 
Note that a interpretation of the 11^ state cannot be excluded from the lattice data. 



However, as we shall see at the end of Section 5.3, other evidence speaks in favour 



of the assignment. Most states are in clear disagreement with the simple model 
expectation up to distances as big as 3 fm where sub-leading terms of the string picture 
are rather small as the differences between dashed and dotted curves show. While this 
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contrasts the findings of Ref. ||213|| for closed strings (torelons) and those of Ref. ||146|| for 
hybrid potentials, investigations of the ground state flux tubes between static sources 
indicate half widths of about 1 fm | |242| , |317| ] . Thus, although 3 fm is big in comparison 
with typical hadronic scales, the amplitude of string fluctuations is still quite large in 
relation to the longitudinal extent. Therefore, in an effective string representation the 
possibility of higher dimensional correction terms to the Nambu-Goto action might have 
to be considered. 

The small distance behaviour exhibits a rich structure too and some states appear 
to try to become degenerate. In particular the change of curvature of the U g potential 
at small r appears puzzling. In the limit, r — > 0, the quarks combine to an octet or a 
singlet colour representation. The octet channel in which the sources explicitly couple 
to gluons should have relevance for the hybrid potentials that differ from the ground 
state by excitations of the gluonic flux tube. One might therefore assume that the 



short distance behaviour |p34], |335|| is determined by the perturbative octet potential 



V a (r) = —l/8V s (r). This is in agreement with the observation that the curvature of 
all potentials (with the exception of Il 9 ) is smaller than and opposite in sign to the 
one of the ground state potential. Note that in the framework of potential NRQCD 
(pNRQCD), the hybrid potentials have also been predicted to follow V Q to leading 
order, up to non-perturbative constants ||335|| . 

We would like to mention that in QED potentials can be classified in exactly the 
same way. Nonetheless, in the deconfined phase, that is realised in nature, the spectrum 
of excitations above the ground state Coulomb potential is continuous since photons of 
arbitrary momentum can be emitted. This is not so in QCD. However, the spectrum of 
QCD potentials will become continuous too above glueball pair0 radiation thresholds 
or, when allowing for light sea quarks, meson pair radiation thresholds. 

5.3 Glueballs, glueballinos and hybrid potentials 

In the limit, r — > 0, the cylindrical symmetry of a (hybrid) potential creation operator 
is enlarged to that of the full rotational group in three dimensions, D^h C 0(3) <8> C 
(or, on a cubic lattice, D ih C Oh <S> C). Irreducible representations of the subgroup 
with spin A can be subduced from irreducible representations of the rotational group 



with spin J > A, as illustrated in Table |5.1| . Note that P = cr v , C = i]a v . Moreover, 



states can be classified as singlets and octets in accord with their local gauge trans- 
formation properties. While a singlet state decouples from the temporal transporters 
within an r = 'Wilson loop" , an octet state couples to a temporal Schwinger line in 
the adjoint representation. In the infinite mass limit, where spin can be neglected, the 
temporal transporter can be interpreted as the propagator of a static gluino, in analogy 
to fundamental lines representing a static quark. Consequently, the octet state is called 
a glueballino or gluelump ||291| , |289| , [294j , |334|| while the singlet state that, neglecting 



quark pair creation, contains nothing but glue represents a glueball. 



19 



Due to momentum conservation radiation of a single glueball is forbidden. 
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Figure 5.4: The gluelump correlation function, Eq. ( |5.2[ ). 



Gluelump masses can be extracted from the decay of the correlation function, 



in Euclidean time. U^{t) denotes an adjoint Schwinger line connecting the point (x, 0) 
with (x, t) and if is a local operator in the adjoint representation. The simplest example 
is, H a oc 5f, where 2tr(H F T a ) = £ fe # fe 2 tr (T b T a ) = H a . This operator corresponds 
to an axial- vector, J PC = from which the -Doo/i. hybrid potentials, U u and £~, can 
be subduced in the limit, r — > 0. The three possible orthogonal choices of the direction 
i of Bf correspond to the dimensionality, 2 J + 1, of the J = 1 representation which is 
identical to the sum of dimensions of the subduced representations, IT U and £~: 2 + 1. 
From Eqs. ( |B.16|) , (|3.7|) and (|3~8|), we obtain lattice definitions of magnetic and electric 
field strength operators, 



that approximate the continuum limit up to 0(a) lattice artefacts [0(a 2 ) in SU (2) gauge 
theory]. In SU(3) gauge theory one would preferably modify the above definitions, 




(5.2) 



(5.3) 



B x ,i — > B' xi = B x>i — Tr(B Xii )l, E x>i — > E' xi = E Xji — Tr(E Xji )l, 



(5.4) 




(5.5) 



(5.6) 
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Table 5.1: What 0(3) ® C representation contains what representations? 









++ 


,1" 


" 2++,3- ... 




0~ 


,1 + 


+ 2~ 3++,... 


K 


o+- 


,1- 


+ ,2+-, 3-+,... 




o-+ 


,1 + 


-,2-+,3+- ... 




1++ 


,1" 


" 2++,2— 


n u 


1+- 


,1" 


+ 2+" 2"+ ... 






,2" 


-,3++,3— 




2+- 


,2" 


+ ,3+-, 3"+,... 



is defined at half-integer values of the lattice time, t/a. Note that IL^y = Xr x 



The correlation function of Eq. ( |5.2j ) is visualised in Figure |5.4| . 1 states can 



be created by operators, Hi oc Ef, or by the operators, Hi oc e^D^B^. The latter 
operator is local in time and would preferably be used in lattice simulations. The five 
operators, DfBf — ^6ijDfB^, couple to 2 states etc.. A table containing continuum 
creation operators for various quantum numbers can be found for instance in Ref. ||335 



The correlation function, Eq. (|5.2|) , can be rewritten in terms of operators in the 
fundamental representation by use of the completeness relation, 

2 T®pT?j 5 = 5 a sS/3 7 - —S a/3 6 7 s, (5.7) 

a 

and the identity (U A ) ab = 2tr (UT a U^T b ). The result reads, 

C(t) = (Tt f \H*U (t)H$rf (tj\ ) , (5.8) 



where the disconnected part, — (Tt H F Ti H F, ^\ vanishes due to, Tr H F = H a TiT a = 0. 
The above correlation function resembles a "hybrid" Wilson loop in the limit, r — > 0. 
In this limit, the Wilson loop can be factorised into singlet and octet components, 

(W 9 (r, £)> = Cl e- m gluelump(«)* + C2e -[ m glu e ball+V E +(na)]t + ( r _> Q), (5.9) 

where on the lattice, V s +(0, a) = V^+(0) + V Be if(a) =0. At r ^> a, V^,+ (r) will approach 
the continuum potential. From the above representation we expect certain groups of 
hybrid potentials to become degenerate with each other as r — ► and to assume the 
mass of the lightest glueball or gluelump within the sector of allowed J PC quantum 
numbers that have overlap with the hybrid string creation operator, W . 

Like static potentials, any gluelump mass will contain a finite contribution and a 
(J PC independent) contribution due to the self-energy of the static sources that will 
diverge in the continuum limit, 

m giueiump (a) = m finite + m se if(a). (5.10) 
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Figure 5.5: The lowest six SU(3) gluelump states, extrapolated to the continuum 



limit EM. 



In order to obtain predictions on (hypothetical) glueballino masses, one has to sub- 
stitute the (unphysical) self-energy by the rest mass of the constituent gluino in some 
appropriate scheme. Keeping this in mind, without additional input, only splittings 
of glueballino masses with respect to the ground state can be determined from lattice 



simulations of glueballino correlation functions. In analogy to Eq. (|4.53| ), we obtain the 
tree level result, 



C A V se ]i{a) 
m sel{ {a) = — 



N 2 - 1 



■V se u(a) > Keif (a) 



(5.11) 



the self-energy associated with the adjoint static source diverges faster than that of the 
two fundamental sources within the static potential. In view of this observation, it is 
clear that in the continuum limit, the glueball within Eq. ( |5.10|) will be the lighter state 
and that the level ordering of the hybrid potentials at zero distance will be determined 
by the glueball spectrum. Increasing the separation a bit such that breaking of the 
rotational symmetry still remains small the generalised Wilson loop will contain a con- 
tribution which resembles the correlation function of a gluelump with a self-energy that 
is reduced as the adjoint source becomes smeared out into two fundamental sources. De- 
pending on the size of gluelump level splittings in relation with the glueball spectrum, it 
is therefore quite possible that at small distances the spectrum of hybrid potentials will 
be guided by the ordering of gluelump levels before, towards r — > 0, the glueballs finally 
take over. Note that if we allow for sea quarks, flavour singlet mesons and meson pairs 
will become lighter than the respective glueball levels and determine the short distance 
behaviour. 
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Figure 5.6: The glueball spectrum of 577(3) gauge theory (from Ref. [|116 



What ordering of gluelump and glueball states do we expect? In the MIT bag 
model P7| for instance the lightest gluonic mode is the TE mode ||32U| | (J p = 1 + ), fol- 
lowed by the TM mode (J p = 1—). Hence, one might expect the axial- vector gluelump 
to be lighter than the vector gluelump. Such concepts have been generalised [|336f by 
assuming that masses of particles increase with the lowest possible dimension of an 
operator with which the state in question can be created. While derivatives, D, have 
dimension m, quark creation operators, q, carry dimension m 3//2 and chromo electro- 
magnetic field operators, E and B, dimension m 2 . Only the 1 H and the 1 gluelumps 
can be created by operators of dimension two; all other states require derivatives or ad- 
ditional fields. Based on this simple picture, one would expect l + ~ to be the lighter 
state since a magnetic operator (B) excites a TE field. The next state would be 1 
(E), followed by 2 (DB) and 2 + ~ (DE) and eventually states containing two deriva- 
tives (3 +_ ,3 ) or two gluonic fields (0 ++ , 2 ++ , CT + , 2~ + ) etc.. Indeed, the gluelump 
spectrum ||334|| of Figure [575] seems to follow this qualitative pattern that has also been 



predicted in Refs. [289, 335 1 



The lowest dimensional operator that can be used to create a glueball has dimension 
four. Here, we would expect the lowest states to be made up from two TE gluons 



65 



(BB), coupling to ++ and 2 ++ , followed by h and 2~ + , containing a TE plus a 
TM excitation, followed by 1 ++ and 3 ++ from dimension five BDB operators (or an 
excited ++ from two TE modes) etc.. However, as is revealed by Figure [5.6| this 
simple picture fails after the first 3-4 states: the l + ~ is too light. The strong coupling 
model [41], in which one would expect the ordering ++ , 2 ++ , l +_ from the perimeter 
of the minimal loop required to create the state in question on the lattice, in contrast, 
fails to predict the low mass of the pseudo-scalar glueball. Of course an abundance 
of alternative qualitative and quantitative pictures of the QCD vacuum exists that 
result in somewhat different expectations. A detailed discussion of such models and the 
underlying assumptions is beyond the scope of the present article. 

From the spectrum of glueballs^] and Table [54] we expect the E+' potential to be 
separated from the ground state by a scalar glueball mass m(0 ++ ) at small distances, 
followed by three degenerate potentials II 9 and A g which will be separated from the 
ground state by m(2 ++ ), £~ separated by m(0 _+ ), another m(0 ++/ ) triplet of potentials 
and a set of II n and £~' states, separated by m(l + ~). In the regime of somewhat bigger 
r, which is dominated by gluelumps, we expect a low, almost degenerate pair of hybrid 
potentials, U u and S~, corresponding to followed by a H g , T,' (1 ) pair and a 
S 9 ,n^, A g (2 ) triplet. Indeed, Figure reveals that the S~ and U u potentials are 
the lowest excitations at small r, and approaching each other. With r — > we would 
expect the levels to cross as the value of S~ will tend towards the ground state potential 
plus a pseudo-scalar glueball mass. Confirmation of this effect, however, requires lattice 
spacings that are sufficiently small to yield a gluelump mass exceeding that of the 
glueball in question^. All the remaining levels are in complete agreement with the 
ordering and degeneracy expectations from the gluelump considerations too, with the 
exception of the A u that comes out to be somewhat higher than its degenerate 2 H 
and IT4 partners. Unfortunately, no data on U' g exists, which we would have expected 
to become degenerate with A g and S~ at small r. 

Lattice simulations ||334|| reveal that at spacings, a -1 > 2 GeV, the sum of the scalar 
glueball mass and the ground state potential at the shortest accessible distance, V s +(a), 
becomes smaller than the mass of the lightest (l + ~) gluelump. In the framework of ef- 
fective field theories (see Section a cut-off on gluon momenta is imposed. We conclude 
that as long as this cut-off does not exceed about 2 GeV hybrid related interactions are 
governed by the spectrum of gluelumps at short distance while when allowing for harder 
gluons, glueball channels will become increasingly important. 



20 When allowing for light sea quarks, due to mixing with flavour singlet mesons, the level ordering 
will be completely different, starting with the pseudo-scalar E~. 

21 Moreover, some hybrid Wilson loops are constructed in such a way that one would expect them to 
better project onto states determined by the gluelump spectrum rather than the glueballs which will 
complicate numerical studies of the expected level crossings at short distance. 
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5.4 Casimir scaling 



It is possible to determine the potential between colour sources not only in the funda- 
mental representation (quarks) but in any representations of the gauge group. We have 
already discussed bound states between static adjoint sources (gluinos) and relativistic 
gluons above. Despite the availability of a wealth of information on fundamental poten- 
tials, only few lattice investigations of forces between sources in higher representations 
of gauge groups, SU(N), exist. Most of these studies have been performed in SU{2) 
gauge theory in three JJ37], ^ |25§, [TH] and four |S3B|, [33g [SI], [230} Q f25§, P7f| 



space-time dimensions. Zero temperature results for SU(3) can be found in Refs. [ 291 
294| , |333| , |342| , |343j , |344| while four-dimensional determinations of Polyakov line corre- 
lators in non-fundamental representation have been performed at finite temperature by 
Bernard |38], |45| for SU{2) and Refs. |H |47|, |4|, |8| for SU{3) gauge theory. We 
have already discussed the SU(2) results of Refs. ||295| , |151j , |296|| and the finite temper- 
ature results in the context of string breaking in Section 4J3 and shall focus on d = 4 
SU (3) zero temperature simulations below. 



Table 5.2: Group factors for SU(3). D is the dimension of the representation, (p, q) are 



the weight factors, z 



exp(27rz/3 
D 



6 
15a 
10 
27 
24 
15s 



and dr, = Cd/Cf denotes ratios of Casimir factors. 



p,q 



1,0 

1,1 
2,0 

2,1 
3,0 
2,2 

3,1 
4,0 



yP-q 



z 
1 

z* 
z 
1 
1 

z* 
z 



p + q 



1 

2 
2 
3 
3 
4 
4 
4 



d 



D 



1 

2.25 

2.5 

4 

4.5 
6 

6.25 
7 



For the static potential in the singlet channel in position space, tree level perturba- 
tion theory yields the resultQ 



V(r,fi) = -C D hoWseifC"), 

r 



(5.12) 



in analogy to Eq. ( f4.45| ). D = 1,3,6,8,10,... labels the representation of SU(3). 
D = 3 corresponds to the fundamental representation, F, and D = 8 to the adjoint 
representation, A. Cd labels the corresponding quadratic Casimir operator, Cp = 
Tr r,T®T® , with the generators fulfilling the commutation relations of Eq. ( B.7 ), 
[T<P , T fe D ] = ifabcXc- Table |5~^ contains all representations D, the corresponding weights 
(p, q) for p + g < 4 and the ratios of Casimir factors, (io = Cd/Cf- In .Sf/^) we have 
CV = 4/3, and z = exp(27ri/3) denotes a third root of 1. 



2 In fact this relation turns out to hold to at least two loops (order a^) |25£ 
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We denote group elements in the fundamental representation by U. The traces of 
U in various representations, Wp> = tiUp,, can easily be worked out, 



W 3 
W 8 

w 15a 

W w 

w 27 
w 15s 



trU, 



(|^3| 2 -1 



- [(tiU) 2 + tiU 2 

tr U* W ( 
1 



6 



trU, 

(tr U) 3 + 3trUti U 2 + 2trU 3 
We, 



trCTWio 
1 



w 6 \ 2 -\w 3 \ 2 , 



24 

+ 8tr£/tr£/ 3 



(tr U) 4 + 6(tr U) 2 tr U 2 + 3(tr U 2 ) 2 
6 tr U 4] 



5.13) 
5.14) 

5.15) 

5.16) 

5.17) 

5.18) 
5.19) 

5.20) 



l> - TjW d- the normalisation of Wp> differs from 



Note the difference, Tt d Ud = 4jtr JJ D = -kW, 
that of the Wilson loop of Eq. ( |4.1| ) by a factor D. Under the replacement, U — > zU, 
W D transforms like, W D — > z p ~ q W D . 

In Section [4.7.2| we have seen that for distances r > 0.6 ro ~ 0.3 fm the fundamental 



potential is well described by the Cornell parametrisation, 



V F (r) = V Q F h o F r. 

r 



(5.21) 



Perturbation theory [Eq. (|5.12j )1 tells us, Vq } d ~ ^dVo,f and en ~ dp,&F- While the 
fundamental potential in pure gauge theories linearly rises ad infinitum, the adjoint 
potential will be screened by gluons and, at sufficiently large distances, decay into 
two disjoint gluelumps. This string breaking has indeed been confirmed in numerical 
studies HI H IH, HH, H • Therefore, strictly speaking, the adjoint string 
tension is zero. In fact, all charges in higher than the fundamental representation will be 
at least partially screened by the background gluons. For instance, 6® 8 = 24© 15a* © 
6 © 3*: in interacting with the glue, the sextet potential obtains a fundamental (3*) 
component. A simple rule, related to the centre of the group, is reflected in Eqs. ( |5.13| ) 
- ( 5.21|) : wherever z p ~ q = 1, the source will be reduced into a singlet component 

it will be screened, up to a residual 



z(z 



at large distances while, wherever z p ~ q 
(anti-)triplet component, i.e. one can easily read off the asymptotic string tension (0 
or ap) from the third column of Table |5.2| , rather than having to multiply and reduce 
representations. As a result, the self adjoint representations, 8 and 27, as well as the 
representation, 10, will be completely screened while in all other representations with 
p + q < 4 a residual fundamental component survives. The same argument, applied 
to SU(2), results in the prediction that all odd-dimensional (bosonic) representations 
are completely screened while all even-dimensional (fermionic) representations will tend 
towards the fundamental string tension at large distances. 
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Figure 5.7: Static potentials between sources in various representations of SU(3) in 
lattice units, a a pa 0.085 fm. 



While the string tension approaches either or op at very large distances, in an in- 
termediate range an approximate linear behaviour is found ||340| , |290| , |291| , |348j , 
342 , |295| , |151| , |343|| , such that one might speculate whether in this region the Casimir 
'i'67 \, op, pa dp><jp, that is exact in two-dimensional QCD, holds. 

13, |g lattice 



scaling hypothesis 



This hypothesis has been challenged by the fact that in all but two 
simulations the expected Casimir slope is under-estimated. Motivated by this obser- 
vation other models have been suggested, like scaling in proportion with the number 
of fundamental flux tubes embedded into the higher representation vortexp 3 ! [p + q in 
SU(3)] ||333| , 513 1, which happens to coincide with Casimir scaling in the large N limit. 
Casimir scaling and flux counting predictions, at least for the lower dimensional repre- 
sentations, are close to each other, such that discriminating between them represents a 
numerical challenge. 

23 Some other alternatives have been suggested in the past. In a bag model calculation, for instance, 
the result, a D = \fCuJCp ap, was obtained |34£ 
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The latest lattice results for SU(3) gauge theory from Ref. ||343|| are displayed in 



Figure |5J in lattice units, a a m r /6 ~ 0.085 fm. Note, that the raw lattice data 
are displayed and no self-energy pieces have been subtracted. The data have been 
obtained on lattices with an anisotropic Wilson action and tiny temporal lattice spacing, 
a" 1 « 4a~ 1 « 24 r^ 1 « 9.5 GeV. The fundamental potential for distances, r > 0.6 r , 
has been fitted to Eq. ( |5.21|) . The expectations on the potentials, Vo(r), which are 
displayed in the Figure, correspond to the resulting fit curve, multiplied by the factors 
do- As one can see, up to distances where the signal disappears into noise or the string 
might break, the data are well described by the Casimir scaling assumption. Since 
this study has been performed on the finest lattice resolution so far, it can very well 
be that the underestimation of the Casimir scaling prediction of previous studies is 
a lattice artefact that will disappear after an extrapolation to the continuum limit. 
Indeed, a close inspection of data for the three-dimensional SU{2) case ||292| , [295| , |151 



shows the tendency that the Casimir scaling expectation is approached from below 
with decreasing lattice spacing. It is still an open question whether Casimir scaling 
only holds approximately or if it is exact for distances smaller than the corresponding 
string breaking scales. 

5.5 Three-body potentials 

Although weak decays turn the phenomenology of hadrons composed of three or more 
heavy quarks experimentally unpromising, predicting properties of such heavy quark 
systems is the starting point for understanding multi-quark bound states from QCD and, 
eventually, nuclear physics. The first steps into the latter direction of including light 
quarks have been done by Michael and Pennanen who investigate systems composed of 
two light and two heavy quarks ||285|| or the even more ambitious study of the uuddss 
-ff-dibaryon, containing six light valence quarks by Wetzorke and collaborators ||350| . 



Forces between three and more static sources are not only interesting to guide the 
phenomenology of multi-quark states and to develop and test the lattice methodology 
required in this context but also for model builders: can multi-quark interactions be 
understood in terms of two-body interactions or have genuine three- and many-body 
effects to be considered? Hadronisation models for instance, which intend to explain 
the formation of hadronic jets in high energy scattering experiments, crucially rely on 
a factorisation hypothesis. 

In the past years the Helsinki group has made extensive investigations of systems 
composed of four static SU(2) sources ||351| , p52| , |353| , [354j , |355|| , where no distinction 



between quarks and anti-quarks exists. These systems, therefore, have the capacity to 
approximate both, meson-meson and baryon-baryon interactions in QCD. The lack of 
difference between baryons and mesons is of course a serious limitation when trying 
to understand multi-quark interactions. For instance, unlike in SU(3) where just two 
different pairings within four quark systems are possible, in SU(2) three different ways 
of dividing the system into two colour singlets are viable: combinatorially, the four 
quark system finds its generalisation in SU(3) systems composed of six quarks; however, 
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geometrically, qqqq systems come closer. Here, we will restrict our discussion to the 
simpler case of three quarks in QCD. 




Figure 5.8: A baryonic Wilson loop. 

In analogy to the standard (mesonic) Wilson loop, in SU(N) (N > 3) gauge theories 
gauge invariant baryonic Wilson loops, W^q, can be defined, for the case of SU(3) the 
binding energy of a system of three static quarks at positions x 1; x 2 , and x 3 (baryonic 
potential, V^) can be extracted in the limit, t —>■ oo. The baryonic Wilson loop is com- 
posed of three staples, U\ % = 1,2,3, whose colour indices are contracted at Euclidean 
times and t by completely antisymmetric tensors, 



W 33 (xi,X2,x 3 ;t) 



-ULuLuL. 



(5.22) 



The definition of the staples, U\ is evident from Figure |5.8| . The spatial parts of the 
baryonic loop will in general be composed of fat or smeared links for enhanced overlap 
with the physical ground state. The contraction of the colour indices can take place at a 
spatial coordinate at time t that differs from that at time 0. Moreover, the contraction 
points do not necessarily have to differ from the static quark positions, Xj. In the limit, 
r i2 — ► 0, two quarks combine to an anti-triplet, 3 (g> 3 = 3* © 6, that interacts with 
the remaining quark at position x 3 : in this limit the baryonic Wilson loop becomes a 
mesonic Wilson loop. 

The colour factor that accompanies the tree level perturbative result reads, 



1 

3! 



, , rpa rpa 



_i Tr(:rn = _CV 



(5.23) 
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Figure 5.9: The star geometry with, ry = J2i r iS, an d the A geometry with, r& 

Si>j r ij- 



With another minus sign due to the different relative orientation of the quark lines with 
respect to a mesonic Wilson loop we arrive at the tree levelp^ relation, 



^3 9 (x 1 ,x 2 ,x 3 ) = ~ [V{r n ) + V(r 23 ) + V(r 31 )} , 



(5.24) 



between baryonic and mesonic potentials. For the non-perturbative long range part, 
two models compete with each other, which we shall refer to as the star (or Y) and 
the A laws. The first model originates from strong coupling and area minimisation 
considerations | |356| , |357] , H3L |358| , |359| |. The solution of the problem of finding the 
shortest connecting path is well known for the case of three points, i.e. for a planar 
geometry; three straight lines emanating from the quarks will meet at an angle of 27r/3 



at a central Steiner point, x# (Figure |5.9|) . Unless one of the angles within the baryonic 
triangle exceeds the value 27r/3, in which case a linear geometry will be preferred, the 
resulting minimal area configuration resembles a Mercedes star shape. In this case, we 
expect the baryonic potential to be described by parameters extracted from a Cornell 
fit, Eq. ( |5.21| ), to the mesonic potential, in the following way, 



F 39 (xi,X 2 ,X 3 ) 



2 V °-2^ + arY ' 



(5.25) 



where we have approximated the terms associated with the short range behaviour by 
the perturbative expectation^} 

In fact, this relation turns out to hold at least to order a 2 s in perturbation theory. 
25 0ne might, however, argue that at least at large r the 1/r term is related to Gaussian string 



fluctuations around the minimal area string world sheet [200, |190| and try a somewhat different ansatz 
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Figure 5.10: Three quark potential obtained on an equilateral triangle, rv,- 
r 23 = 7*31, in SU(3) gauge theory with Wilson action at (3 — 6.0 |[360|| . 



ri2 



The competing contender is the A law [|361| , [362 

V3,(xi, x 2 , x 3 ) ss -V - | — + | r A : 



(5.26) 



although obviously, r& > ry, in this case each static quark line is shared by two surfaces 



which, depending on the underlying model ||362|| , can result in a pre-factor, 1/2, to 
avoid over-counting. Since ta/2 < ry, this might then be the dominant configuration. 
Obviously, whenever the three quarks belong to a straight line, the two models yield 
identical predictions. The biggest difference is encountered for the case of an equilateral 
triangle where the predictions disagree by about 15 %, 



r Y 



2_rA 
2 ' 



(5.27) 



Only very few lattice results [|363| , |364| , [365|1 on baryonic potential existed so far, with sta- 
tistical errors too big to rule out either possibility. In an as yet unpublished study ||360|| . 
however, clear evidence in support of the A law has been found. The result for an 
equilateral triangle in lattice units is displayed in Figure |5.10| for SU (3) gauge theory 
with Wilson action at (3 = 6.0 (a ^ 0.094 fm): the data perfectly agree with the sim- 
ple expectation, V^ q = (3/2)V r g g(ri 2 ), which, like the ratios between potentials between 
charges in different representations of the gauge group presented above, happens to 
coincide with tree level (and higher order) perturbation theory, Eq. (|5.24j) . 

It is interesting to notice that phenomenological fits of the baryon spectrum, for 



instance in the framework of relativised quark models [59] in which the Y law is assumed, 
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yield a string tension that is reduced by about 20 % ||366|| , compared to the corresponding 
mesonic result. This is fairly consistent with the lattice results for the configuration of 
the equilateral triangle presented above. We conclude that while the agreement of the 
lattice data with the A law is appealing the question whether the A law holds or an 
approximate Y law with a reduced string tension is satisfied cannot conclusively be 
answered until additional source geometries have been investigated. 



6 Relativistic corrections 

In this Section we attempt to bridge the gap between QCD dynamics of heavy quark 
bound states and potential models. We will sketch the derivation of a quantum me- 
chanical Hamiltonian, containing the static potential as well as semi-relativistic cor- 
rection terms. To leading order this has been pioneered by Wilson, Brown and Weis- 
berger |], |175|| some 20 years ago. As soon as the approach was generalised to higher 



orders |22| |23], [25[ in the inverse heavy quark mass, m~ , or, better, relative heavy 
quark velocity, v, certain inconsistencies appeared between the non-perturbatively de- 
rived general form of the interaction and a direct perturbative evaluation ||367|| of the 
potential between two heavy quark sources at order a 2 s /m 2 . 

A lot of progress in the understanding of effective theories, in particular in the 
matching of low energy theories to QCD has been achieved since then, and the problem 



is now understood [28] and removed. Motivated by these developments, we choose to 
start our discussion from non-relativistic QCD (NRQCD) in the continuum and on the 
lattice, before we address relativistic corrections to the heavy quark potential. Special 
emphasis is put on the matching problem. We shall also see that the validity of the 
adiabatic approximation is very closely tied to that of the non-relativistic expansion. 
Finally, lattice results on the heavy quark interaction will be presented. 

6.1 NRQCD 
6.1.1 The problem 

We wish to consider mesonic bound states that contain two heavy quarks, namely the 
J/ip, T and B c quarkonia families. Typical binding energies, A, turn out to be a 
few hundred MeV, similar to systems that are entirely composed of light constituent 
quarks. The quark mass, m, however, is much larger. This difference in scales results in 
complications when evaluating physical properties. In a standard lattice computation 
for instance one has to work at lattice cut-offs, a <C m _1 , in order to resolve the heavy 
quark while at the same time the box size has to be kept sufficiently large to resolve 
the scales that are relevant for the dynamics of the bound state like the binding energy, 
L a a ^> A . This results in a prohibitively large number of lattice sites that seems 
physically unnecessary since the scale, m 3> A, appears to be rather irrelevant for the 
quarkonia level splittings (cf. Section |2.2| ). Indeed, closer inspection shows that only the 
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Table 6.1: Hard, soft and ultra-soft scale estimates. 





J/1> 


T 


tt 


e + e 


m 


1.4 GeV 


4.7 GeV 


175 GeV 


511 keV 


mv 


0.7 GeV 


1.3 GeV 


45 GeV 


3.7 keV 


mv 2 


0.4 GeV 


0.4 GeV 


12 GeV 


0.027 keV 


V 


0.5 


0.29 


0.26 


0.007 



temporal lattice spacing, a T <C m _1 , is limited by the quark mass. The computational 
effort becomes tolerable when an anisotropy, a T ja a ~ A/m, is introduced. 

The two scale problem can even be turned into a virtue within an effective field theory 
formalism. The strategy would be to integrate out the ultra-violet behaviour at scales, 
jiiRim, into local Wilson coefficients of an effective low energy action that encodes the 
information relevant for bound state properties. Heavy quark effective theory [|368| , |369| . 
|370|| (HQET) for instance is a very effective framework for the calculation of properties 
of systems containing one heavy quark. The strategy is to write down an effective action 
that approximates QCD to a given power v in m~ l . In general, the effective Lagrangian 
will then contain all operators of dimensions smaller than or equal to that of m u+4 . 
The coefficients that accompany these terms can be determined by matching on-shell 
Green functions, calculated in the effective theory, to those calculated in QCD, in the 
ultra-violet. This can be done for example in perturbation theory which is supposed 
to be applicable as long as, fi ^> Aqcd- The tree level matching coefficients can be 
obtained by formally expanding the Dirac Lagrangian in terms of m~ x . 

Although NRQCD [33|, that applies to systems containing two heavy quarks, is 
somewhat more involved it has in fact been formulated earlier than HQET. The power 
counting scheme required for quarkonia differs from the one used in heavy-light systems. 
This is related to the fact that in the lowest order HQET Lagrangian, heavy quarks 
with non-vanishing relative velocity decouple from each other. In order to allow for 
interactions, a kinetic term, p 2 /2m, is required that causes changes of the relative quark 
velocity, v. Therefore, the lowest order effective Lagrangian depends explicitly on the 
quark mass, m, in a way that cannot be absorbed into simple field redefinitions: the 
HQET power counting is obscured and a different expansion parameter is required. 

As an alternative it has been suggested to expand the effective Lagrangian in terms 
of the quark velocity, v. One consequence is that in NRQCD a hierarchy of scales, 
m ^> mv ^> mv 2 ^> . . ., is introduced. The binding energy, A, is of the order of 
the ultra-soft scale, mv 2 , while the typical three-momenta exchanged, mv, are soft. 
The hard scale, m, is integrated out into matching coefficients, Cj(/i/m, a 8 ), at a scale 
m > /i > mv. With the hierarchy of scales comes the possibility of a hierarchy of 
effective theories: after integrating out the soft scale, mv, another effective field theory, 
potential NRQCD (pNRQCD), can be formulated ^TTJ, g35 |. 

In Table |6.1|, estimates of the three scales for the charmonium, bottomonium and 
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Figure 6.1: Determination of the relevant coupling for hypothetical tt bound states. 



(unstable and therefore hypothetical) toponium ground states from a potential calcula- 
tion are listed. For comparison, we include the corresponding estimates for positronium. 
The top quark is so heavy that for our estimate only a Coulomb like potential needs to 
be considered. In this case, V(r) = —Cpan/r, and one easily obtains from the virial 
theorem [Eq. with v = 1], (V) = -2(T) = 2E h using E x = -mC 2 F a 2 R /A, 

WJ. = C *a? R , y—L=-pnC FaR . (6.1) 

In Figure |Q| , 2/(Cj?7«i)(r _1 ) and are plotted as functions of r. The upper rif — 5 
curve has been obtained from the input value, o.jjg{mz) = 0.1214(31), from e + e~ 
experiments at LEP and SLAC J52], by use of Eqs. ( |4.71| ) and ( [4.751) . The quenched, 



rif = 0, estimate has been calculated by use of the lattice result of Ref. |pL03|| , displayed 
in Eq. ( |4.76|) . The intersects correspond to the values, a^(17GeV) ~ 0.145, and, 
a^(23 GeV) ~ 0.20. The estimates quoted in Table [TT] were obtained using the latter 
{rif = 5) result. The matching coefficients between QCD and NRQCD are calculable 
in perturbation theory as long as m 3> Aqcd while perturbative matching between 
NRQCD and pNRQCD can be performed whenever mv ^> Aqcd- However, given the 
numbers in the Table, a reliable perturbative determination of the matching coefficients 
for charmonia states appears to be doubtful while for top quarks even perturbative 
pNRQCD should be applicable, up to power corrections ||274| , |373| , [33| . 
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6.1.2 The NRQCD Lagrangian and power counting 

In order to derive an effective field theory that includes a kinetic term in its leading 
order Lagrangian, we introduce a second dimension [v] = c, in addition to the mass, 
and expand the Lagrangian formally in terms of ||374|, |375|| 1/c. As a result, time is 



measured in different units than space: t = x 4 /c, d t = cd 4 and D t = cD 4 = d t — iA 4 . 
The spatial covariant derivative reads, Di = <9j — -Ai. The kinetic term, tr F^F^, has 
dimensionP] m 4 c 6 . We define, 



S 



-Jd 3 x JdtC, (6.2) 



where, 



C = ^tr F^F^ + q[y ■ Dc + j 4 D t + mc 2 ]q. (6.3) 



1 

In order to make the leading order NRQCD Lagrange density independent of c, we have 
rescaled C by an overall factor, c [Eq. ( |6.2| )], which is compensated in the fermionic part 
by rescaling the fermion fields, q and q, by factors, y/c. 

At tree level a classical derivation of the NRQCD Lagrangian from the Dirac La- 
grangian is possible by means of the Foldy-Wouthuysen-Tani (FWT) rotation [ |376] , |377|1 : 



one starts from the Dirac basis in which 74 is diagonal and then, order by order in 
p/(mc), the fermion kernel is iteratively block diagonalised into separate non-interacting 
quark and anti-quark sectors. The first order transformation takes the form, 

q -> exp (j^^j q, (6-4) 

which happens to be the correct expression to all orders in the free field case. Subse- 
quently, the rest mass can be removed by rescaling the (anti-)fermion fields with the 
factor, exp(±mc 2 t). For details see e.g. Ref. J35|. We choose the decomposition, 




(6.5) 



of the Dirac spinor in the FWT basis into quark and anti-quark Pauli spinors. 

The resulting effective Lagrangian to order c -2 for the two-particle sector reads, 



with 



+ C x = [D 4 + H i> ]iP- X ] [Da - H[]x (6.6) 



D 2 ga-B (D 2 ) 2 
H * = q c F~ n ^—3-^ ( 6 - 7 ) 

. c j(D -E-E-B) | ^-(DAE-EAD) | 1 3 

8m 2 c 2 8m 2 c 2 c 3 



Note that [t] = l/(mc 2 ), [xi\ = l/(mc) and, gE t = -i[D i ,D 4 ] 1 gB l = -ieyfcc [Dj, Dk\- 
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Table 6.2: NRQCD operators and dimension. 



operator 


dimension 


effect 


description 


C 


4 6 

m c 


4 5 

m f 


Lagrange density 


v> 


m 3//2 c 2 


(mt> ) 3 / 2 


quark creation operator 


X 


m 3/2 c 2 


mv 2 


anti-quark creation operator 


A 


mc 2 


covariant time derivative 


D 


mc 


mv 


covariant spatial derivative 


<?E 


i ^ 
m c 


2 3 

m v 


chromo-electric field 


<?B 


2 3 

m c 


2 4 

m v 


chromo-magnetic field 




1 


V 


strong coupling "constant" 



Note that we are using the conventions of Eq. ( |B.16| ) to relate the field strength tensor 
and electric/magnetic fields. The well known Fermi term of the Pauli equation that is 
responsible for the hyperfine splittings in atomic physics is accompanied by a matching 
coefficient, Cp, the Darwin term by Cr> and the spin-orbit (Thomas) term by c$- The 
normalisation within Eq. ( Jb.7]) is such that in the free field limit, Cp — Cr> — Cs — 1. 
The kinetic term defines the mass, m^ = m + 5m, where Sm accounts for the difference 
in the self-energy subtractions between effective theory and QCD. In the classical limit, 
5m = 0. The coefficient of the relativistic correction to the kinetic energy is fixed by 
Lorentz [or, in Euclidean space, 0(4)] symmetry. If this is broken, as for example on 
the lattice, it can also obtain a non-trivial value ||378 . 

In Table |6.2j , we list the naive dimensions of various operators as they result from 
the above equation. By considering the fermionic part of the action and writing down 
the equations of motion in Coulomb gauge, phenomenological scaling laws have been 
derived |3^] that somewhat differ from the power counting in c but should closely resem- 
ble the relative numerical importance of a given operator with respect to the resulting 
quarkonium spectrum^]. According to the analysis of Ref. plj , the coupling g 2 is ex- 
pected to scale in proportion to the velocity [which is the case for a Coulomb potential, 
Eq. (|6.1|)], making the rescaling of the quark fields by factors, yfc, in Eq. ( |6.3| ) super- 
ficial. The results of Ref. [34| are included in the third column of the Table. In what 
follows, we will distinguish between "dimension" in terms of c of a given operator and 
the "effect" on energy levels of a bound state in terms of v. The above order c~ 2 La- 
grangian, without radiative corrections (a s cx v), corresponds to order v 4 in terms of 
these (original) NRQCD power counting rules. 



We find the v power classification of Ref. |34| useful for phenomenological purposes. 
However, we believe that for a consistent construction of an effective field theory, formal 
expansion in terms of a dimensionful parameter, c, is more illuminating ||374| | . As long 



as it is not possible to cleanly disentangle soft 



mv , 



and ultra-soft 



[mv 



degrees of 



27 



To add to the confusion, yet another set of counting rules that arises from a multipolc expansion 



of the gluon field has recently been suggested [379|. 



78 



freedom, each operator will receive additional contributions that are sub-leading in 
v; the velocity scaling arguments are not exact but have to be interpreted at leading 
order. Moreover, the effective v will depend on both, the state under consideration and 
the operator in question. While an expansion in terms of c -1 can be performed on the 
quark-gluon level, the velocity size classification is based on bound state properties. For 
instance, we find the Fermi term, cr-B, of Eq. ( |6.7| ) to be of order c _1 while according to 
the velocity classification it has relative size v 2 . Within a bound state, the spin variable 
within this term has to be saturated by a second spin, such that its leading effect on the 
energy levels is suppressed by an additional power of c _1 , in accord with the velocity 
size counting. Taking such bound state arguments into account, it appears favourable 
to always truncate the Lagrangian at an even order in c _1 . The Darwin and Thomas 
terms have the same dimension (c -2 and v 2 ) in both counting schemes. In what follows, 
we will assign the orders 1 and v 2 in the c _1 and v power counting schemes, respectively, 
to the lowest order Lagrangian. 

For completeness of the effective Lagrangian we have to consider the two-particle 
sector: 

d ■ d 
= !f — oV'V'XXH ^—sW^X^X (6-8) 

m^,m x c A m^,m x c A & 

Note that by means of a Fiertz transformation an alternative basis can be chosen, 

4 X = — ^-s^xxV + — ^-jV'VxxW (6-9) 

+ -^- 2 ^T*xX ] T a ^ + -^- 2 ^T a a X X ] T a aij + i(9(c 7 ). 

m^,m x c z m^m x c z c A 



The coefficients are related to each other [ 5SD | 



d ss = ^[-d c ss -3d c sv -C F d c vs -3C F d c v ], (6.10) 

d sv = ±j[-d c ss + d c sv -C F d c vs + C F d c v ], (6.11) 

d vs = ^[-2Ndt s -6Nd c sv + d c vs + M c v ], (6.12) 

d vv = ±[-2Nd c s3 + 2Nd c sv + d c vs -d c v ]. (6.13) 

For m$ 7^ m x the tree level coefficients are all zero while in the equal mass case the 
coefficients in the standard basis are proportional to a s , due to an annihilation diagram 
that results in, d c vv = —na s H — -. According to the v power counting rules the first two 
terms are only suppressed by a factor v, relative to the leading kinetic term. However, 
the matching coefficients guarantee a suppression by an additional factor of a s ~ v, 
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in the equal mass case and of a 2 ~ v 2 for ^ m x . Therefore, the effect of the 
contact terms is of combined orders v A and v 5 , for equal and non-equal quark flavours, 
respectively. 



Finally, we consider corrections to the gauge action ||381 

h 

2 



g 2 £-YM — —toF^F^ 



~ \^ + ^h]^F^ 2 F^ (6.14) 
~ 2(^ + ^)trF, l [D l ,D J ]F Jll + - 3 0(c 7 ). 

We have adopted the notation, D 2 = D^D^. The tree level values are, &2,v> = ^2, x — 
bz,ip = &3,x = 0- ^ turns out that the radiative corrections to the tree level value, b\ = 1, 

in 



compensate the effect of the heavy quarks on the running of the QCD coupling ||380 1 
the effective theory, between the matching scale, /x, and the quark masses, and m x . 
Heavy quark loops are subsequently explicitly reintroduced into the NRQCD Lagrangian 
via the terms containing derivatives that are proportional to 62 and 63. As long as we 
are only concerned with the quenched approximation to QCD the corrections to the 
gauge action should be ignored. In an un-quenched world they should, however, be 
included for consistency. It is clear though that such heavy-quark un-quenching effects 
are numerically tiny in comparison to the error one makes when ignoring light sea quarks 
for example. 

6.1.3 Matching NRQCD to QCD 

We wish to construct an effective theory which is applicable to gluon momenta, q < /i, 
where mv < \x < m, and which reproduces QCD up to corrections that are of higher 
order in terms of the expansion parameter (m -1 in the case of HQET and c _1 in the case 
of NRQCD). Differences between the effective theory and the correct QCD behaviour 
that would otherwise arise for momenta, q > fi, have to be compensated for by an 
adequate choice of the Wilson coefficients, Cj(/x/m, a s ), dj(/x/m, a s ) and 6j(/x/m, a s ), 
that encode the high energy behaviour. The full relativistic Poincare symmetry [or 0(4) 
plus translations for the Euclidean space-time conventions adopted here] is not evident 
from the NRQCD Lagrangian of Eqs. ( |6.6|) - ( |6.14| ) in which only Galilean invariance is 



explicitly manifest. By imposing the full four- dimensional invariance to the respective 
order of the non-relativistic expansion, matching coefficients that accompany operators 
of different dimensions become related. This re-parametrisation invariance [[38 2| , |383|| is 
discussed in Refs. | |384| , |381|| . One result is the relation, 

cs(/J,/m, a s ) = 2c F (/i/m, a s ) - 1, (6.15) 

which can be derived by imposing invariance under an infinitesimal Lorentz boost, 
v — > v + 5v. Furthermore, the relativistic dispersion relation implies that the matching 
coefficients accompanying the kinetic terms are unity. 



80 



Those coefficients that are not determined by fundamental symmetries can be ob- 
tained by matching amplitudes, calculated in the effective theory, to their QCD coun- 
terparts. This can be done non-perturbatively on the lattice or to a given order in a s 
in perturbation theory. In the first case one would calculate a set of quantities that are 
sensitive to the choice of the matching coefficients in lattice NRQCD at finite values of 
the lattice spacings, a a < ir/q, for a range of quark masses a~ Xr K > m > a^ 1 and, ide- 
ally, match them to their continuum QCD counterparts. The difficulty of obtaining the 
continuum QCD result, which requires simulations at small lattice spacings, a T <C m _1 , 
has turned out to be prohibitive so far. The exception is the mass renormalisation, 
Sm, that can be fixed by demanding that the rest mass of an T meson must equal its 
kinetic mass. This is done by comparing finite momentum T masses with the expected 
dispersion relation^ |JU], 150 , |386|| . The only other attempt into this direction was an 
estimation of the coefficients of the order a s correction terms to Cp = 1 and 5m = 
from small volume simulations 13871, 



As an alternative to matching to continuum QCD, one could in principle treat all 
coefficients as free parameters and fix them by demanding T splittings, determined in 
lattice NRQCD, to match experimental input values. This, of course, would severely 
limit the predictive power of (NR)QCD calculations^]. Moreover, it is hard to combine 
experimental input and lattice NRQCD in a conceptually clean way; experiment has 
2 + 1 + 1 flavours of sea quarks of the right physical masses built in while lattice QCD 
calculations in general require extrapolations to the physical sea quark masses and, 
eventually, the relevant number of sea quark flavours. 

The general procedure of matching an effective theory to QCD is to start with the 
lowest dimensional operators, the dimension of which we assume to be n in terms of 
an expansion parameter, A, and to determine their Wilson coefficients in one or an- 
other scheme. Since the theory only has to reproduce QCD to the given order, n, of 
the expansion, the coefficients are ambiguous: corrections of order A can always be 
added. In the next step one would examine the set of operators of the next avail- 
able dimension, n + 1. These terms will not only undergo mixing with each other 
under renormalisation group transformations but also with lower dimensional opera- 
tors. The resulting set of coefficients to this order will depend on the conventions used 
to determine the lower dimensional ones: order A terms added to the coefficients that 
accompany dimension n operators have to be cancelled by operators of dimension n+ 1. 
This freedom of re-shuffling power corrections between ultra-violet Wilson coefficients 
and infra-red operators of higher dimension is nothing but the well known renormalon 
ambiguity ||26 1| , [262| , [263] , 264] . In conclusion, any matching scheme can be used but it 



2N 



Sm has also been computed to one loop perturbation theory in one version of lattice NRQCD [385]. 
On the lattice, where 0(3) rotational symmetry is broken, the (non-trivial) coefficients accompanying 
the kinetic terms can in principle be fixed by imposing the continuum dispersion relation. 

29 For a prediction of B meson properties at the 1 % level, it appears to be sufficient to consider 
the order mT 1 HQET/NRQCD Lagrangian 

mm. To this 

order, the only parameter that requires 
continuum QCD input is cf, such that the reduction in predictive power by using T fine structure 
splittings as an input is not great. 
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has to be employed consistently. 




Figure 6.2: Continuum inspired estimates of the matching coefficient, cp(fi), as a func- 
tion of the lattice spacing for bottom quarks. The upper curves correspond to /i = ir/a, 
the lower ones to /i = 1/a. 



In Refs. |[381| , [380|1 it has been argued that despite of the fact that the leading order 



NRQCD Lagrangian differs from its HQET counterpart by the kinetic term it incor- 
porates, HQET delivers a viable prescription for determining the NRQCD matching 
coefficients, at least up to order m~ 3 . In HQET the Fermi coefficient, Cp, is known to 
two loops and all other coefficients at the one loop level in the MS scheme of dimen- 
sional regularisation. We display the results in Appendix [E]. As discussed above, the 
coefficients are specific to the prescription used. Other regularisation schemes or ways 
of organising the expansion, e.g. in powers of c _1 , will in general yield different results. 

Unfortunately, no lattice NRQCD perturbation theory results exist for the matching 
coefficients, Cp and cp>- Large contributions from lattice tadpole diagrams | |12U| | in 
general result in big renormalisations between non-spectral quantities, calculated by use 
of lattice regularisation, with respect to continuum schemes such as the MS scheme. 
This is for instance reflected in the ratio [ 391 , |392j , |393| , |394 
the pure gauge Wilson action. 



A^/Ai « 28.81A L , for 



It has been suggested ||1 19| , |120|| to (partially) cancel tadpole contributions by divid- 
ing each lattice link that appears in a given operator by the fourth root of the measured 
expectation value of a plaquette, Up = Re (Tr U X;/XU ) . Other prescriptions using the 
average link in Landau gauge or expressions containing the logarithm of the plaquette 
have been suggested as alternatives ||120|| . It is argued that the diagrams that result in 
large renormalisations also cause the plaquette (or the average gauge fixed link) to sub- 
stantially deviate from the free field expectation of unity. Moreover, such (ultra-violet) 
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a'VGeV 



Figure 6.3: The same as Figure |6.2| for charm quarks 



renormalisation effects might commute with the infra-red physics of interest. Based on 
these ideas, the following replacement has been suggested ||119|| , (3S = PmpS^, with 
Sir — S/Up. This yields the relation, c*mf = cupjUp. Another popular choice of an "im- 
proved" coupling is ||120| , ofnal = — liaUp/(iiCp). From the perturbative expansion 
of the plaquette expectation value [|395| , |39(j|| , one finds ps 2.63A M f ~ 4.19A F nal 
for rif = 0; indeed, the coefficients appearing in the one loop perturbative matching 
between ajpj(ir/a) and q;mf(o) or ofnaiX^) are much smaller than for the bare lattice 
coupling. 

After obtaining a "tadpole improved" version of lattice NRQCD by following the 
above recipe, one might hope that the running of the coefficients with the quark mass 
closely resembles that of continuum NRQCD. Unlike dimensional regularisation, the 
lattice imposes a hard cut-off, ir/a, on gluon momenta, q, such that it is not entirely 
clear what matching scale corresponds to the /i of the MS formulae of Appendix [Ej. 
It is reasonable, however, to assume ||120 |, n/a > [J, > 1/a. In Figure |6.2| , we display 
the resulting estimates for cp for the bottom quark as a function of the inverse lattice 
spacing, based on Eq. ( |E.1| ). The widths of the two bands correspond to the above scale 
uncertainty, oi-^g as a function of the scale has been obtained by running down the 
quenched result ||103|| of Eq. ( f4.76|) from a high energy scale by means of the four loop 
P function, Eq. ( $72^ . 

Within the region, 1.5 GeV < a -1 < 3 GeV, cp can easily deviate from the tree 
level value by as much as 15 %. Moreover, the two loop result significantly deviates 



from the one loop prediction, indicating a slow convergence of the perturbative series. 
In Figure |6.3| , it is convincingly demonstrated that for lattice resolutions better than 
1 GeV a perturbative estimation of cp for charmonia is unreliable. On the other hand 
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Figure 6.4: Continuum inspired estimates of cd(/z), as a function of the lattice spacing 
for charm and bottom quarks. The lower curves correspond to fi = ir/a, the upper ones 
to ji = 1/a. 



lattice spacings, a" 1 <C 2 GeV, are too big to sample the relevant bound state dynamics 
and would result in huge scaling violations. Finally, in Figure |6.4j we plot our estimates, 
Eq. ( E.2|) , of the cut-off dependence of the Darwin coefficients, cp, for both, bottom 
and charm quarks. We find cd to vary much more with the quark mass than cp. 
Unfortunately, no two loop calculation for this quantity is available. 

In conclusion, the perturbative calculation exhibits a significant dependence of the 
coefficients on the quark mass (or lattice spacing). While at a given lattice spacing cp 
decreases with increasing quark mass, cjj shows the opposite behaviour. Perturbation 
theory seems to be slowly convergent. Moreover, in general power corrections can 
contribute to the coefficients. 



6.2 Lattice NRQCD 

It is straight forward to discretise the Lagrangian, Eqs. (|6\6|) - (|6.8| ), and to simulate it 
directly on a lattice. Let us start with the leading order continuum NRQCD Lagrangian, 



£nrqcd,v 2 — W 



D D ' ^ 

4 2ra^ / 



(6.16) 



Note that from now on we use, c = 1. We define the heavy quark propagator, K = 
(ipip '). K is the direct product of a 2 x 2 matrix acting on the Pauli spinor space, a 3 x 3 
matrix acting on colour space and L\L T matrix acting on space-time. From 
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Eq. ( |6.16 ) it follows that the evolution of K with time is governed by the Hamiltonian, 



D 2 

Ho = -r— : (6.17) 
-d 4 K = (igA 4 + H )K. (6.18) 

By formally solving the above differential equation, we obtain the evolution equation, 

rs(t)=x ( r rt+a 1 •) 

K(x,t + a) = Y DsVlexp - dt' {igA 4 + H ) \K(y,t), (6.19) 

V Js(0)=y I L Jt J J 



y 



where we have assumed the sum over all paths to be appropriately normalised. The 
initial condition reads, 

K(x)\ X4=0 = 5 3 ( X ). (6.20) 

Note that we have suppressed the dependence of the propagator on the source point, 
K(x) = K(x,y = 0). 

A natural discretisation of Eq. (|6.19 ) is |35] , 



K(t + „) = (i _ a -^±^Jut i (l - 2^) \(t). (6.21) 

We have omitted the dependence on the spatial coordinates from the above equation. 
The temporal link, ll} 4 , is diagonal in space. The covariant Laplacian within H (t) can 
be written as, 

Dx y W = a- 2 £ + Ul^A-^y ~ 2S *y] . ( 6 - 22 ) 



up to 0(a 2 ) lattice artefacts. For the naive n — 1 discretisation, the evolution equation 
might become numerically unstable as 1 — aH Q becomes negative for momenta larger 
than the quark mass; lighter quarks try to travel faster than they are allowed by the 
evolution equation. Introducing the stabilisation parameter, n, improves the spatial 
propagation and relaxes this criterion to m&x(aHo) < n. In the free field case, the 
maximal eigenvalue of the Laplacian, Eq. ( |6.22| ), is J2i max PiPi = 3a -2 , such that 
ma > 3/(2n) has to be maintained. When switching on interactions, the factor 3/2 is 
reduced somewhat. 

Working with an anisotropy, a T < a a , offers an alternative to introducing the param- 
eter, n. In this case, the free field stability criterion relaxes to, ma a > (3/2)a T /a a . It is 
amusing to see that in lattice NRQCD simulations, discretisation effects become more 
pronounced in light quark propagators, rather than for heavy quarks as in relativistic 
lattice QCD. While in the latter case, heavier quarks can be realised by reducing a T , in 
NRQCD lighter quarks require smaller a T (or larger n). Of course one would not rely 
on results obtained for quark masses, m < a' 1 , as the non- relativistic expansion breaks 
down for a cut-off on gluon momenta larger than the quark mass. On the other hand 
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one would also not want to simulate quarks much heavier than the lattice resolution 
to keep the scale, mv, separated from the lattice cut-off, m <C (a CT t>) _1 . Otherwise, 
the matching coefficients between lattice NRQCD and QCD would explode and their 
behaviour could no longer reliably be estimated. 

The above evolution equation approximates the continuum equation only up to 
0(a T ) lattice artefacts. These can be removed at tree level by the substitution, H — > 
H [1 + (a T /4n)Ho] (or reduced by increasing n). 0(v 4 ) correction terms, 5H V 4, can be 
included too, 



K(t + a) = 1-^1 1-^ (6.23) 



x 




Details can be found in Ref. ||35| . Although a method to incorporate the four fermion 
terms of Eq. ( |6.8| ) is suggested in this Reference too, these have not been included into 
any lattice simulation so far. Typically, tadpole improvement is employed in NRQCD 
simulations, i.e. link variables are divided by factors Up 4 ' or equivalent quantities that 
approach unity in the continuum limit. 

The NRQCD evolution equation has also been applied to the heavy quark within 



heavy-light systems [ |389|| . H does not only consist of the static propagator but also 
incorporates the kinetic term, while the Fermi term, that is of the same order in m , 
appears within 5H. The main advantages of this procedure over a naive discretisation 
of HQET lie in smaller wave function renormalisations and in a reduction of statistical 
fluctuations. Both effects are related to the use of a propagator that samples gauge 
fields over an extended spatial region. The disadvantages in applying lattice "NRQCD" 
with HQET like power counting to heavy-light mesons is a loss in conceptual clarity as 
the wave function renormalisation depends on the expansion parameter, m~ l , in a way 
that cannot be absorbed into multiplicative field redefinitions. 

By contracting quark and anti-quark propagators with suitable combinations of 
gauge transporters and Pauli matrices, particular s lj states can be realised whose ground 
state masses can be extracted from the asymptotic decay of two-particle Green func- 
tions in Euclidean time in the usual way. Like in all direct spectrum evaluations, radial 
excitations present a major problem. Thus, it is a tremendous achievement that the 
3S as well as the 2P states have been determined, with statistical errors of about 



100 MeV ISSBl, IToTI IT33, H531, |3m Precision results exist for 2S, IP and IS states. 



6.3 The potential approach 

6.3.1 Deriving a bound state Hamiltonian 

We wish to derive a Hamiltonian that governs the evolution of a quarkonium state from 
the order c -2 (or v A ) NRQCD Lagrangian of Eqs. (6.G) - (6.14) that is formulated on 



the quark-gluon level. As a first step in this direction, we calculate a heavy quark 



St. 



propagator in a representation that will turn out to be suitable for our purpose. The 
time evolution of the Pauli propagator, K, is controlled by the equation, 



dJC = HiK, 



(6.24) 



where the Hamiltonian, 

H 1 =m + igA 4 + (6.25) 

can be read off from Eq. (|6.7|) . Unlike in Eq. ( |6.18| ) we decide not to eliminate the heavy 
quark rest mass, m = m^ — 5m, and not to rescale the (anti-)fermion fields by factors, 
exp(±mt). For the initial condition, 

K(x,y)\ X4=m = 6 3 (x-y), (6.26) 

Eq. (|6.24j ) can be formally solved by summing over all possible paths connecting y with 

z(x4)=x 



X 



K(x,y) 



dt [pz - ^!(z,p)] 



(6.27) 



Dz Dp exp 

Mv4)=y 

where the dot denotes a derivative with respect to the time coordinate. The correct 
normalisation is assumed to be included into the definitions of Dz and Dp. 

We can now combine two such propagators into a generalised (fluctuating) rectan- 
gular Wilson loop, 



Giffij(r,T) = (Tr [UoKtiiuuX! )U T K rj (x 2 ,y 2 )} ) , 



(6.28) 



where the indices i,j and represent the spins of the initial and final states. Note 
that G, unlike the argument of the expectation value, is real in Euclidean space-time. 
For the case of quark and anti-quark having different masses, two different propagators, 
and K x , must be used within the above formula. We denote the temporal extent 
by, t = 2/1,4 — x\a = 2/2,4 — £2,4, and the spatial separation by r = y 2 — yi = x 2 — xi. 
The situation is visualised in Figure |6.5| . For the sake of simplicity, we switch to leading 
order NRQCD with equal quark masses. In this case, 



H\ = m + igA^ + 



P 



2m,,, 



(6.29) 



To lowest order in v/c, the exponent within Eq. (|6.27 ) can be approximated by the 
value it takes along the shortest path |26|, Thus, 



G(r, r) = exp 



/ dty 



m 



P 



2m„ 



(W(r,r)). 



(6.30) 



In higher orders of the v/c expansion fluctuations of the propagators around the classical 
paths have to be taken into account that result in additional terms. From the spectral 
decomposition of the Wilson loop, 



(W(r,r)) oc exp[-Fo(r)r] 



oo), 



(6.31) 
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Figure 6.5: The four point function of Eq. ( |6.28D 



we arrive at, 



with 



It 



G = HG, 



p- 



H = 2m + t- +V (r), 



(6.32) 



(6.33) 



in the limit of large r: the result is a Schrodinger equation that governs the evolution 
of a quark anti-quark state in a gluonic background whose average effect is contained 
in the static potential, Vq. The validity of the instantaneous approximation is tied to 
that of the naive quark model: if quarkonium states can be completely classified by the 
quantum numbers of the constituent quarks, then the spectrum and wave functions can 
be obtained by solving the quantum mechanical equation, 



Hi[) n u 3 (r) = E n iip n n 3 (r). 



(6.34) 



We have discussed above that m will in general differ from the "kinetic" mass of the 
quark, = m(/i) + 5m(fi). Furthermore, this difference will depend on the matching 



88 



scale, fi. In Section f|]5] we have also seen that the potential, Vq{t, fi) = Vo(r) + K e if(A*), 
can be factorised into a physical and a self-energy part. This observation results in the 
relation, 

Sm(ji') = 6m(ji) + i [V S M - V scU (fi)] , (6.35) 

i.e. Sm(fi) diverges as fi — > oo. 

In lattice NRQCD, it is straight forward to calculate masses of quarkonia states, 
Ey(p), projected onto non- vanishing momentum, p. We use the convention that Ey is 
the sum of the bare quark masses, 2m, and the energy shift due to the interaction terms 
of the Hamiltonian. By requiring the correct dispersion relation to the given order of 
the expansion, 

StW-B t (0) = ^--^ ( +... ), (6.36) 

the T rest mass, my, can be determined. The mass shift, then, is given by, 28m = 
Ey(0) — my. Within the potential approach, the zero point energy at first appears to 
be difficult to determine in a similar way. However, in principle it should be possible 
to calculate potentials governing finite momentum quarkonia states too, by starting the 
derivation from a boosted NRQCD Lagrangian. 



6.3.2 Relativistic corrections 



The Hamiltonian Eqs. (|6.31|) - ( 6.33 ) was first obtained in Ref. [ 175 1 in a systematic 



way from continuum QCD where the static Dirac equation, (74-D4 + m) q = 0, is solved 
by a Schwinger line times a factor, e~ mt , after projecting onto quark and anti-quark 
states. Starting from QCD, Eichten, Feinberg and Gromes derived spin dependent 
correction terms [f^, |25j (see also the article by Peskin [§4}]). Finally, Brambilla 
and collaborators (BBP) |26], B7I] found an additional relativistic correction term to 
the central potential that had previously been ignored and added further velocity (or 
momentum) dependent terms by taking fluctuations of the heavy quark propagators into 
account. In general, apart from the one-particle Lagrangians, Ci,i = ip,Xi [Eqs. (|6.6|) 
- fl6.7p ] the two-particle Greens function receives contact term contributions from the 
two-particle sector Lagrangian of Eq. ( |6.8| ) f28|j , C^ x . Taking these into account too, the 
complete result to this order in c , with the NRQCD matching coefficients included [E8, 



p9| . 398 1, in the centre of mass frame (p = Pi = — P2 and L = = L 2 ), for m x > m.2 
is, 

H = 21 m; - bmi + 



+ V(r,p,L,S 1 ,S 2 ), (6.37) 

where the potential, 

V(r, p, L, S x , S 2 ) = V (r) +V C + V SD (r, L, S x , S 2 ) + V MD (r, p), (6.38) 
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contains corrections to the central potential (C) as well as spin dependent (SD) and 
momentum dependent (MD) corrections. V (r) denotes the static potential while, 



V c (r) = 



d. 



17111712 



-AnC F a s 8*{r) 



+ E 



,0 



+ 



1 1 

mi m 2 



Wo(r) + V 2 Vf (r)] + cf\ 2 V a B {r) 
V f (r), 



(6.39) 



Vsv(r, L, Si, S 2 ) 



+ 
+ 
+ 
+ 
+ 



2r 



S 1 + S 2 l c + K 2 / (^) 
m\irt2 r 

\ml m\) r 
S 1 _-S 2 _ L c_^(r) 

mim 2 r 

Si^l (i) (2) D T/ / x 

mim 2 

S l • S 2 ^(1)^(2) 

3mim 2 



(6.40) 



and 



with 



Vmd(f,p) = 



17111712 



{Pi,Pj, [$ijVb( r ) - R ij V c(r)}} Weyl 



1 



Rij 

Srrii 

Si) 

C F,D 



C± = 

d s = 
d„ = 



Wcyl 



r 2 3 ' 
5m(mi,fi), 

CF,D{mi,H), 

c±(m 1 , m 2 , n) = - (c ( P ± cg } ) , 

— — I [d ss (m 1 ,m 2 , £t) + C F d vs (m 1 ,m 2 , //)] , 

AirC F a s 

— — I [4„(mi, m 2 , //) + C F d vv (m 1 , m 2 , . 

AnC F a s 



(6.41) 



(6.42) 

(6.43) 
(6.44) 

(6.45) 
(6.46) 
(6.47) 
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The symbol {a, b, c}w ey i = \{ a , {b,c}} denotes Weyl ordering of the three arguments. 
Note that in the equal mass case, that has been considered in Ref. [f23|, where c_ assumes 
its tree level value, c_ = 0, two of the spin-orbit terms vanish. The term proportional to 
V, in Eq. ( |6.39|) has been identified very recently |}2 and in principle additional 1/m 2 



corrections to Vq should exist [ 32(1 , albeit to higher order in the c _1 power counting than 
order c~ 2 considered above. 

The last term of Eq. (|6.40|) has been written in a somewhat suggestive way that is 
motivated by the expectation, V^r) ~ 8irCFa s ^ 3 ( r )- c s has been eliminated from the 
above formulae by using the re-parametrisation invariance relation, Eq. ( |3.15| ). Note 
that neither d c vv or d c vs nor d ss or d sv contribute to d s or d v . This means that even in 
the equal mass case, where d c vv = —na s + • ■ ■, d s and d v are of order a s . The one loop 
results in the MS scheme are displayed in Eqs. (|E.22| ) and ( |E.23| ) of Appendix [E]. 

Vo, V h V 2 Vf, V 2 V^, V{,...,V4 and V&, . . . , V e can be computed from lattice cor- 
relation functions (in Euclidean time) of Wilson loop like operators. The functions 
V{, . . . , Vi are related to spin-orbit and spin-spin interactions. The MD potential gives 
rise to correction terms of the form -L 2 , 4s L 2 , -p 2 , and S 3 (r), and the correction to 
the central potential includes the expected Darwin term, V 2 Vb, as well as V 2 V^ and 
V 2 K r 



a 



6.3.3 Scale dependence 

The SD potentials, V{, . . . , V4, as well as V 2 V^ and V 2 V^ depend on the matching 
scale, /i. The potentials Vo, V, as well as V&, . . . , V e can contain additive, [i dependent 
self energy contributions; however, their derivatives are scale independent^]. 

Due to Lorentz invariance, certain pairs of potentials are related to the static po- 
tential by the Gromes ||399|1 and BBP [27fl relations, 



V^;r)-V((fx;r) = Vtf(r), (6.48) 
V b (r;ri + 2V d {r;») = ~F '(r) - ±V {r; »), (6.49) 
V c (r) + 2V e (r) = -^(r), (6.50) 

such that three potentials, e.g. V[, Vd and V e can be eliminated from the Hamilto- 
nian. Note that Eq. (|6.48|) implies Eq. (|6.15|) . Given the structure of the Hamiltonian, 
Eqs. (|6371 ) - O), and the Gromes relation, Eq. (gTg ), we can deduce |DD|, g9j the fol- 
lowing relations between potentials, evaluated at cut-off scales fi and \i! ', by demanding]^] 

30 Vq is a spectral quantity while V, and the MD potentials V&, . . . , V e originate from the terms D4 
and D 2 /(2m) of the NRQCD action that are protected by reparametrisation invariance. Therefore, 
these potentials do not undergo multiplicative renormalisation. 

31 The S function within Vc represents a problem: no spin- and momentum-independent counter term 
is known that has the right mass dependence to cancel the running of the coefficient, d s . This might 
hint at further, not yet discovered, relations. 
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dH/dhin = 0, 

W a V;r) = -^{c D (riVX E (^r) + [c D ( f i)-c D (»')]V 2 V (r) 





+ c?M V 2 V*{n- r) + cW) W a B (//; r)j , 


(6.51) 






1 c F (//) 

C F (//) 




(6.52) 


Vfo';r) 






(6.53) 








(6.54) 


V^';r) 


= cU^)i cUfi)v ^ r) 








- 12[d v (p) - 


^(/i')]47rC F a s 5 3 (r)}. 


(6.55) 



Since the potentials, appearing in the above relations, do not depend on the quark 
mass, the ratios CF 5 £)(m, fi)/cF,D(m, must not depend on m. Therefore, the matching 
coefficients can always be factorised into two separate functions, Cj(m, /i) = /j(m)^ l ~ 1 (^). 



6.3.4 Integrating out gluons 

We have managed to separate the time dependence of the interaction into coefficient 
functions of various interaction terms, Vi, which we shall call the potentials. These 
potentials can be computed as expectation values in presence of a gauge field back- 

ground mm mm, 

= 2/fim fcft((E(0,0)-E(0,t)))^, (6.56) 

J 

= 2g 2 lirn Tci* ((B(0, 0) -B(0,t))) w , (6.57) 

J 

2 

V,(r) = - 9 - lim fdtt((E(0,0)-E(0,t)» c w , (6.58) 

where the superscript V denotes the connected part 32 , 

{{AB)) C W = ((AB)) W - ((A)) W ((B)) W . (6.59) 



rB I 



For the SD potentials one finds 

!± 

r 



V({r) = e ijk g 2 \im[ T dtt((B i (0,0)E j (0,t))) w , (6.60) 

J 



32 Both, electric and magnetic fields transform oddly under charge conjugation. Therefore, in SU(2) 
gauge theory, where all traces are real, ((E)) w — ((&))w = 0- Under PC transformations the electric 
field transforms evenly. However, the magnetic field has PC = —1. Therefore, ((B)) w — still holds 
for SU(3) gauge theory. However, components of ((E)) w that are not orthogonal to r do not have to 
vanish (cf. Table [□]). 



92 



V'(v) 



X -e i3k g 2 Um £dt t((Bi(0, 0)E j (r, £))) 
2g 2 lim f T dt [((BiiO^BjfaWw 



8 a 



«B(O,0)-B(r,t)» 



y 4 (R) = 2g 2 lim [ T dt((B(O,0)-B(r,t))) w , 

J 

where Rij is defined in Eq. Q6.42 ). Finally, the MD potentials are, 
H(r) = 



RijV c (r) 



VJr) 



RijV e (r) 



- l -g 2 lim J^dtt 2 ((E(O,0) -E(r,t)» 

/lim fd« 2 [((^(O,0)^(r,t))) 
o 



c 



c 



^((E(0,0)-E(r,t)))^] 



6* 



g 2 lim / ^i 2 ((E(0,0)-E(0,t)» c w , 



g 2 lim / dtt 2 [((Ei(0, 0)^(0, t)» 



(6.61) 
(6.62) 

(6.63) 

(6.64) 
(6.65) 

(6.66) 
(6.67) 



- ^((E(0,0)-E(0,t)))^]. 

While V , Vfe, . . . , V e have the dimension m, V,, V{ and V{ have dimension m 2 and V 3 , 
V4, V 2 Vq E and V 2 V^ have dimension m 3 . 

Throughout the previous equations, the expectation value, ((FiF 2 ))w, is defined as, 

_ (TrV[exp(igf BC dx ll A ll )F 1 F 2 ]) 



((FiF 2 )) 



W(C) 



(6.68) 



(Tr V [exp (ig J sc dx^A^)}) 

where 5C represents a closed path [the contour of a Wilson loop, W(r,T), T > r]. The 
nominators of Eq. ( |6.68|) that are required to compute the potentials are depicted in 
Figure |61| The correlators appearing within the coefficient functions of the spin-orbit 
potentials, V{ and V 2 \ involve electric and magnetic fields, the latter originating from 
the angular movement. Correlators between two magnetic fields are required in the 
spin-spin potentials, V3 and V4, which arise from interactions between the two Fermi 
terms oc g Si ■'B/rrii. The corrections to the central potential, V 2 V^ and V 2 ^ 5 , involve 
electric-electric and magnetic-magnetic interactions, respectively, while V, and all MD 
corrections involve two electric field insertions. The latter arise from re-expressing 
derivatives acting on the static propagators in terms of field strength insertions. 

In principle similar results that would include Wilson loops with more than two 
field strength insertions can be obtained from the order c -4 (or v 6 ) NRQCD Lagrangian 
of Ref. P5 |. This tedious work has not been done yet since the dominant sources of 
error at present are the uncertainties of the matching coefficients and certain transition 
matrix elements (cf. Sections 772 , 6.5.2 and 6.5.4 ), rather than higher order relativistic 
corrections. 
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Figure 6.6: The nominators of the arguments of the integrals within Eqs. (|6.56|) - ( |6.67|) . 
6.3.5 The potentials as perturbations 



In Refs. ||401| , p^l spectral decompositions of the above potentials have been derived. 
The results can be written as follows, 



V 3A (r 
V*V a E > B (r 

V,{r 
V b , c , d ,e{r 



^ n [ ) Jo ^ A7.fr) 



n>0 



n>0 



n>0 7 ° 



D {E,B) [ 



r 



n>o A7(r) 
V £> (1,2) M f T dtte- AVn{r)t = V jP "' 2) ^- > 

T D' (r) f T dtte~ AVnir)t = V 



n>0 



/•t +2 r\(b,c,d,e)( r \ 

V- n (b,c,d,e)( r ) / Af l _ p -^V n {r)t _ \p Ijj 

^ " ( } A, 2 " £ [A7(r)P : 



n>0 



(6.69) 
(6.70) 
(6.71) 
(6.72) 
(6.73) 



where AV n (r) = V n (r) — Vq(t) denotes the difference between the nth hybrid excitation 
and the ground state potential. The coefficients, D n (r), are real parts of products 
of two transition amplitudes and can easily be read off from Eqs. fl6.56| ) - ( |6.67| ). For 
instance, in the case of V4, one obtains, 



D*(r) = 2g 2 Re [($ ri0 |B(0)|$ r , n )($ r , n |B(r)|$ r , )] 



(6.74) 
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where |$ r ,n) denotes the nth excitation of a quark anti-quark state at separation, r, and 
the states are thought to be normalised, (^ r ,i\^r,i) = 1- Note that, = —4:D' n = 6D^. 

Physically, the above result can be interpreted as follows | |402| , [401 ]: at time the 
spin of the quark at position, 0, interacts with the background glue and excites the 
gluonic vortex until, at time r a second interaction with the spin of the anti-quark at 
r takes place that returns the flux tube into its ground state: the non-perturbative 



analogue of a gluon exchange! From Table |5.1| one can read off that in general the 
intermediate state will be a superposition of excitations within the S~ and U u channels 
in the particular cases of V 3 , V4 and V 2 V^ . 

We add the term proportional to V 2 V^ of Eq. ( |6.39[ ) to the two terms proportional 
to V3 and V4 of Eq. (|6.40|) . The result reads, 



Vjr) 



(6.75) 



+ 




W ». 



By inserting the spectral decomposition with the correct coefficients, D n , determined 
from Eqs. (|6.57| ), (|6.62| ) and ( |6.63|) , into Eq. ( |6.75| ) one obtains, 



VJr) 



V l - 

to VJ?) - V (r) 



(6.76) 



x 



r.O 



' (1) pSx-B(0) , (2) gS r B(r) 

Cp "T Cp 



mi 



m 2 



where we have exploited the fact that B 2 = 4(S ■ B) 2 . The result is exactly the energy 
shift one would have expected in second order perturbation theory from the Fermi terms 
for quark and anti-quark within Eq. ( |6.7| ), 



with 
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B(x) 



<5 3 (x) + 



cyga 2 ■ B(x) 



5 3 



x - r 



(6.77) 



(6.78) 



2mi 2m2 
where Oi = 2Sj. 

Other potentials in their spectral representation can be interpreted as perturba- 
tions too. However, relating these to the NRQCD Lagrangian requires somewhat 
more involved formal manipulations. From the considerations above it is obvious 
that the formalism cannot readily be applied to spin dependent interactions of hy- 
brid quarkonia where the ground state would appear as an intermediate state: since 
exp[— (V n — V m )t] diverges with t for m > n; the matrix elements corresponding to 
Eq. ( |6.77| ), for an external state, |3> r ,m)j cannot be obtained from a simple time integral 
over double bracket expectation values. 
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6.4 Model expectations 



We discuss expectations for the potentials and the resulting Hamiltonian, and discuss 
the Lorentz structure of the effective interaction kernel. 



6.4.1 The potentials 



We will present simple model expectations for the above potentials. The double bracket 
expectation values of colour field operators can be obtained from infinitesimal deforma- 
tions of a generalised, non-static Wilson loop 33 [ |358| , 



i'U\X)Fpa{y)))w 



5 2 \n(W) 



6S flu (x)5S pa (y)' 



(6.79) 



If the functional dependence of the Wilson loop expectation value on its contour is 
known, the above formula can be applied to calculate the corresponding long distance 
behaviour of the potentials. This has been done for the stochastic vacuum model (SVM), 
dual QCD and the area law assumption in Refs. IHO, 6l|. A variety of predictions on 



SD and MD potentials exists in the literature that are based on effective modified one 
gluon exchanges or Bethe Salpeter kernels. Ref. ||403|| represents a recent example^. 
Given these different suggestions, lattice results with a precision that is sufficient to 
discriminate between them are highly desirable. 

combined with tree 



Here, we shall only discuss the area law expectations |26|, [27 



level perturbation theory and constraints from the renormalisation group mixing be- 
tween the potentials [ffil , Eqs. (|6.51|) - ( |6.55D , 



Vo{r;ii) 

V 2 V a E {r;n) 
V 2 Vf(r;/i) 
V{{r;n) 

V 3 (r;v) 

V 4 (r;fi) 
V,(r) 

V b (r; fi) 



c E a (a*) - 

Cf(A*), 
h{fi) 



■ - + ar, 
r 

2a + &(//) 



a. 



h(fj) 



r' 



8n[e-h(/j,)}5 3 (r), 
0, 

^ i \ 2 e a 

CbW + »-- a r > 
3r 9 



33 For the definition of the functional derivative acting on a Wilson loop with respect 
element, S^(x), see e.g. Ref. [19S|. 



(6.80) 

(6.81) 
(6.82) 
(6.83) 

(6.84) 

(6.85) 

(6.86) 
(6.87) 

(6.88) 

to a surface 



:S4 



Note, however, that their result is incompatible with Eqs. (3.64) - (S.67). 
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V e (r) = ~--^r, (6.89) 
2 r o 

V d (r;/i) = C d (jj,)-^r, (6.90) 

V e (r) = -|r, (6.91) 

withQ e — /i(/i) ~ Cfoi s . The above formulae conform with the Gromes and BBP 
relations, Eqs. (|6.48|) - ( |6.50|) . While Vq, V, and the MD potentials do not undergo mul- 



tiplicative renormalisation, Vq, V& and V d still contain additive self-energy contributions 
that will diverge as /x — > oo and whose /x dependence has to be cancelled by the quark 
mass shifts^, Srrii. The constants, V^iKaO, C,f (/x), C,f (/x), Cb(/x) and Cd(/x) as well as all 
terms proportional to e originate from perturbation theory, while all terms proportional 
to the string tension, a, are due to the area law ansatz, with the exception of V 2 V^ . 
We allow for terms proportional to h in V[ and V% that are thought to originate from 
the mixing between these two potentials under renormalisation group transformations, 
Eq. (|6.52p . In perturbation theory as well as in vector exchange models, one obtains, 
Vz(r) = V^/r — V2 and V4 = 2V 2 V-2. Therefore, replacing Cfol s by e — h within these 
potentials appears to be reasonable. However, we remark that further corrections to 
this ansatz must exist since the scaling behaviours under /x — *• /x' of V2 [Eq. ( |6.53| )1, V3 



[Eq. ( |6.54D 1 and V4 [Eq. (|6.55|) 1 are incompatible with each other. Finally, the expec- 



tation, 2<t + 6(/x), within V 2 Vf is motivated by the lattice results to be presented in 
Section |6]6| as well as by dual QCD and SVM calculations |K| . One would expect 
additional 5-like contributions to V 2 V^ and V 2 V^ from Eq. ( 6.51]) , which we ignore for 
the moment. 



Interestingly, by adding a perimeter term to the Wilson loop area law [^, |2j 
one obtains a non- vanishing Cd = —V se \f/4, which agrees with the expectation from 
perturbation theory. However, the perimeter term does not contribute to C^, or 
Cb- In continuum perturbation theory as well as in lattice perturbation theory in the 
infinite volume limit, one obtains the tree level results |2l 



C b (fM) = 0, C d (ji) = ~V ae M, (6.92) 

where the lattice perturbation theory result for V se if(a) with the Wilson action is given 
in Eq. (|4.53|) . By using the lattice field definitions of Eqs. ( |5.3| ) - (|5.6|), we obtain the 
lattice perturbation theory results [ |j| , 

Cf (a) = -C F a s a- 3 x 7.91084 . . . , (6.93) 
Cf(a) = C F a s a~ 3 x 14.89413 ... , (6.94) 

for the other two self-energy contributions. 

35 In one loop perturbation theory one obtains Q, V,(r) = —CFCAOL 2 s /(^r 2 ). 

36 We ignore the possibility of self energy contributions to V,, V c and V e that vanish at least to lowest 
order perturbation theory in the parametrisation. In lattice determinations of V c and V e these have 
indeed been found to agree with zero within errors 0] . 
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6.4.2 The Hamiltonian 



The Hamiltonian that results from the ansatz Eqs. (|6.80|) - (16.91 ), in the equal mass 
case, m = mi = m 2 , takes the form, 



with 



H 
H 
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Rij S 1 5*2 



1) -1(1 + 1) -8(8 + 1)], 

6(L ■ S) 2 + 3L ■ S - 2s(s + 1)1(1 + 1) 
" 6(2/ - l)(2/ + 3) ' 



(6.95) 



(6.96) 



(6.97) 
(6.98) 

(6.99) 



(6.100) 

(6.101) 
(6.102) 

(6.103) 

(6.104) 

For a discussion of the non-equal mass case we refer the reader to Ref. |29]. The 
parametrisations of the potentials that enter the above Hamiltonian, can of course be 
improved in several ways, for example ||398| , [404|| by including the known one loop 
perturbative results for the spin dependent terms ||405| , |367|| and the two loop result for 
the static potential ||260| , |258|| . Note that all terms containing the low-energy parameter, 
a, are independent of the matching scale, \i 3> y/a. 

We have eliminated Cd from the above Hamiltonian by use of Eq. Q6.49 ), V^eif = 
— 2Cb — 4:Cd- The subscripts of the correction terms, 5 Hi, do not necessarily relate to the 
potentials of origin. H contains contributions from Vo, Vc as well as from Vmd while 5Hs 
contains terms due to Vc and Vmd- We have used the relation 27r(5 3 (r)) = — i(r~ 3 r • p) 
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to cast a term that appears within Vmd into a 5 function. The radial Schrodinger 
equation, Eq. (|A.3| ), can be solved numerically for H and, subsequently, the 6 Hi terms 
can conveniently be treated as perturbations. 
We substitute, 



. Keif + 4C 6 

m = m H (6.105) 



into H . To order m 2 this yields, 



2 

H = C + 2rh + ^ + V(r), (6.106) 
m 

V(r) = -- + ar, (6.107) 
r 

C = _25m-4C 6 + i ^( Cl? Cf + c|Cf), (6.108) 
3cd& + 2cr 



12m 2 



(6.109) 



the static quark self-energy shift, Keif, is eliminated from the Hamiltonian. This was first 
noticed in Ref. ||27| . The remaining scale dependence of 5m(n) has to compensate that 
of the sum of the (small) term Cb(n), which vanishes in tree level perturbation theory, 
and the term containing and C^, which is suppressed by a factor mT 2 . Moreover, 

and Ca have different relative signs, such that partial cancellations occur. Of course, 
the above substitution is only valid for quark masses, m 3> V^if oc a s fi. This relation, 
however, is automatically fulfilled for matching scales, // < m. In conclusion: the mass 
shift, 5m, which is related to the wave function renormalisation, becomes reduced as 
relativistic corrections are taken into account. 

Some of the correction terms are well known from atomic physics, others are specific 
to non-Abelian gauge theories. One piece of 5H$ [Eq. Q6.98Q 1 as well as the term pro- 
portional to b/r within Eq. ( |S.96| ) stem from the Darwin interaction. A string whose 



energy density, a, is carried by a constant longitudinal electric field |^2| , gives rise to 
the (classical) orbit-orbit interaction term, — <r/(6r)L 2 , that appears within Eq. ( |6.99|) . 
5HsT> [Eq. (|6.100|) 1 contains a spin-orbit (Thomas) interaction term that, unlike its QED 
counterpart, only falls off like r -1 at large distances. In addition, it contains two spin- 
spin interaction terms that take very much the same form as in QED, the first of which 
does not affect S waves and the second of which only affects S waves to the order in a s 
considered above. 



6.4.3 The Lorentz structure of the effective interaction 

The general form of a Hamiltonian governing relativistic two-particle bound states has 
been derived within the Bethe-Salpeter formalism (see e.g. | |227f and references therein), 
under the assumption that the interaction kernel only depends on the transfer mo- 
mentum, q 2 : The momentum space kernel can be decomposed into the five Lorentz 
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invariants, 



I = V s 1 <8> 1 + Vv 7m ® 7m + Vr-<V <g> cr^ + 7m75 <8> 7m7s + Vp 75 <8> 75, (6.110) 

where the form factors, V, only depend on q 2 . In the QED case, within the ladder 
approximation, only Vy assumes a non-trivial value and the resulting Hamiltonian has 



the Breit-Fermi form, well known from atomic physics. In the most general case [[227 
the equal mass Hamiltonian reads, 



V{r) = V v (r) + V s (r)+A[V T (r)-V A (r)}S 1 -S 2 
+ 4(VMr) + ^[3K(r)-Vi(r)] 



m 2 4 2r 



+ T 



+ 



Si • S2 



2V 2 VV(r) + V 2 Vp(r) 



(6.111) 



where we have ignored momentum dependent terms as well as possible m~ 2 corrections 
to V A {r) and Vr(r). 

Since QCD interactions are spin-independent to leading order, Va{t) = Vr(r) must 
be satisfied. Moreover, in comparing the above formula with the potential of Eqs. (|6.38|) 
with tree level matching constants, q = 1, di = 0, one finds, 



(Mi 



Vo 

V-2 



V v + V s , 

-a-v)v s , 

V v + vVs, 

Vy-V'p 



2VVy + V 2 Vp, 



(6.112) 
(6.113) 
(6.114) 

(6.115) 
(6.116) 



where we have also used the relation, Eq. ( |6.48 ). There are indications that the linear 
term, err, within V should be purely scalar since vector type potentials are thought to 
rise at most logarithmically in r ||406|| . The Darwin term appearing within Eq. (|6.39|) 
implies that 6 = 0, i.e. any scalar contribution to V 2 Vo has to be cancelled by the 
combination, V 2 Vq E + V^Vj 3 . It is clear that the picture becomes more involved when 
the matching coefficients assume non-trivial values. Moreover, the assumption that the 
interaction kernel only depends on the momentum transfer does not necessarily apply. 



6.5 Beyond the adiabatic approximation 

We shall briefly discuss the interrelation between local potentials, sum rules and the 
stochastic vacuum model. Following this, we shall describe pNRQCD which is a sys- 
tematic and conceptionally attractive approach to quarkonia bound state problems. 
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Subsequently, we will discuss some consequences that arise from including MD poten- 
tials. Finally, we incorporate hybrid states and transitions between different gluonic 
excitations of the string into the potential approach. 



6.5.1 Are potentials enough? 



The local potential picture of heavy quark bound states has often been challenged. 
Voloshin [|407|1 and Leutwyler [|408|1 for instance investigated the effect of the gluon 
condensate on quarkonia levels and found a dependence proportional to n G {a s F 2 ), on 
the principal quantum number, n. From this they concluded that this effect could not be 
reproduced by a local potential. However, a term growing that rapidly would certainly 
dominate the spectrum, if not for n = 2, then for n = 3, in contradiction to experiment. 
In this light, it appears questionable whether all non-perturbative physics relevant for 
excited state quarkonia can be approximated by the gluon condensate alone or if other 
infra-red scales play a role. The gluon condensate does not result in a linear contribution 
to the static potential but will only add a short distance term, proportional ||409|| to 
r 2 , to the perturbative result. Thus, the gluon condensate alone is not sufficient for an 



understanding of non-perturbative physics at large (as well as small, cf. Section |4.8| ) 
distances. Based on somewhat different arguments this has also been pointed out in 
Refs. HTOl, Bill. 



One instructive extension of the sum rule approach is the stochastic vacuum model 
(SVM) by Dosch and Simonov [|R Btl in which non-local condensates, i.e. correlators 
of field strength tensors at different space-time points, 



D{x) = (a s F 1 (x)U A (x)F 2 (0)), 



(6.117) 



are introduced. Fi symbolise linear combinations of electric or magnetic fields. Calcu- 
lating a Wilson loop in this approach indeed yields a linear contribution to the static 
potential at large distances pi] . In order to achieve gauge invariance of the correla- 
tion function, the adjoint Schwinger line, U A , has been included into the definition, 
Eq. ( |6 . 1 1 7| ) . Note that the above non-local condensate resembles the gluelump corre- 
lator of Eq. ( [5.2| ). It is not entirely clear how to cancel the self-energy contribution 
that is due to the Schwinger line and how to interprete the possible path dependence 
of the result. Putting these problems aside for the moment, lattice determinations of 
such correlators by use of two different methods exist [ 412 , 413|, 414 1. 

The correlator will decay exponentially for large Euclidean separations, D(x) oc 
exp(— \x\/Tq) — > oo), with the gluon correlation time, Tq, being a second dimen- 
sionful infra-red scale. Let us further introduce the characteristic time scale, associated 
with a quark in a bound state, T ni oc A~* oc l/{mv ^). One can now distinguish between 



two limits. In the case, Tq 3> T n i, the non-local condensate can be well approximated 
by a local condensate. Therefore, the Leutwyler |[408|| result is reproduced and no local 
potential that describes the long distance behaviour can be found. This is not too sur- 
prising, though, as one would expect the adiabatic approximation to be violated if the 
characteristic time scale of the gluon dynamics becomes larger than that associated to 
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the heavy quarks. On the other hand, for gluons harder than the bound state energies, 
T G <C T n i, the effect of the non-local condensate cannot be neglected and under cer- 
tain additional model assumptions one indeed finds level splittings to scale like 415], 
AE n i oc Tc(r 2 ). This would imply the local potential itself to be proportional to r 2 at 
small distances, in contradiction to the lattice results but in agreement with sum rule 
expectations^ on the static potential ||409|| . However, one would not expect the SVM 
to reproduce the correct behaviour for distances, r < Tq, anyway. 



6.5.2 Potential NRQCD 

A more systematic approach to the bound state problem is potential NRQCD (pN- 
RQCD) |J7l|, [372], |33j|, the QCD generalisation of pNRQED |3j, |T§, |T7j in which 
on top of the NRQCD Lagrangian, an expansion in terms of the quark separation, 
r oc l/(mv), is performed. The remaining colour fields are living at the centre of 
mass coordinate, 0. By means of a multipole expansion, A^ijc^t) can be obtained from 
A^(0,t) and derivatives thereof. The resulting Lagrangian is p85| . 
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O f r • EO + O f Or • EH + 



(6.118) 



where V s (r), V a (r), Ve(t) and V#(r) represent (infinitely many) matching coefficients 
that have to be determined by some prescription. Apart from r the coefficients depend 
on the scale fi and, to higher orders of the expansion, spins and momenta. Since all 
r dependence has been separated from the interaction terms, these can be factorised 
according to their properties under local gauge transformations. V s (r) and V (r) can 
be identified with the singlet and octet potentials of Section [4.5| in the case that no 
relevant physical scale exists between mv and mv 2 ; S is the colour singlet contribution 
to the wave function while O represents the colour octet part. 

Interestingly, in the situation, A Aqcd, a non-perturbative r 2 contribution to V s 



is obtained, in agreement with Ref. |409|1 . For details we refer the reader to Ref. ||335 
We also remark that in Ref. ||379|| vNRQCD is introduced which is based on a similar 
multipole expansion in momentum, rather than in position space. 

While in pNRQCD local and non-local terms are clearly separated, unfortunately, 
it is not clear how to arrange for such a factorisation in lattice simulations. Moreover, 
once the matching coefficients, Vi(r), are determined, all remaining dynamics are ultra- 
soft, requiring lattice resolutions, a -1 , of order mv or smaller. This would result in 

37 The static potential differs from the interaction potential between mov ing quarks of finite mass. 
Sum rules predict the latter to be proportional to r 3 at short distances [ 408 . 
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intolerably large discretisation errors, unless one is interested in top quarks. However, 
the form of the pNRQCD Lagrangian with its transitions between singlet and octet 
states is quite instructive. 



6.5.3 Consequences of momentum dependence 

We will briefly discuss an effect that is sometimes mistaken as a violation of the adiabatic 
approximation: let us assume that the spectrum, E N , N = {nll 3 }, and Coulomb gauge 
wave functions, ^jv(r) = (r|-0jv), of a quarkonium bound state are known. In this 
case one might attempt to determine the interaction potential from the Schrodinger 
equation, 

H\tP n ) = E n \tP n ). (6.119) 
In the non-relativistic case, we have dr/dt = p/fiR- Therefore, 

[H,t] = -£ = -— p. (6.120) 
at ii R 

Let us consider a Hamiltonian of the form, 

H=^ + Hi. (6.121) 

From the canonical commutation relation, [p, r] = —i, and Eq. fl6.120| ), one can easily 
see that, 

[#*,r] = 0, (6.122) 

i.e. the interaction term, = Vo(r), is only a function of the distance and does 
not depend on the momentum. In this case, the potential can be obtained, wherever 

Mr) t 1 0. 

Vo(r) = E N - -^<^f>, (6.123) 

where we have assumed rotational symmetry. Note that Vq does not depend on the 
state \4>n) under consideration! 

To higher orders of the non-relativistic expansion not only spins and angular mo- 
mentum have to be included into the set of canonical coordinates but also Eq. ( 6.122| ) 
will in general be violated: the interaction Hamiltonian contains the explicitly momen- 
tum dependent terms of Eq. ( |6.41| ). Ignoring SD terms as well as the correction to the 
kinetic energy to keep the expressions simple, we have, 

H^Vo^ + Vu^p). (6.124) 

Naively applying Eq. (|6.123|) will result in the effective interaction potential (due to 
being forced to depend only on the position variables) to change with the state under 
consideration, 

V N (r) = V a (r) + + (6.125) 

{r\y>N) 
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Table 6.3: Combinations of spins and angular momenta that can couple to J = 1 



k 


s 


A 


S 


1 


y+ 


P 







P 







P 


1 




P 


1 




D 


1 




D 


1 




D 


1 





this dependence of Vn{t) on the state has nothing to do with the Lamb shift of QED 
since the (MD) potential, V(r, p) = Vo(r) + Vmd(^, p), of Eq. (|6.124j ) does of course not 
dependent on the quantum numbers N. 

6.5.4 What is the effect of hybrid states? 



From the discussion of Section |5]2| it is clear that gluonic excitations can play an impor- 
tant role in bound state problems. In general, the total angular momentum will be the 
sum of the angular momentum due to the relative movement of the quarks within the 
bound state, L = r A p, and the spin of the gluons, S 5 : K = L + S g . Af = S g r denotes 
the projection of the gluon spin onto the inter-quark axis and A = |Aj |. K 2 has eigenval- 
ues, k(k + 1), k > A. K will couple to the quark spin to give the total spin of the state, 
J = K + S. We also recall that the gluonic string could be classified with respect to, 
A^", where r] denotes the combined parity under charge inversion and reflection about 
the midpoint of the axis, and a v denotes the symmetry under reflection with respect to 
a plane, containing the axis. £ states with a v = ±1 fall into two different irreducible 
representations of the relevant symmetry group, D^, while irreducible representations, 
A > 1, which are two-dimensional, contain both a v parities: |A±) = (|A) ± |— A)), 
with a v \A±) = ±|A±). The resulting (hybrid-) quarkonium state has the symmetries, 

P = a , u (_l) fc + A + 1 ) C = a v r](-l) k+A+s . (6.126) 
In general, many possibilities exist to realise a given 

jPC 

assignment. In Table 16.3, 



we illustrate this by listing all combinations that yield a vector, J PC = 1 . Note that 
even without considering hybrids, the state can either be an S (k = 0) wave or a D 
(k = 2) wave. In a direct lattice NRQCD simulation of the spectrum, all the above 
combinations will share the same 1 ground state and none of the quantum numbers, 
s, k, A, are strictly conserved. However, we shall see that mixing between S and D waves 
for instance is likely to be small, such that almost pure S or D states, that can be created 
by different almost orthogonal operators, should still be distinguishable. In the potential 
approach mixing effects have been completely neglected so far and they may matter, 
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at least for high radial excitations. Dipole transitions are suppressed by order c -1 in 
the NRQCD velocity expansion while quadrupole transitions are accompanied by pre- 
factors, c -2 . Dipole induced mixing effects will be suppressed by order c -2 with respect 
to the leading order NRQCD Lagrangian and should, therefore, be included into an 
order c -2 spectrum calculation. k 3 will not be affected by magnetic dipole transitions, 
however, S3 and 77 are changed. Magnetic transitions also alter the representation: 
the 3 Si state in the Table can mix with hybrid X P\ states, which contain a flux tube in 
the S~ or in the Ii u representation. 

Electric dipole transitions cannot affect S3 or r] but change k%. As the Table reveals, 
only the mixing of 3 Si states with 3 Pi T, g states is possible in this case. We shall, 
however, see that the corresponding transition amplitude vanishes identically. Either 
a quadrupole transition or two separate dipole transitions connect the S and D wave 
E+ states. Therefore, mixing effects between these channels only have to be considered 
from order c -4 onwards. 



xi=(0,x) ^ X2=(r,x) 




r 



Figure 6.7: The four point function, G a b- 



In the derivation of the Schrodinger type bound state equation, Eq. ( |6.34| ), from 
Eq. (|6.32|) we have assumed that quarkonia can be completely classified by the quantum 
numbers of the constituent quarks. If this is not the case, the two-particle Green 
function, G, of Eq. ( |6.28| ) and Figure ^]5| has to be generalised to the G a b of Figure |6~7] . 
where the indices, a and b, run over all excitations that will contribute to the J PC of 
interest. To account for energy level shifts of S wave vector mesons, T(nS), to order 
c~ 2 , clearly only a,b = Il u are relevant. All other channels decouple to this 
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order in c 1 . The Hamiltonian, H, acts on G a b, 

-^-Gac = Y, H -bG bc , (6.127) 
at b 

and the resulting Schrodinger equation reads, 

Y,H a ^ h njP c{v) = E nJP cr njP c{r). (6.128) 

b 

The normalisation is such that, J^ai^Pa 1 4>a) = 1- Note that the state vector, (V^jpc), 
now contains information about gluonic excitations too. 

The 0(1) Hamiltonian is diagonal in the space of hybrid excitations and, to this 
order, K, S and A° v are separately conserved. However, to 0(c~ 2 ), off-diagonal elements 
appear and the direction of ip will change with time. To compute the off-diagonal 
elements of H we introduce W a b, Wilson loops where the spatial transporter at t = 
is in representation b of and at t = r in representation a. The orthogonality 
of states within different representations of implies, (W a b) = 8ab(W a ). We now 
intend to relate the generalised four point function, G a b, to the expectation value of 
W a b, with appropriate colour field insertions on the temporal lines. Let us first consider 
the corrections from fluctuations around the static propagator, x 2 (t) = r + v(t — r/2). 
The expectation value of the perturbed Wilson loop, W^ b , can be related to that of the 
static Wilson loop, 



(W: b ) = (W ab )+wg f T/2 dtt(E(r,t)W ab ) : (6.129) 

J-t/2 



the integral vanishes, unless the expectation value is negative under time reversal, CP = 
T = —1, in which case the correction matrix element itself disappears. We conclude 
that to the lowest non-trivial order, electrically mediated transitions between different 
hybrid potentials do not exist. 

Next, we consider magnetic transitions. The relevant perturbation term, Aif ss (x), 
is given in Eq. ( |6.78| ). In analogy to Eq. fl6.77Q , we obtain, 



(a\AH ss \b) 

Hah= {(a\a)(b\b)) l/r ( } 

We have introduced the denominator, such that \a) and \b) do not need to be normalised. 
In the equal mass case the matrix element can be expressed in terms of Wilson loops 
in the following way, where we have exploited the fact that B is even under parity 
inversions and, S = Si + S 2 , 

H ab (r) = SE^v^r), (6.131) 
m 

V ab ,i{r) = g lim((B i (0,r/2))) ab , (6.132) 
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where, 

//z?\\ (Tr(U ab F)) . qq\ 

mu = WSm^- ( ' 

U a b is a path ordered product of SU(N) matrices, starting from and ending at the 
space-time position of F, with Tr U a b = W ab . Note that V a b,i = H<M have dimension m 2 . 

If we are interested in corrections to states only, it is sufficient to consider the 
leading order diagonal elements, 

2 

H 0>a = 2(m - 5m) + + K(r), (6.134) 

where 14 denotes the respective hybrid potential. We can start from the unperturbed 
(diagonal) Hamiltonian, H , and determine the spectrum in all hybrid channels, 

H , a ^N a = E^ a i^ a . (6.135) 

Subsequently, the order c~ 2 corrections to the levels can be determined in pertur- 
bation theory, the corrections to the diagonal part, H 01] + in first (and, for spin-spin 
interactions as well as MD corrections, second) order, the corrections due to mixing 
with hybrids, AE™ X , in second order, 



AE™*= J2 



rpO,a t-,0, S 9 



(6.136) 



Note that radial excitations like 3S and AS whose energy levels are close to those of 
hybrid states, will be more strongly affected by the mixing than IS or 2S states, that 
are separated from the hybrids by substantial energy gaps. Also note that although the 
above equation very much resembles the general form of Eq. Q6.77 ), in Eq. (|6.136|) static 



hybrid state creation operators are substituted by wave functions of quarkonia bound 
states, and hybrid potentials by (r-independent) quarkonia energy levels. 

In QED similar mixing effects between |e + e~) states and |e + e~7) states exist ||416|, 



417| , [418|| . In QCD such effects are naively enhanced by factors, a s v 2 : /(af s v 2 + e -), with 



respect to QED, however, the denominator of Eq. (|6.136|) guarantees an additional 



suppression; the lowest hybrid level is well separated from the ground state and the 
spectrum of hybrid potentials is discrete, rather than continuous. In addition to tran- 
sitions between the ground state string and hybrid excitation, glueball creation can be 
considered. However, with masses of 3 - 4 GeV, the vector and axial-vector glueballs 
will only play a minor role while the scalar glueball will only enter the scenario at order 
c -4 , when quadrupole transitions have to be considered. In the case of QCD with sea 
quarks, additional flavour singlet meson channels open up, however, these particles are 
rather heavy too. Another possibility is the (OZI suppressed) radiation of three 7rs. 
The main change with respect to the quenched approximation is related to the spectra 
of static potentials at large r, where string breaking becomes possible. This will give 
rise to mixing effects with BB states. 
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It is interesting to observe that Eq. ( |6.128 ), which corresponds to the Lagrangian, 

+ ^k H *t>^ + ^-u H *t^ +•••)' (6 - 137) 

somewhat resembles the general form of the pNRQCD Lagrangian, Eq. (|6.118| ). In our 
case, ■0 E + replaces the singlet wave function, S, while the octet finds its analogue in 
various hybrids. An important difference is that, unlike in Eq. ( |6.118| ), the leading 
order mixing elements contain magnetic fields while electric contributions proportional 
to, r ■ E, have been found to vanish. Of course in higher orders of pNRQCD similar 
magnetic terms will appear too. 

The potential approach not only allows all sorts of effects to be systematically in- 
corporated but also enables the determination of many quantities that are not directly 
observable, for example the spectra of would-be hybrid states and the mixing matrix 
elements between these states and quark model states. This information is hidden in a 
direct lattice simulation. The results can readily be translated into languages commonly 
used in the context of the quark model and flux tube extensions thereof and put other- 
wise only heuristically defined concepts onto a firm basis. It also becomes obvious that 
the heavy quark interaction potential will only converge towards the static potential in 
the limit v/c — > 0, rather than m — ► oo as one naively might have assumed, ignoring 
the kinetic term in the NRQCD Lagrangian. However, unlike in heavy-light systems, 
v/c is not proportional to m but v/c oc a^r) [Eq. |S.1| : the desired limit v/c — > will 
be approached logarithmically slowly as the spatial extent r of the bound state wave 
function vanishes. This freezing of v/c as a function of the quark mass m at large m is 
also illustrated by the estimates in the last row of Table |6.1| . 

6.6 Lattice determinations of the potentials 



The potentials, Eqs. ( |6.56|) - ( |6.67|) , are given in a form in which they can be easily 



evaluated on the lattice. Spin dependent potentials have been computed in SU(2) gauge 
theory g0|, |T|, [420], |00], gOjJ , SU(3) gauge theory [|21], |22|, [63, |2(| [424], [k], |S 



and in exploratory studies of QCD with sea quarks ||425| , 426 ]. In Refs. [ 400 , [HJ the 
momentum dependent corrections in SU(2) and SU(3) gauge theories, respectively, have 
been considered too. The correction to the central potential, V,, of Eq. (|6.58|) [ [32]] as 
well as the transition potentials, V a b,i, of Eq. ( |6.132| ), however, have not been calculated 
so far. 

6.6.1 The method 

The simplest discretisations of magnetic and electric field insertions, gB and gE, are the 
clover leaf definitions of Eqs. ( [5.3|) - ( |5.6|) . Alternative discretisations have been inves- 
tigated in the first lattice study |[402|| of spin dependent potentials. Since the temporal 
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lattice extent is always finite, the limit, r — > oo, cannot be performed exactly. Moreover, 
the arguments of the integrals within Eqs. ( |6.56| ) - ( |6.67|) can only be obtained on a 
discrete set of t values. Spectral representations of the potentials, Eqs. ( |6.69| ) - ( |6.73D , 
however, are extremely useful to guide and control interpolations and extrapolations as 
well as in improving the lattice operators used. 



.At, 



At 2 
-♦ — 



• ♦ 



Figure 6.8: Lattice definition of the operator in the argument of the nominator of 
Eq. ( |6.68 ) for the example of F\ being an electric and F 2 being a magnetic field. 



Relative statistical errors explode with the temporal extent of a Wilson loop, T > r, 
while the distances between the field strength insertions and the spatial closures of 
the loop, At 1 and At2, determine the degree of pollution from excited states. There- 
fore, adapting the size of the Wilson loop within the double bracket expectation value, 
Eq. ( |6.68i ), to the distance between the two field insertions, t, T(t) = Ati + At 2 + t, 
turns out to be the optimal choice in terms of statistical errors as well as in terms of a 



fast convergence to the asymptotic limit of interest [|401|| . The resulting lattice operator 
is depicted in Figure |6l| In addition to keeping Ati large, the overlap with the ground 
state can be enhanced by smearing the spatial connections within the Wilson loop (cf. 
Section |4.7.1| ). In the first lattice studies | |421| the integrals, Eqs. ( |6.56 ) - ( |6.67[ ), were 
replaced by discrete lattice sums. By parameterising the arguments as continuous func- 
tions of t [ |330| , prior to the integration, discretisation errors can be reduced and 
the effects of the region of large t (where statistical errors dominate the signal) can 
still be incorporated. If the hybrid potentials are known, the exponents of such multi- 
exponential fits to Eqs. ( |6.69| ) - (|6.73|) can be determined independently ||330| , |419| . |420| 
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6.6.2 Matching to the continuum 



In all lattice studies, based on naive discretisations of the continuum expressions, the 
potentials V^', V3 and V4 have been found to be much smaller than one would have 
expected from perturbative arguments or quarkonia phenomenology. In Ref. [4~02] this 
has been attributed to the anomalous dimension of the magnetic moment while in 
Ref. ||421|| this has been interpreted as a lattice artefact. As we shall see, both suggestions 
are true in parts. In particular the difference, V 2 ' — V{, has been found to be a factor of 
three to four times smaller | |419| , |423| , [420 1 than the inter-quark force Vq, in violation of 



the Gromes relation, Eq. ( 6.48 



Nowadays, we know that such behaviour is caused by large renormalisations between 
lattice operators and their continuum counterparts p9[j . In addition, the matching co- 
efficients between NRQCD and QCD, discussed in Section |6.1| , will affect quarkonium 
spectrum predictions. We can separately perform two matchings: lattice NRQCD to 
continuum NRQCD and continuum NRQCD to QCD. In Ref. |[419|| a procedure remi- 



niscent of "tadpole improvement" ||120|| has been suggested to reduce the former renor- 
malisation factor: the lattice operators are improved by dividing out factors, Up, from 
the double bracket correlation functions. This prescription does not affect the contin- 
uum limit and still the leading order lattice artefacts are proportional to a 2 . However, 
in perturbation theory all lattice specific one loop self-interactions of the field inser- 
tions are cancelled. This procedure can be refined by the Huntley-Michael (HM) con- 
struction [[420|| , in which additional un- wanted higher order graphs cancel too. This 
becomes possible by taking the relative position of the field insertions with respect to 
the Wilson loop into account. This HM scheme has been employed in the simulations of 
Refs. gOg, [IJT], fZ}J, and as a result the Gromes and BBP relations Eqs. flB^gp - (|63D|) 
are found to be respected within the achieved numerical accuracy of a few per cent. 
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Figure 6.9: The potential V 2 V^ , together with a fit curve, V 2 V^ = — (2a + b)/r. 
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It has been suggested ||423|| to fix the lattice renormalisation factors non-pertur- 
batively from the Gromes relation at distances, r ^> a a , where rotational symmetry 
is effectively restored on the lattice. From Eqs. ( |6.69| ) - ( p773[ ) it is evident that the 



relativistic corrections to the static potential, rather than being spectral quantities 
themselves, are proportional to amplitudes, £)W(r), which will undergo renormalisation. 
Ratios of these amplitudes for different n, however, should approach the continuum 
ratios, up to order a 2 scaling violations. Let us define the renormalisation constants, 
Z B (r) and Z E (r), 

($ r , |B(0)|$ r , n > = Z B (r)($ ri0 |B L (0)|$ r , n ) (6.138) 
($ r ,o|E(0)|$ r , n ) = Z E (r)($ r , |E L (0)|$ r , n >, (6.139) 

where n should be chosen such that the corresponding amplitude does not vanish. From 
considerations analogous to Eq. (|6.74j) it is obvious that V[ and V 2 ', measured on the 
lattice, have to be multiplied by factors ZeZb, V3, V4 and V 2 V^ by factors Z 2 B and all 
other potentials by Z E to make contact with the potentials in a continuum scheme. 




Figure 6.10: The spin-orbit potential V{, with a fit, —V{ = a + h/r 2 . 

The definitions of Z E and Z E are ambiguous; any term that vanishes at least like a 2 
can be added. Since the left hand sides of Eqs. ( |6.138| ) and fl6.139|) are approached by 
the right hand sides in the continuum limit, Z E and Z E are r independent, up to order a 2 
lattice artefacts, and can be defined from the value at r = r , for instance. Therefore, 
we only have to distinguish between four independent renormalisation factors, Zb,±, 
Zb,\\i Ze,± and Ze,\\, where _L refers to a component orthogonal to the inter-quark axis 
and || parallel to the axis. By demanding the Gromes relation, Eq. ( |6.48| ), to hold 
for r ~ To, different linear combinations of products between Z E and Z& components 
can be determined. By varying the direction of r, the three combinations, Zb,\\Ze,±, 
Zb,±Ze,\\ an d Z B ,±Z E ,±, can be fixed. From the BBP relations, Eqs. ( p.49|) and ( |6.50|) , 
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all Z E Z E products can be over-determined. Therefore, a completely non-perturbative 
evaluation of the renormalisations required to restore the continuum Lorentz symmetry 
is viable. From the rotational symmetry of the relativistic correction potentials, ob- 



served in Ref. |p9|| , one can conclude Zb,\\ ~ Z B> ± as well as Z E u w Z Ei ±. In fact, up 
to the inherent order a 2 ambiguity, one would expect such (approximate) equalities if 
one considers that the renormalisation between lattice and continuum NRQCD is an 
ultra-violet effect and, therefore, should be primarily related to properties of the local 
field strength insertions themselves, rather than to their interaction with the ultra-soft 
background of bound state gluons. Moreover, using the same argument, on an isotropic 
lattice, a a = a T , one would expect, Z E rs Z b . 



6.6.3 Results 



We conclude this section by reviewing the lattice results obtained in the most concise 
and precise study so far p9fl . The SU(3) potentials have been determined by use of the 
quenched Wilson action on isotropic lattices at (3 = 6.0 and (3 = 6.2, that correspond to 
lattice spacings, a" 1 « 2.14 GeV and a" 1 ~ 2.94 GeV, respectively. In this reference, the 
HM renormalisation procedure |[420|| has been employed. Subsequently, the continuum 
Gromes and BBP relations were found to be satisfied within the statistical accuracy of 
the study. 




Figure 6.11: The spin-orbit potential V 2 ', in comparison with the expectation, Eq. (|6.84 ). 



In Figure |6.9| , we display the result for V 2 V r Q E , together with a fit of the form of 
Eq. ( |6.81| ). V 2 V^ was found to be consistent with a constant. In the Figure, we have 
subtracted the fitted self-energy constants, C^a, from the data points. The resulting 
cpCf'a ~ —1 seemed to cancel c 2 F C^a « 1 almost perfectly. At (3 = 6.0 and (3 = 6.2, 
values, b = (1.13±0.45)er and b = (1.83±0.61)<r, have been found, respectively. The sign 
of the difference, although not statistically significant, coincides with the expectation 
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that the matching coefficient Cd decreases with the lattice spacing [Figure |6.4|| . Note 
that a value, b ^ 0, is incompatible with Eq. (|6.111|) that results from the assumption 
that the form factors, Vi, within the interaction kernel, Eq. ( |6.11(J| ), only depend on the 
momentum transfer, q 2 . 

The observation, 3cd& + 2ct > 0, means that besides the 5 like Darwin term, another 
1/r like mass- (and, therefore, flavour-) dependent correction to the central potential, 
with a coefficient of approximate size, — 2<j/(4m 2 ), exists. However, this correction, 
together with an additional —a /(6m 2 ) term from the MD potentials, yields an increase 
in the effective Coulomb coefficient of the Cornell potential of less than 2.5 % for bot- 
tomonium. In the case of charmonium the situation is less clear: the uncertainty in 
can result in an increase of the effective Coulomb coefficient of anything from 8 % to 
18 %. The effective Coulomb coupling within the static potential will weaken at short 
distances as soon as one goes beyond the tree level inspired Cornell parametrisation. 
This is, however, not the case for the coefficient of the mass dependent corrections pro- 
portional to, V 2 (Vo — Vo iPer t) ~ 2a /r whose relative weight will, thus, increase at very 
short distances. Considering the discussion of the potential at very short distances in 



Section [4.8|, such contributions could turn out to be more important than one would 



have assumed, guided by the Cornell parametrisation alone. 
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Figure 6.12: The spin-spin potential, V3, in comparison with the expectation, Eq. ( |6.85| ). 

In Figure |6.10| > the long range spin-orbit potential V{ is displayed, together with a fit 
of the form Eq. ( |6.83D , where the string tension has been taken from a fit to the central 
potential. The values, h = 0.071 ± 0.013 and h = 0.065 ± 0.009, have been found at 
the two lattice spacings, respectively. Therefore, the dimensionless parameter h turns 
out to be somewhat bigger than one fifth of the Coulomb coefficient, e ~ 0.3. Since cp 
increases with decreasing a, we expect h to decrease slightly as a function of the lattice 
spacing. We can ask ourselves at what lattice spacing, a' « 1/V, we would expect h 
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to assume its (unmixed) value, h = 0. From Eq. fl6.52| ), we can derive the relation, 
cp(ijf) = [(e — h)/e) cp(/x) ~ 0.78 cp(/x): a decrease of Cp by more than 20 % is required 
which, as can be seen from Figure |6.2| , will correspond to a scale (much) smaller than 
1 GeV. 

In Figures |6.11| and |6.12| , we display K,' and V 3 , together with the model expectations 
of Eqs. (|6.84j) and ( |6.85| ). After having determined e from the static potential and h 
from V[ there are no free parameters in the function displayed. Excellent agreement 
between the data and the predictions is found. V{ does not contain any long distance 
contribution and therefore can be identified with the vector potential, Vy, within models 
that are based on an interaction kernel that only depends on the momentum transfer, 
i.e. f] = 0, within Eq. ( JUlg) . The fact that V 3 « V{jr - K," implies [Eq. (pT5|) 1 



Vp(r) « cr and, therefore, Vp ~ 0. 



cr 
> 




Figure 6.13: The spin-spin potential, V4, in lattice units at = 6.2, in comparison with 
the expectation, Eq. (|6.86|) , and a one loop perturbative improvement. 



In Figure '6.1'S , we display the potential V4, determined at (3 = 6.2 in lattice units. 
We decided not to plot the potential in physical units since the behaviour expected from 
Eq. (|6.86|) is a 6 function. Hence, the result will be cut-off and discretisation dependent. 



For the clover leaf definition of the magnetic fields, employed in the study, the lattice 
S function has been calculated (indicated as "tree level" in the plot). Indeed, the data 
are described well by this expectation. A one loop improved version brings the data 
even more in line with the expectation. 

The errors in the lattice determination of the relativistic correction potentials are 
much bigger than those on the static potential of Figure |4.2| . However, one should keep in 
mind that the effect of these terms on the spectrum is suppressed by factors of v 2 / c 2 with 
respect to the static potential; even an error as big as 10 % on a 10 % correction term is 
completely tolerable for most phenomenological purposes as the induced uncertainty of 
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Figure 6.14: The potential, Vd, together with the curve Vd(r) = — a/9r. 



a few MeV on the T spectrum will still be smaller than the effect of neglecting higher 
order relativistic or radiative corrections. In general, operators involving electric field 
insertions result in stronger statistical fluctuations than magnetic fields. Therefore, V3 
and V4 are the most precisely determined potentials, followed by V{ and V{ while V 2 V^ 
as well as the MD corrections are subject to big statistical uncertainties. 

In the case of the MD potentials, this is particularly disappointing as the expecta- 
tions, Eqs. ( |6.88D - ( |6.91| ), all contain a long range part and are all numerically small, 
in comparison with the other potentials. This means of course that these potentials are 
not of prime phenomeno logical interest. However, being dominated by non-perturbative 
effects, they are needed to discriminate between competing predictions arising from dif- 
ferent assumptions on the QCD vacuum |60| , ^1] . As an example of a MD potential, 



Vrf(r) is depicted in Figure |6.14j , together with the expectation, Eq. ( |6.90| ). The fitted 
self-energies Cd have been subtracted from the data sets. Note that the BBP relation, 
Eq. ( |6.49| ), has been confirmed to hold for the self-energies, 2Cb + ^Cd + K, e if = 0. 
Clearly, further improved numerical simulations are required to arrive at definite con- 
clusions about the functional form of the MD potentials. 



7 Application to the quarkonium spectrum 

After having determined the potentials, quarkonia spectra and wave functions can read- 
ily be predicted, within the limitations of the non-relativistic and adiabatic approxima- 
tions. Vice versa, quarkonium spectra can in principle be used as an input to fix 
parameters that have not yet been determined accurately, in particular the match- 
ing coefficients appearing within the effective action. The same values could then be 
taken in lattice NRQCD studies or HQET calculations of heavy-light bound states and 
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their decay matrix elements. In particular, the S and P state fine splittings react in a 
very sensitive way towards variations of these coefficients. Unfortunately, the r/& whose 
splitting with respect to T states would yield the cleanest information has not been dis- 
covered yet. Moreover, the fine structure as well as decay rates, that are proportional 
to the wave function (or in the case of P states, its derivative) at the origin, probe the 
heavy quark interaction at very small distances. Here we will restrict our discussion 
to spectrum determinations and estimations of the systematic errors inherent in order 
v 4 (or c~ 2 ) continuum and lattice NRQCD as well as uncertainties from neglecting sea 
quarks. 



7.1 Solving the Schrodinger equation 

Once the interaction potentials are determined, the Schrodinger equation, Eq. (|2.8|) , can 



or m 



be solved numerically on any personal computer, either on a discrete lattice 
the continuum [12, 0, [IT]]. In the latter case, one would start by integrating the radial 
equation, Eq. ( |A.3| ) or Eq. ( |A.11| ), for the Hamiltonian H of Eqs. ( |6.96| ) and 1U6| ) , 



H \nlh) = E° nl \nlh). 



(7.1) 



Subsequently, the 1/m 2 corrections, Eqs. ( |6.98|) - (|6.100|) , can be treated as perturba- 
tions, 

Enju = E° nl + ± 5>ZZ 3 |^(r, J, I, s, p)\nll 3 ). (7.2) 



rrr 



By use of the identities [[U], [30 



(nll 3 \f(r)p 2 \nll 3 ) 
An(nll 3 \5 3 {r)\nll 3 ) 



fh (nils f(f 

I fnlls (0)1 2 

ma 

7T 



l + -(nll 3 \r- 2 \nlh)) 



nib- 



(7.3) 
(7.4) 

..,1. 



all perturbations can readily be computed from expectation values, (r a ) , a 
Note that fh, V(r) and e are defined in Eqs. (|6.105|) - (|6.109| ). 

The Hamiltonian of Eqs. ( |6.96| ) - ( |6.100| ) originates from the parametrisations, 



Eqs. Q6.8CQ - (|6.91|) , of Vo(r) - V e (r) that are in qualitative agreement with the lat- 



tice data. These lattice inspired parametrisations for intermediate and large distances 
can of course be combined with perturbative short range expectations ||367| , |258| , |398|| . 
for example along the lines of Ref . PU| , for the purpose of a phenomenologically more 
accurate description of bottomonia states. This is certainly worthwhile doing as soon 
as lattice results on the transition matrix elements discussed in Section |6.5.4 become 
available and some of the NRQCD matching coefficients, in particular cf, have been 
determined in a non-perturbative way. Anticipating such future results we can, how- 
ever, use the available lattice data to estimate systematic uncertainties that are due to 
radiative and relativistic corrections as well as neglecting sea quarks. 
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Figure 7.1: The bottomonium spectrum from the lattice potentials 



Before discussing such effects, we reproduce the bottomonium spectrum obtained 
in Ref. |23] from the lattice potentials in Figure |7TT| . The displayed spectrum has been 
obtained from the parameter values, 



e 
h 
b 
C 



r, 



i 



o 



0.32, 
0.065, 

1.81a = 1.81(1.65 -e) r \ 
0, 

4.676 GeV, 
406 MeV, 



(7.5) 
(7.6) 
(7.7) 
(7.8) 
(7.9) 
(7.10) 



with one loop matching coefficients, cp, d v and cp. While e, h and b have been computed 
entirely on the lattice, the quark mass, rhb, and the scale, ro, have been determined from 
a fit to the experimental spectrum. 
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The Figure illustrates the precision to which experiment can at present be repro- 
duced, without recourse to phenomenological input other than that required to fix the 
quark mass and the scale. It is not a priori clear whether the average deviations of 
almost 20 MeV are dominantly caused by relativistic and radiative corrections or due 
to ultra-soft gluons that have not yet been incorporated into the potential approach. 
Uncertainties resulting from the statistical errors on the potentials as well as differ- 
ences between data sets obtained at (3 — 6.0 and (3 = 6.2 cannot be resolved on the 
scale of the plot. In addition to results based on the fit parameters, extracted from the 
quenched simulations (e = 0.32), results for a stronger Coulomb coupling, e = 0.40, 
are displayed. We intend to model the changes that one might expect when includ- 
ing sea quarks by this latter choice of e. It is amusing to notice that, when ignoring 
the mass dependence of the matching coefficients, q = 1, all ratios of splittings come 



out to be consistent |30] with those determined in direct lattice NRQCD simulations, 
indicating that higher order relativistic corrections as well as effects due to ultra-soft 
gluons do not play a prominent role, at least for the lowest few levels; all differences 



between published lattice NRQCD results and the spectrum of Figure |T7T| are entirely 
due to different prescriptions for assigning a physical scale to the lattice results and a 
different choice of the matching coefficients, q. In the potential case, overall agreement 
with experiment has been optimised while in lattice NRQCD usually the most precisely 
determined 2 3 Si — l 3 Si or 1 3 P — l 3 Si splittings are taken as the only input. 

7.2 Systematic uncertainties 

Having a Hamiltonian representation of the bound state problem at hand it is straight 
forward to investigate how the spectrum changes when the input parameters are varied. 
For instance, fine structure splittings are to first approximation proportional to the 
matching coefficient, cf.. Such effects are discussed in detail in Ref. [0O| . Here, we briefly 
summarise the main results and discuss the uncertainties common to the potential and 
the lattice NRQCD approaches. In addition, the effect of neglecting sea quarks is 
investigated and finite volume effects for lattice studies of T properties are estimated. 

Unfortunately, no precision results on the corrections to the static potential in 
QCD with sea quarks exist. However, the static potential has been determined ac- 
curately ||246| , [255| , |167| , [248| , |247|| for n/ = 2 (cf. Section |4.7.4| ), and an increase of the 



effective Coulomb coefficient e by 16 % to 22 % has been detected [|167|, [248|1 for quark 



masses, m u = m d > m s /3 (cf. Figure ( |4.5| ). 
The ratio, 

iZ= m 2 , Sl -m 1 , Sl > 
m^jp - m 1 3 Sl 

reacts in a very sensitive way towards quenching. The potentials yield R ~ 1.38 which 
(while in perfect agreement with quenched lattice NRQCD [|386|| ) disagrees with experi- 
ment, R ~ 1.28. The dependence of R on e, keeping the parameters r^ 1 and m fixed, is 
displayed in Figure |7]2|. Values, e « 0.4, appear necessary for the real world with three 
active sea quark flavours to reproduce this ratio. Keeping in mind that two sea quarks 
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DC 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 

Figure 7.2: R = Am 2 sis / ^ m iPis as a function of my 1 for various values of the Coulomb 
coupling, e. The upper curves correspond to the lowest order Hamiltonian, the lower 
ones incorporate relativistic corrections. 

resulted in the effective Coulomb strength to increase by about 20 %, such an increase 
by 30 % indeed appears to be very reasonable. Our suggested sizes of quenching effects 
will be based on this estimate. From the Figure it is also obvious that while order v 4 
(c~ 2 ) effects on this ratio are small around the bottom mass, relativistic corrections 
explode in an uncontrolled way towards the charm: while (v 2 )y ~ 0.1 for bottomonia, 
( v2 )j/ip ~ 0.4 is not exactly a small expansion parameter anymore. 

Relativistic 0(v 4 ) correction terms affect spin averaged 2S — IS splittings by about 
2.5 % and 11 % for bottomonium and charmonium, respectively; the corresponding 
numbers for the 1 3 P — IS* splittings are 4 % and 8 %. No experimental values for S 
T states are available since pseudo-scalar rji, mesons are not yet discovered. The T 



Table 7.1: Relativistic and radiative corrections to T splittings. 



splitting 


C(v 4 ) 


0(v 6 ) 


radiative 


2S-1S 


13 MeV 


2 MeV 


5 MeV 


2 3 S 1 — l 3 Si 


17 MeV 


2 MeV 


6 MeV 


l^P-lS 


24 MeV 


1.5 MeV 


2 MeV 


13P - l 3 Si 


12 MeV 


1.5 MeV 


5 MeV 
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Table 7.2: Relativistic and radiative corrections to J/ip splittings. 



splitting 


0{v A ) 


0(v 6 ) 


radiative 


2S-IS 


65 MeV 


25 MeV 


15 MeV 


1 3 P - IS 


35 MeV 


15 MeV 


15 MeV 



Table 7.3: Relativistic and radiative corrections to fine structure splittings. 



family 


0(v 6 ) 


radiative 


quenching 


T 


10 % 


25 % 


35 % 




40 % 


70 % 


30 % 



2 3 Si — l 3 Si and 1 3 P — l 3 Si splittings that are therefore at present of greater interest 
become reduced by another 1 % and 2.5 % due to spin-spin interactions when switching 
on relativistic corrections. In Table 771, we display estimates of the effect of even 
higher order relativistic correction terms on various T splittings as well as the size of 
error induced by ignoring the mass dependence of the matching coefficients between 
QCD and NRQCD when simulating the theory at lattice spacings, 1.5 GeV < a" 1 < 
3 GeV. In Table [T72"| , the corresponding results for J/ip states are displayed. To set the 
scale: Am^sis ~ 580 MeV, Amipis ~ 430 MeV. Note that the order v A corrections 
have been calculated while order v 6 and radiative corrections are estimates only. Within 
order v A NRQCD, radiative corrections to the matching coefficients, that are of size 
a s log(m/ fi)/m 2 , dominate over relativistic correction terms that are accompanied by 
factors 1/m 3 , at least for bottomonia. 

In Table [fTJ we summarise the estimates of the uncertainties of the fine structure 
splittings. Since we only have results from the lowest order at which the splittings can 
occur, the relative sizes of the relativistic corrections can only roughly be estimated to 
be of order v 2 . We did not try to assign quenching errors to individual spin averaged 
splittings. Only mass ratios, and not the overall scale, can be determined from the 
QCD Lagrangian. Therefore, assigning a quenching error to an individual mass is 
highly subjective since the result will depend on the experimental input quantity used 
to fix the lattice spacing. One finds different scale determinations to scatter by up to 
20 % within the quenched approximation which should be interpreted as the overall 
systematic uncertainty. In contradiction to this philosophy, estimates on quenching 
errors are given for the fine structure splittings. These are explicitly proportional to 
the Coulomb coupling. The quenching error estimates have to be interpreted as typical 
changes of the size of fine structure splittings with respect to spin averaged splittings. 
Including sea quarks will result in an increase of such ratios. The effect of radiative 
corrections goes in the same direction, this is obvious from the continuum two loop 
inspired estimate of Figure ET2l Besides quenching, the latter uncertainty again seems 
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to be the dominant source of error. 

Indeed, using tree level matching coefficients, one underestimates P wave fine struc- 
ture splittings for e = 0.40 by almost a factor two [30|, compared to experiment. How- 
ever, for the ratio, 



m Xb2 m X bi 

tlv 



i-FS 



(7.12) 

"Xbl m XbO 

from which the dominant radiative correction cancels, one obtains Rps ~ 0.56 which 
has to be compared to the experimental value, Rps ~ 0.66. By incorporating running 
coupling effects into the parametrisation of V3 and calculating order v 6 effects, it should 
be possible to further improve the agreement. 
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Figure 7.3: Difference between T levels and their infinite volume values as a function 
of the lattice extent, La. 



The potential approach not only offers an intuitive and transparent representation 
of quarkonia bound state properties in a continuum context but it can also guide lattice 
simulations. By numerically solving the Schrodinger equation on a three-dimensional 
torus for instance finite size effects can be estimated. This has been done in Ref. |30t| , 
and the main result is displayed in Figure [773 
volume limit for n 



While the approach to the infinite 



1 states is monotonous, this is not so for radial excitations. Some 
states, in particular the 3S, show a non-trivial behaviour that results in infinite volume 
extrapolations from data obtained at lattices with a a L a < 2 fm to become uncontrolled. 
The relevant symmetry group on a torus (as well as on a discrete lattice) is Oh, rather 
than 0(3). The five-dimensional continuum 0(3) D wave representation splits up into 
the two-dimensional Oh representation, E, and the three-dimensional representation, 
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T 2 . It is amusing to see that rotational symmetry is not only broken for finite lattice 
spacing but also at any finite volume, with 1T 2 approaching the continuum ID state 
from above and IE approaching it from below. If one aims at finite size effects below 
3 MeV, a lattice extent L a a rs 1.5 fm seems to suffice for IS* and IP states while for 
25*, 2P and ID states, L a a 2 fm is required and for 35" or 3P even 2.5 fm become 
necessary. On a 1.5 fm lattice for instance one would underestimate the 3S level by 
more than 50 MeV. 

8 Conclusions 

QCD contains a rich spectrum of purely gluonic excitations. Glueballs and torelons 
that are colour singlet states can be realised, as well as glueballinos that transform 
according to the adjoint representation of the gauge group. Chromo electro-magnetic 
flux states between static colour sources in the fundamental or in higher representations 
of the gauge group can be constructed that have non-local gauge transformation prop- 
erties. Besides mesonic potentials and hybrid excitations thereof, baryonic three-body 
potentials and even more complicated situations can be investigated. 

All these excitations can be accessed in lattice simulations with much more ease than 
properties of states containing fermionic constituents. While the lattice reveals many 
interesting and non-trivial aspects of QCD in some cases it is hard to detect effects 
that very obviously do exist like the breaking of the hadronic string. Lattice results 
are extremely useful to test and improve models of low energy QCD. Moreover, phe- 
nomena like the Casimir scaling found between potentials between charges in different 
representations of the gauge group or the A law of baryonic potentials, provide some 
insight into hidden aspects of the dynamics of the theory. Such results give reason for 
optimism that once out of the chaos there might arise understanding. On the other 
hand QCD exhibits a complex vacuum structure, and even in the allegedly perturbative 
short distance domain non-perturbative effects seem to play a role in some cases. 

QCD predicts the quark model as a classification scheme of hadronic states to be 
incomplete. Nonetheless, it seems to do quite well; in particular those gluonic excitations 
that would make a difference, come out to be quite heavy. The quark model can of course 
be improved by incorporating the known gluonic excitations. In doing so, one would 
expect the lightest quark-gluon spin exotica to be vectors, J PC = 1~ + . This result 
can be systematically derived for heavy quark bound states. However, direct lattice 
simulations show that it also applies to the light meson sector. 

We have demonstrated that in a non-relativistic situation, it is possible to factorise 
gluonic effects from the slower dynamics of the quarks. This adiabatic approximation 
is violated when ultra-soft gluons are radiated, i.e. when the nature of the bound state 
changes during the interaction time. However, such effects can be incorporated into 
the potential formulation by enlarging the basis of states onto which the Hamiltonian 
acts. Moreover, the validity of the adiabatic approximation is tied to that of the non- 
relativistic approximation in so far as the transition matrix elements are suppressed by 
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powers of the velocity, v. 

Within the adiabatic framework, valence gluons that accompany the quarks in the 
form of hybrid excitations of the flux tube and sea gluons, whose average effect is 
parameterised in terms of interaction potentials, can be distinguished from each other. 
To lowest order of the relativistic expansion pure quark model quarkonia and quark- 
gluon hybrids exist, which then undergo mixing with each other as higher orders of the 
relativistic expansion are incorporated. Hybrid mesons become a well defined concept 
in the potential approach and translation into the variables used for instance in flux 
tube models is straight forward. 

It has been shown that potential models can be systematically derived from QCD. 
An understanding of effective field theory methods turned out to be essential for this 
step. The resulting Hamiltonian representation of the bound state problem in terms of 
functions of canonical variables offers a very intuitive and transparent representation 
of quarkonium physics. It highlights parallels as well as differences to well understood 
atomic physics. 

In view of phenomenological applications, a non-perturbative determination of the 
matching coefficients between QCD and lattice NRQCD is urgent. As an alternative 
to lattice NRQCD and the lattice potential approach, quarkonia properties can also be 
calculated from relativistic quarks by introducing an anisotropy, a T <C a a . However, 
the potential approach is unique in its capability to access high radial excitations and 
to determine wave functions. From a non-perturbative determination of the matching 
coefficients, simulations of heavy-light systems would benefit too. Another challenge 
is to generalise the results presented to heavy-light systems, i.e. to achieve a similar 
factorisation into sea (gluon and quark) effects and valence (quark and gluon) effects, 
within an expansion in terms of the inverse heavy quark mass, vrT x . 
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A The radial Schrodinger equation 

For a rotationally symmetric potential, the standard substitution, 

^ 3 (x) = ^r„ 3 (M), 



(A.l) 



into the Schrodinger equation, Eq. ( fT5| 

V 

results in the radial equation, 



+ V(r) 



nils \ 



nllz ' 



(A.2) 



v<ni(r) + 2p 



R 



E, 



nl 



V(r) 



1(1 + 1) 
2fi R r 2 



u n i(r) = 0, 



(A.3) 



with u(0) = 0, it'(0) = ^(O). In order to understand the dependence of the spectrum 
on the underlying potential, we discuss power law parametrisations, 



V{r) = \r l 



■1,1,2,.... 



The virial theorem implies, 



(T) = E — (V) 



2 + v 



E. 



(A.4) 



(A.5) 



where T = p 2 / (2/i#) denotes the kinetic energy and (T) = (ip n u 3 \T\ip nU3 ). For simplicity 
we have omitted the quantum numbers from Eq. (|A.5| ). The average relative velocity 
of the heavy quarks within the bound state can easily be determined from Eq. ( |A.5| ), 



V 



fi R 2 + u 



-E, 



whereas 



(O = \(V) = \^~E. 
A A 2 + v 



(A.6) 



(A.7) 



Eq. (|A.3| ) can be reformulated in terms of dimensionless variables by a simple scale 
transformation, 



p = (2n R \\\)*+» r, 
e = 2fi R (2fi R \X\y^ E, 
w(p) = u(r). 



As a result, Eq. ( |A.3| ) reads, 



<i + 



e n i - sign(A)p — 



P 1 



W nt = 0. 



(A.8) 
(A.9) 
(A.10) 



(A.11) 
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The primes now represent derivatives with respect to the argument p. The dependence 
of an energy splitting AE and a length scale I on the coupling strength A and reduced 
mass [Mr is evident from Eqs. ( |A.9| ) and ( |A.8| ), respectively, 

AEoc |A|^/V^, I oc (\\\fi B y^ . (A.12) 

For negative powers z/, level spacings decrease with increasing quark mass, while for 
positive exponents, the opposite is the case. 

A logarithmic potential, V(r) = C ln(r/ro), constitutes the limiting case between 
positive and negative v. Indeed, for such a parametrisation one obtains a velocity (v 2 ) = 
C / jj,R as well as quark mass independent splittings, AE [224] while I oc tqC 1 ^ / jJ 1 ^ 2 - For 



a Coulomb potential [y = — 1) we obtain, AE oc X 2 fiR and I oc 1/(A/xr) while a linear 
potential (z/ = 1) yields AE oc A 2//3 //iR 3 and / oc I/^/irA) 1 / 3 . For a detailed discussion 
of the connection between spectrum and potential, we refer the reader to an excellent 
review article by Quigg and Rosner 16| . 



B Euclidean Field Theory 

We summarise the conventions and notations used in this article. We start by translating 
some Euclidean space-time objects into Minkowski space-time (superscript M) with 
metric r\ = diag(l, —1, —1, —1) for reference: 

di = df = d l > M , d A = -id^ 1 = -id°' M , (B.2) 

7, = i^ M = -il?\ 74 = 7°' M = 7o M (B.3) 

Ai = —A i,M = Af 1 , A 4 = -iA°> M = -iA™, (B.4) 

Fli = iF i0 > M = -i F M, F ij = F i i> M = F M, (B.5) 

B = B M , E = —iFi M . (B.6) 

The above conventions conform to the anti-commutation relations {7^,71/} = 25^ for 
the Dirac 7-matrices. A^ denotes the electro-magnetic four-potential, F^ v the Maxwell 
field strength tensor and E and B its components, the electric and magnetic fields. 
While in Minkowski notation Lorentz indices assume values fi = 0, 1, 2, 3, in Euclidean 
notation, they run from 1 to 4. 

We denote the generators of the SU(N) group by T a , with a — 1, ... , Na, Na = 
N 2 — 1. They fulfil the commutation relations, 

[T a ,T b }=if abc T c , (B.7) 

with f abc being real, totally antisymmetric structure constants. In SU(2) the generators 
can be represented in terms of Pauli matrices, T l = a l /2, in SU(3) by Gell-Mann 
matrices, T a = X a /2. The vector potential lives in the Lie algebra, 
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while local gauge transformations generate the Lie group, 

Q(x) = exp[iu a (x)T a ] G SU(N). (B.9) 

Under a gauge transformation, Q, the field transforms in the adjoint representation, 

Ap^A* = Q{A ll -id$$. (B.10) 

Fields that are in the fundamental representation of SU(N), e.g. Dirac spinors, q(x), 
transform like, 

q(x) — > Q(x)q(x), q(x) — > q(x)Q'(x). (B-ll) 
The Dirac fermionic Lagrangian in Euclidian space reads, 

C f = 5(7,^ + m)q, (B.12) 

with the covariant derivative, 

D ll = d ll + iA ll , (B.13) 
while the Euclidean Yang-Mills Lagrangian, 

Cym = ^tr {Fy^FfJ) = ^tr (F™F^> M ) = -C^ M , (B.14) 

can be constructed from the field strength tensor, 

F^ = -i[D^ D v ) = d^A v - d v A„ + i[Ap, A u ). (B.15) 



The relative minus sign within Eq. ( [B . 1 4j ) with respect to the Minkowski version implies 
that solutions of the classical equations of motion minimise the action, which is bounded 
from below. Therefore, quantum fluctuations are suppressed with respect to classical 
solutions by factors, e~ ss , within the path integral measure. 

The phase and normalisation of the field strength tensor above is chosen such that, 

gBi = -eijkFj k , gEi = Fj 4 , (B.16) 

correspond to the chromo-magnetic and electric fields, respectively, g = ^Ana s denotes 
the strong coupling "constant". Note that our definition of the electric field, E = — iE M , 
differs by a phase i from some text book conventions, resulting in, EfEf > 0. The gauge 
action expressed in terms of the colour fields reads, 

Sym = \J d 4 x {E?E? + BfBi). (B.17) 
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C The perturbative /3-function 

The one, two and three loop coefficients of the /3-function, Eq. (|3.22|) , 



(5(a s 



da* 



din fx 2 



(C.l) 



have been calculated in Refs. [p7] P8| , |429| , p0| , plj and |432| , p3| , respectively, in 
the modified minimal subtraction MS scheme ||434| , |435| , [436|| of dimensional regulari- 
sation [)437|| . The n/ flavour results for an SU(N) gauge group read, 



A) 



V 3 3 V ^~ 



4vr' 



V 3 



A/5 
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2857 
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N 
N 



1 



V 54 
11 \ 



16vr 2 ' 
187 1 



36 AN 2 



n .f 



(C.2) 
(C.3) 
(C.4) 



64vr 3 ' 



while the four-loop coefficient has been calculated in Ref. | }438j| . The latter refer- 
ence also contains the coefficients for all compact semi-simple Lie gauge groups. The 
conversion between MS scheme couplings and the bare lattice coupling for Wilson glu- 
onic and fermionic action is know to two loops ||391| , |392| , |393| , [394j , f439j , [440|| . The 
numerical pure gauge result for the f3 function coefficient (3% reads ||441| , [442j ], 



Pi 



„ 2143\ 1 

-366. 2iV 3 + 1433.8iV -. 

N J 64vr 3 



while for SU(3) with rif flavours of Wilson fermions one obtains ||440 

1 



ft: 



-6299.9 - 1067/2/ + 59.89n 



fj 64vr 3 ' 

Translating between one scheme and another is straight forward: from 
a'(ii) = a{n) + ci« 2 (/i) + c 2 a 3 (/i) + c 3 a 4 (/i) H , 



one obtains, 



Ao 
A' 



A) j — Al, 

A ~ Cifa + (C 2 - cl)P , 

fa - 2 Cl (3 2 + clfr + 2(c 3 - 3c lC2 + 2c 3 )/? , 

Ae ci/(2/3 )_ 



(C.5) 



(C.6) 



(C.7) 



(C.8) 
(C.9) 
(CIO) 
(C.ll) 
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D The centre symmetry 



On a torus, a global symmetry is associated with each compactified space-time 
direction, /z, besides the invariance of the action and path integral measure under local 
gauge transformations. Zn denotes the set of the N Nth. roots of unity and z 6 Zn C 
SU(N). This means, 

[z,U Xttt ] = 0. (D.l) 
Multiplying all links crossing a hypersurface perpendicular to the fi direction by a factor 

r T nz _ j z U XjV V x M = 0, v = n , v 

-» u x>v - j ^ otherwise ' 1 j 

leaves traces of closed loops of link variables with trivial winding number around the 
boundary in //-direction invariant: since every such loop crosses every hypersurface an 
even number of times, all factors z that are collected when crossing in the positive 
direction are cancelled by the z* factors collected from negative crossings. In particular, 
this argument applies to all pure gauge SU(N) actions which are linear combinations of 
traces of such loops. The fermionic part of the action, containing a covariant derivative, 
however, explicitly violates this Z^ invariance. 



D.l The Polyakov line and deconfinement 



The position of the hypersurface of Eq. ( p.2| ) can be moved by means of ordinary gauge 
transformations. Therefore, in infinite volume, it can be sent to infinity and the centre 
symmetry will be in no way different from an ordinary (large) gauge transformation. 
On the torus, however, the surface can still be moved around but not removed. The 
simplest object that is sensitive to the centre symmetry is the Polyakov line, 



PU) = Tr {T 



I raL T 

exp I i J dx^A^x) 




(D.3) 



a loop encircling the temporal boundary. T denotes time ordering of the argument. 
Obviously, under a centre transformations with respect to the 4-direction, P z = zP. 
In the pure gauge case, where the centre symmetry is a symmetry of the action and 
the path integral measure, this meansf^], (P) = on any finite spatial volume. In the 
infinite volume limit there is, however, the possibility of spontaneously breaking global 
symmetries and, indeed, similar to Ising and Potts spin models, at high temperatures, 
the Zn symmetry is broken. 

The expectation value of the Polyakov loop can be related to the free energy of an 
isolated static colour source R H, 



I CP) I 



x 



-0F g 



oo 



(D.4) 



38 Due to translational invariance, the expectation value of a Polyakov line does not depend on the 
position, x. 
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where the inverse temperature (3 = T _1 = aL T , that should not be confused with the 
inverse Yang-Mills coupling, is related to the temporal lattice extent. The vanishing 
expectation value of the Polyakov line observed in low temperature lattice simulations 
implies an infinite free energy of an isolated quark and, therefore, confinement. Vice 
versa, above a critical temperature > T c in SU(N) gauge theories, the expectation 
value will move into the direction of one of the Nth roots of unity, implying a finite 
free energy and the possibility to eventually find isolated quarks. The case of QCD 
with sea quarks is interesting in so far as the centre symmetry is explicitly broken: an 
isolated quark comes along with only a finite free energy penalty. However, the energy 
required to isolate a quark is still sufficiently high to create a quark anti-quark pair 
out of the vacuum. Therefore, despite the fact that chromo-electric strings between 
opposite charges can break, the theory is still effectively confining colour sources at zero 
temperature. 



D.2 Torelons 

It is not only worth considering Polyakov lines wrapping around the lattice in the 
temporal direction but also to discuss their analogue, which we will call the Wilson line, 
that encircles a spatial lattice direction i, 

1 / aL a \ 

^) = 7/£ Tr I n u ^)A . (d.5) 



We have already included the projection onto zero momentum into the definition by 
summing over all spatial points. Note that the above sum over the component Xi yields 
L a identical contributions. In principle, a projection onto any momentum orthogonal 
to the direction of the Wilson line is possible. From the correlation function, 

(Re Lj(t)Re Lj(0)) oc e~ mTt (t -> oo), (D.6) 

the mass my of a torelon ||212| , [443|| can be extracted, an excitation that only exists 
on the torus and that corresponds to a colour flux tube wrapping around a periodic 
boundary |[211 . 



While for small spatial extents, aL a , the centre symmetry of the classical Lagrangian 
with respect to spatial directions can be dynamically broken, analogous to the finite 
temperature case, for sufficiently large aL a 3> T~ , the centre symmetry implies forQ 

iV > 3, 

= (LityLtiO)) = (Re Lj(t)ReLj(0)) - (ImL 4 (t)ImL,(0)). (D.7) 

Note that the imaginary part, i(ReLj(£)ImLj(0) + ImLj(t)ReLj(0)), of the correlation 
function vanishes by charge invariance. From the above equality, it follows that, 

(ReLj(t)ReLj(O)) = (ImL f (t)Im £<(())) = ^(Li(t)L*(0)). (D.8) 



:l<> 



For SU(2), U = L* is real. 
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Moreover, centre symmetry yields, 

(Re Li(t)Re Lj(0)) = 5 ij (ReL i {t)ReL i (0)). (D.9) 

The above two equations imply that all correlation functions between linear combina- 
tions of imaginary or real parts of Wilson lines are proportional to each other. There- 
fore, all torelon states that correspond to one unit of flux are degenerate. In addition to 
torelons corresponding to one unit of flux, torelons wrapping several times around differ- 
ent boundary directions can be constructed and labelled according to n = (m, n 2 , 713), 
n\ > n 2 > n 3 , rii = 0, 1, . . . 

j n maxi where N/2 — 1 < n max <J N/2 since the centre 
symmetry implies that states with winding numbers N ±n are indistinguishable from 
n wrappings. As soon as fermions are included into the action, the symmetry is 
broken and torelons corresponding to different representations of the cubic group Oh®C 
will in general assume different masses. 

In the limit of large aL a the situation of a closed flux tube encircling a periodic 
boundary becomes indistinguishable from a flux tube with fixed ends, created between 
point-like charge and anti-charge at infinite separation. Therefore, in pure gauge theories 
the energy stored per unit length will become identical to the string tension, a, the 
infinite distance slope of the static potential: 

ttit — > <jaL a (aL a — > 00). (D.10) 



E Matching NRQCD to QCD 

In this Appendix, we display results on the matching coefficients between NRQCD, 



Eqs. - (§TJ), and QCD, calculated in the MS scheme 40 O g3§ g3§ 



cf 



c d 

dss 
d sv 

dns 



+ 



a s (m) 


70 

2/% | 


_a s (fi) _ 





13 

1 + -^a s (m) 

07T 



[a s (m) - a s (n)] 



a s (m 



a 8 (ji) 
4 aPM 



^ 308 , 

+ Ti7 <l 



<x s (n) 



!37Q 
120o 



9 mf — 

4 a 2 s (n) 



2 / . m 2 1 \ 2 / , m i 1 



9 rrii ~~ m 2 

1 a 2 8 (fJ>) 

2 ml — ml 



m x 

9 mim 2 m — , 
m 2 



2 1 m 2 , 1 2 ( 1 m l 

m, In 1 — I — m 2 In — 

V n 61 V n 



6 



••Note that the result for cd derived in Ref. [|28| turned out to be incorrect ]444 



(E.l) 

(E.2) 

(E.3) 
(E.4) 



130 



6 



■m\rri2 
1 a$(p) 



2 ml 



5 mi 
-m\mi In — 
3 m 2 



m 



. /i 3 / 



+ 3 



mi 



m 2 3 
ln_+ 2 



m G 



In 



mi 



1 - 



Qtt 

1207T ' 

720tt ' 



ln^I + ln^ 



3 % 



(E.5) 
(E.6) 



(E.7) 

(E.8) 
(E.9) 



where 70 = 6/(47r) and 71 = (68 — 52n//6) /(167r 2 ) are the first two coefficients of the 
quark mass anomalous dimension function. The above values for the dj's only apply to 
the non-equal mass case. In the equal mass case, one encounters additional contributions 
that are due to annihilation diagrams^] ||380|| : 
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In addition there are the non-annihilation contributions of Eqs. ( E.3|) - 

for mi = m.2, 



that yield 
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41 We ignore imaginary parts within d° s and d% 8 . Such contributions, however, appear in the matching 
calculation and are related to the fact that deep inelastic QCD cross sections cannot be obtained 
correctly within NRQCD. 
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Note that a renormalisation group improved result for d vv that agrees with the above 
equations has also been derived in Ref. |[28|| . 

Within the potentials of Eqs. (|Og ), (|QDQ , (ggg ) and ( gMj ), the matching coeffi- 
cients, 
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[-4( m b + 4( m i- m 2, /x)] , 



(E.18) 
(E.19) 
(E.20) 
(E.21) 



are required. For the equal mass case, we obtain from Eqs. (E.10|) - (|E.17|) and (|6.10|) 
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